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1. INTRODUCTION 


In this volume we describe the model for the major and minor 
neutral species in the ambient atmosphere and the ambient ionosphere 
[ROSCOE Model 1]. The overall model consists of 11 subroutines of 


which three are major subroutines: 


a. ATMOSU provides the major neutral species and the gen- 


eral properties of the ambient atmosphere, 
b. SPCMIN provides the minor neutral species, and 


c. IONOSU provides the ambient ionized species and the gen- 


eral properties of the ionosphere. 


For simplicity in presentation, we have adopted flexible defini- 
tions of which species are major and which are minor. It is anticipated 
that the meaning will always be clear to the reader in the context of the 


usage. 


The overall inputs, some intermediate outputs, and final outputs 
for Model 1 are given in Table 1. 


A flow diagram of the 11 subroutines, with their driver routine 
for development and test problem, is given in Fig. 1. A brief, simpli- 


fied description of the working of the 11 subroutines follows. 


The subroutine ATMOSU is initialized on a call to ATMOSU(1, 
120.) to set up needed parameters and to evaluate the solar-flux- 


dependent Fourier coefficients used in computing the time-dependent 


values of T (the variable controlling the temperature gradient at the 


lower boundary (120 km) of the high-altitude model) and T,, (the 


Table 1. Inputs, Intermediate Outputs, and Final Outputs 
for Major and Minor Neutral Species and 
Ionosphere for Ambient Conditions 
(ROSCOE Model 1) 


INPUTS 


Time: Year, month, day, and zone time 


Place: Geographic colatitude and longitude; altitude 


SOME INTERMEDIATE OUTPUTS 


Time: 


Solar-Cycle Property: 
Minor Species: 


FINAL OUTPUTS 
Neutral Species: 


Ionized Species (2 90 km): 
Atmospheric Properties: 


Ionospheric Properties: 


Universal time, Julian day number, 
local (apparent) time, index for day 
or night 


Solar flux at 10.7 cm 


Fit parameters for day and night density 
profiles 


N2Q; 02, O, Ar, He, COg, N, NO, NOg, 
Oo( As), O3, H9O 


e, oO", M* 


Pressure, density, density scale height, 
and (gas) temperature 


Electron (and Ng vibration) temperature, 
effective ion-pair production rate 
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exospheric temperature). In this call the values of the time (HL, hours), 
the 10.7-cm solar flux (SBAR), and the day-or-night parameter (IDORN) 
are determined by a series of calls from ATMOSU to five auxiliary sub- 
routines (ZTTOUT, JULIAN, SOLCYC, SOLORB, and SOLZEN) and are 
passed to ATMOSU through ATMOUP Common. After an initialization 
call from ATMOSU to SPCMIN(1, ZH), daytime and nighttime fit parame- 


ters are determined for O and CO, and an initialization call is made to 


IONOSU(1, ZH). During the initialization of SPCMIN, six calls to FITTER 
and seven (direct) calls to SOLVE are made to determine the fit coeffi- 
cients for the day and night profiles of the minor smecies NO, N, 
us 
O4( A)» Oz, HO, and NO,. 
The working of the above-mentioned five auxiliary routines is 


as follows: 


a. Subroutine ZTTOUT, receiving from TIME Common the 
input parameters year (IYRS), month (IMONS), day (IDAYS), and zone 
time (ZT) at east longitude PLON, returns to TIME Common the year, 


month, day, and mean or universal time (UT) at Greenwich. 


b. Subroutine JULIAN, called with the input parameters of 
year (IYRS), month (IMONS), and day (IDAYS), returns the Julian day 
number at the first of the year (YRFJ), the Julian date for vernal equinox 
(VEQJ), and the Julian day number on the day of interest (DAYJ). 


c. Subroutine SOLCYC, called with DAYJ, computes the aver- 
age 10.7-cm solar flux (SBAR), an input to ATMOSU through ATMOUP 


Common. 


d. Subroutine SOLORB, called with YRFJ, VEQJ, and DAYJ 
and receiving UT from TIME Common, computes the Greenwich appar- 
ent time GAT, placed in TIME Common, and returns the north latitude 
(SOLLAT) and east longitude (SOLLON) of the subsolar point. 


e. Subroutine SOLZEN, called with SOLLAT and SOLLON and 
receiving PLAT, PLON, and GAT from TIME Common, computes IDORN 
and HL, inputs to ATMOSU through ATMOUP Common. 


Subroutine FITTER, called from both ATMOSU and SPCMIN 
with values Y(I) of the dependent variable at NPTS values of the inde- 
pendent variable X(I), the degree NO of the polynomial used as the fitting 
function, an index IKIND denoting whether it is the dependent variable 
itself or its natural logarithm that is to be fitted, and an index ISIGN de- 
noting negative or positive exponents in the polynomial, returns the 


polynomial coefficients determined by the method of least squares. 


Subroutine SOLVE, called from Subroutines ATMOSU, SPCMIN, 
and FITTER with elements A(I,J) of a matrix of constant coefficients, 


returns the solutions of NO simultaneous linear algebraic equations. 


The three major subroutines are ready for use after they have 
been initialized. On subsequent calls to ATMOSU(2, ZH), with ZH the 
altitude in kilometers, ATMOSU uses ATMOUP Common to return the 
pressure (PP), the mass density (RHO), the temperature (TT), the num- 
ber densities of six species (SNI(I), I=1,6), and the density scale height 
(HRHO). On subsequent calls to SPCMIN(2, ZH), ATMOUP Common is 
used to return the number densities of the six minor species (SNI(I), 
I=7,8, 13-16). On subsequent calls to IONOSU(2, ZH), ATMOUP Common 
is used to return the number densities of the three charged species 
(SNI(I), I=9-11) and the electron (and Ny vibration) temperature (SNI(12)) 
and IONOUP Common is used to return these same quantities (with dif- 


ferent names) and the effective ion-production rate (QDEF). 
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2. AMBIENT ATMOSPHERE AND MAJOR NEUTRAL SPECIES 


The main subroutine for the ambient atmosphere and the major 
neutral species is ATMOSU. It is based on the subroutine ATMOS de- 
veloped by R. W. Lowen [Lo-73a]. For the convenience of the reader 
we have reproduced Lo-73a here as Section 2.1. Comments, revisions, 
and extensions to Lo-73a are given in Section 2.2. To facilitate finding 
the location in Lo-73a to which a comment or change in Section 2. 2 
applies, we have added an encircled letter (keyed to Section 2.2) at 
the appropriate location in the margin of the reproduced paper. The 
correspondence between symbols in Lo-73a (and in the revisions and 


extensions) and their Fortran names in ATMOSU is given in Section 2.3. 


See Fig. 2 for a simplified flow diagram of ATMOSU and Table 2 


for a summary of inputs and outputs for ATMOSU. 


2.1 AN AMBIENT ATMOSPHERE MODEL FOR ROSCOE 


(On pages 15-27) 
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Call auxiliary routines ZTTOUT, 
JULIAN, SOLCYC, SOLORB, and 
SOLZEN to get HL, SBAR, and 
IDORN. Use day profile of [O] in 
SPCMIN and call FITTER to com- 
pute least-squares polynomial fit 
coefficients for fpay- Compute 
Fourier coefficients for tempera- 
ture gradient variable and the 
exoatmospheric temperature, 
parabolic coefficients for density 
scale height in 110- to 120-km 
transition region, coefficients for 
O(night) profile, coefficients for 
COQ2 profile in 100- to 120-km 
transition region, N2 at 230-km 
altitude for N-profile initializa- 
tion, and atmospheric properties 
at 90-km for IONOSU 
initialization. 


Subroutine 
ATMOSU 


Initialization 
Call 
? 


Return 


HIGH ALTITUDE MODEL 


Computes pressure, density, 
temperature, density scale 
height, and Ng, Og, O, Ar, He, 
and CO2 concentrations. 


Is 
Altitude, 
ZH, above 
120 km 


Return 


LOW ALTITUDE MODEL 


Computes pressure, density, 
temperature, density scale 
height, and Ng, Og, O, Ar, He, 
and CO9 concentrations. 


Return 


Fig. 2. Flow Diagram of ATMOSU Subroutine. 


11 


Table 2. Summary of ATMOSU Input/Output Variables. 


INPUT VARIABLES 
Argument List 
JJ 


ZH 


ATMOUP Common 


HL 
SBAR 
IDORN 


TIME Common 
IYRS 


IMONS 


IDAYS 
ZT 


PLAT 


PLON 


ALTODN Common 


NALTOD 


ALTKM(47) 


Calculation flag 


If ie = 1: calculate initialization parameters 
JJ = 2: calculate atmospheric properties. 


Altitude of interest (km). 


Local time (hrs). 
Average 10.7-cm solar flux fre"=* W/on" Hz) |. 


Parameter for day or night. If COSCHI is the 
cosine of the zenith angle of the sun at point P, 
IDORN is 1 for daytime, i.e., IF(COSCHI. GE. 
0.0), and is -1 for nighttime, i.e., 
IF(COSCHI. LT. 0. 0). 


Number of the year in the 1900's (e.g., 1974 
becomes 74) at east longitude PLON. 


Number of the month (e.g., February becomes 
2) at east longitude PLON. 


Day of the month at east longitude PLON. 


Zone time for the 15-degree longitude interval 
containing PLON (decimal hours). 


North latitude of point P (say, grid origin) 
(radians). 


East longitude of point P (say, grid origin). 


Number of altitudes at which the daytime 
O-values are specified as data in SPCMIN. 


The array of altitudes at which minor species 
are specified as data in SPCMIN. 


Se See eee 


ODA Y (27) 
ONITE(18) 


CO2(25) 


OUTPUT VARIABLES 


ATMOUP Common 


PP 
RHO 
rr 
SNI(1) 
SNI(2) 
SNI(3) 
SNI(4) 
SNI(5) 
SNI(6) 
HRHO 
FEHSEQ 


ALTODN Common 


S1Z2N 


Table 2. (Continued). 


The daytime O-values specified as data in 
SPCMIN. 


The nighttime O-values specified as data in 
SPCMIN. 


The COo-values specified as data in SPCMIN. 


Pressure (dynes/cm?) 

Density (e/em®*) 

Temperature (°K) 

No concentration A /fem”) 

O, concentration (1/em*) 
) 


Ar concentration (1/em”) 


O concentration (1/em 


He concentration jem”) 
CO, concentration (i/em”) 
Density scale height (km) 


Fractional error in hydrostatic equilibrium 


Ng density at 230-km altitude for use in 
N-density initialization in SPCMIN 
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Lo-73a 


AN AMBIENT ATMOSPHERE MODEL FOR ROSCOE 


R. W. Lowen 
Science Applications, Inc., La Jolla, California 92037 


ABSTRACT 


An initial computerized model of the ambient 
atmosphere has been developed for the ROSCOE code. The 
model reproduces within a few per cent all the data of the 
CIRA 65 atmospheres, extended smoothly to sea level in 
conformity with the U. S. Standard Atmosphere. The routine . 
is fast and requires only 120 cards, yet has a sound physical () 
basis, so that many continuous derivatives are provided and 
so that there is reason to believe the technique can be ex- 
tended to newer atmospheric data as they become available. 


i. INTRODUCTION 


Science Applications, Inc. (SAI), has the responsibility for the devel- 
opment of the next-generation radar and optical systems code ROSCOE, 
with the General Research Corporation handling the systems aspects, under 
sponsorship of the Defense Nuclear Agency. It is intended that this devel- 
opment shall draw upon the phenomenology experience of the entire com- 
munity. SAI must therefore take every opportunity to keep the community 
familiar with the status of ROSCOE development, and to solicit contributions, 
comments, and criticisms. Here we report on the status of one of the first 
phenomenology modules to be available in preliminary form, namely a model 


of the ambient atmosphere. 


2. REQUIREMENTS OF THE MODEL 


There are certain general requirements placed on all models for 
ROSCOE, such as high speed, low storage need, high accuracy, generality, 


high physics content, and smoothness. It is not always possible to achieve 


* . . . . . . 
See Section 2,2 (p. 29ff) containing notes describing revisions. 
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all of these, but they are goals. Then there are specific requirements of 
the ambient atmosphere model. The ambient atmosphere routine is required 
to provide certain properties of the atmosphere when entered with values of 
altitude (2 0), latitude, local time, day number, and year (or equivalent 
parameter describing phase in the solar cycle). The desired outputs are 
species densities and temperatures, and various quantities that are deter- 
mined by these, such as mass density, pressure, scale height, mean molec- 
ular weight, and so on. Outputs that are specifically excluded because they 
will be considered separately are the ambient ionosphere and ambient winds. 


Also excluded is tropospheric weather. 


As for the general requirements on speed, accuracy, and storage, 
the first is not likely to be severe; it would be difficult to write an atmos- 
phere routine that contributed seriously to overall ROSCOE running time, 
although called frequently. The accuracy requirement is first of all to re- 
produce atmospheric properties within their experimental uncertainties. 
This is not difficult, and turns out not to be the determining factor; rather, 
if the code purports to reproduce some atmosphere model in general use, it 
must do so with sufficient fidelity that the discrepancies produce no Signifi- 
cant differences from other codings of the same model. We have interpreted 
this to mean that discrepancies should be no more than a few per cent in any 


quantity that is important. 


The storage requirement, however, is a serious one, and all but 
eliminates certain kinds of models from consideration. Because atmos- 
pheric properties depend not only upon altitude but also (in various altitude 
regimes) on time of day, season, latitude, and solar cycle, and because so 
many output quantities are wanted, simple interpolation from tabular data 
would require very extensive tables (of course, not all would have to be in 


the machine at any one time). 


The requirement of generality means primarily that we want to choose 
an adequate form for the preliminary model, although the details may change 
as newer data become available. This can be achieved provided the model 
has a sufficiently high physics content. Finally, the requirement of smooth- 
ness means that it is desirable for the output quantities to possess as many 


continuous derivatives as possible. 


3. THE DATA BASE 


It is expected that new CIRA (COSPAR International Reference At- () 
mosphere) atmospheres will be issued shortly. Champion is preparing the 
high-altitude models, using the ideas of Jacchia's 3} static diffusion models. 


(3) 


Groves*’ is preparing the lower-altitude models. When these data become 
available, they will presumably represent a sufficient improvement over the 
data presently available that one will want to use them. At present, however, 
one can only proceed with what is in hand now, and try to select a model 


structure that can be readily adjusted to accept the new data. 


The older atmospheric data come from the CIRA 65 models and the 
U. S. Standard Atmosphere”) with the later models ‘®) of latitudinal and 
seasonal variations. Data on minor species are not provided in these sources, 
and generally have to be collected from the literature and collated; recent 


(7) 


and provide many references. 


reviews by Strobel’ © and by Shimazaki and ee have been very helpful 


Besides these data sets, the literature contains some examples of 
attempts at simplified modeling of high-altitude atmospheres. Deriving from 


an observation by Bates”) 


that the equations describing an atmosphere in 

diffusive separation can be solved exactly analytically for a certain form of 
temperature profile, these models have been fitted to satellite data by Stein 
and Walker 9), to Jacchia atmospheres by Walker!) and to CIRA 65 by 


Nisbet(22), 
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4, THE PRELIMINARY MODEL 


A structure for the Ambient Atmosphere Model has been chosen that 
(4, 5) 

data 
well, and that seems likely to be able to accommodate newer data as they 


rests on a sound physical basis, that fits the currently-available 


come along. It is convenient to divide the discussion into three sub-topics, 
Major Species at Low Altitude, Major Species at High Altitude, and Minor 
Species. 


4.1 Major Species at Low Altitude (0 < z < 120 km) 


For the preliminary ambient atmosphere model we have elected to 
follow CIRA 65 (4) in dividing the altitude regime into a low-altitude regime, 
O< zs 120 km, wherein the major and inert species are thoroughly mixed so 
that fractional concentrations are (almost) altitude-independent, and a high- 
altitude regime where diffusive separation prevails (120 kms z). Under 
these circumstances, the pertinent physics at low altitude comprises 
(a) hydrostatic equilibrium, (b) the perfect gas law, (c) the law of partial 
pressures, and (d) perfect mixing. Assumption (d) begins to fail above about 
80 km altitude, because of solar dissociation of Oo, and somewhat higher be- 
cause of increasing diffusion. The model can still be preserved by using as 


the defining quantity ''molecular-scale temperature," 


Tu = M,T/M , (1) 
where 
T = the true kinetic temperature mK), 
M = the mean molecular weight, 


M, = the mean molecular weight at sea level = 28.96 g/mole. 


This quantity is specified in both the CIRA 65 and U. S. Standard Atmospheres 


as a piecewise-linear profile, which permits the remaining equations to be 


integrated analytically. There are two difficulties with this procedure. First, 
it does not satisfy our preference for an atmospheric model with continuous 
derivatives. Second, CIRA 65 is not defined below 30 km altitude; U. S. 
Standard is, but differs, although the two are fairly close near 40 or 50 


kilometers. 


To get around these two difficulties we have arbitrarily selected a 
profile of the quantity g/ Tw where g is the gravitational acceleration, that 
agrees with Ref. (5) below 30 km altitude, agrees with Ref. (4) above 50 km 
altitude, and more or less agrees with both between 30 km and 50 km alti- 
tudes. The chosen profile is shown in Fig. 1, along with the profiles defin- 


ing the two atmospheres. 


We have next fitted this profile as a least-squares polynomial in 


terms of the altitude, z(km). An eleventh-degree polynomial, 


1 

z* Z (2) 
M_~ k=0 

fits to within a fraction of one per cent. The governing equations then yield 


the pressure, 


| 


5 Z 7} 
2 M 10 Z , i 
p (dynes/cm”) = p, exp}- ——R— peal =| (c) 


4) Zz 
M, __10 g(z) ga 
Py ©XP{- R Fé Ty) zy) 9% 


. 10° Sb Be ks ‘ 
ait | laces - aaa m1 * ; (3) 
k=0 
where 
p= 1.01325 x 10° dynes/cem” (sea-level pressure), 
R = 8.3143 x 10" erg/g-mole-°K (universal gas constant). 


] 
| 
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The Altitude Profile of Gravitational Acceleration 
Divided by Molecular-Scale Temperature Adopted 
to Define the Preliminary Ambient Atmosphere 
Model for z < 120 km. 


Because of biasing problems in this fitting procedure, the fit (3) does 
not provide an entirely adequate fit, especially at the important "join" alti- 


tude of 120 km. This can be fixed by multiplying the result (3) by a small 


correction, 
p(x) = pg(z) expl9.4144 x 1078 22-833) (4) (e) 
where 
Pz = the value given by Eq. (3). 
With the auxiliary equations 
g=¢,RL'R, +2 , (5) 


> FFE) ° 


n= @L/|() my «) 
n, = pL/M, , (8) 
ee (9) 
(1) 
[A] = 0.009n, , (10) 7) 
[O,] = 1.211n,-n, (11) 
fo] = 2n, (Me = 3) , (12) ( 
5 


[He] = 4.625«10°°n, , (13) (i ) 


this completes the low-altitude model, once M/M, is specified. Here the 


symbols are: 


g sea-level gravitational acceleration = 980. 621 cm/sec” 


oO 


Avogadro's number, 

the acceleration of gravity (em/sec") at altitude z, 
the radius of a spherical earth = 6367. 65 km, 
mass density (g/em”), 

kinetic temperature (°K), 

total number density fom"); 

total number density if no dissociation (em™), 


= number density of Ny (em™%), 


= number density of O, (em), 


number density of O fom), 
(i) [A] = number density of A (cm7°), 
[He] = number density of He (em™?), 


Note, please, that helium is included here among the ''major species,'' where 
it clearly does not belong, because (a) it is necessary to provide a value at 
the join altitude of 120 km, and (b) its physics in the low atmosphere is 


probably more like that of the major species than like that of the minor ones. 


The other problem facing the model has to do with M/M,. The CIRA 
(4) 


generated in the difficult region between 80 and 120 km altitudes, where 


65 description’ “ is too vague to permit one to understand how the model was 
solar dissociation of 05 begins and where the importance of diffusion in- 
creases. Results consistent with CIRA 65 are obtained if one follows a pro- 


cedure of (a) first calculating [O,] as at lower altitudes, and (b) then 


dissociating some of the [09], so that M/M, decreases. This gives the 


equations 


[0,] + 5[0] = 0.211 n, (14) 


[Nj] +[A]+[0,]+[O] =n, (15) 


from which Eqs. (11) and (12) follow. Our model is completed by a fit to 
M/M, of the form 


M/M, = 1/(1+f) , (16) 


where f is the quantity given by the eighth-degree polynomial fit 


| 

a k . | 

f= aes fi, Z : (17) (k) 
k=0 

the data actually fitted were values of M,[O]/2Lp, with [O] taken from Ref. (4) 

above 80 km, from Ref, (8) for 40 < z < 80 km, and artificially made to go to (1) 


zero at z=0. 


4.2 Major Species at High Altitude (120 km < z) 

For the preliminary high-altitude model of the ambient atmosphere 
we have elected to adopt the Stein and Walker model as adapted by (m) E 
Nisbet ‘12 to CIRA 65. Nisbet has given Fourier-series fits to the local-time i 
variation of exospheric temperature and to the temperature gradient at 120 km 
altitude, with coefficients that depend on the value of the 10.7 cm solar flux. (n) 
Using these, and with starting values at z = 120 km provided by the low- 
altitude model, one has the analytic solutions for species densities, | 
l+aty, iF 


ae (18) } 


n.(z) = n, (120) exp/-ry.t . 
i i vs Na-ae 75 
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SS 


where 


(2 - 120) R? 
= t= R, +a) (R, + 120 (geopotential altitude)* , (19) 
(p) ¥, = m,e/tkT_ , (20) 


thermal-diffusion coefficient 


R 
i 


-0.4 for He, 0.0 for all other species considered, 


a= (8+ Tiga. s (21) 


T., = exospheric temperature (°K) ; 


Tio9 = temperature at 120-km altitude Cx) , 


n, (120) = density of it) species at 120-km altitude (em >} : 
k = Boltzmann's constant (erg/°K) , 


=1 
F2(P,-Tya9) T/A), 156. + (22) 


th 


m, = mass of the i species (g/particle) 


From the results of (18)-(22), any other quantity of interest can be found, 


e.g., from the auxiliary equations 


p=) nm, . (23) 
i 


(q) *An incorrect equation for € is given in both Refs. (10) and (11). 


M = / > n, ete. (25) (4 


In the current model we have provided for No, Oo, O, A, and He. 
Because H has been neglected, the mass density given will be too low at ex- 
tremely high altitudes. It is believed that this is an unimportant defect for 


present purposes; H can readily be added if it matters. 


4.3 Minor Species 


Besides the major species that contribute most to the overall number 
density, there are a number of minor species that are important to IR emis- 
sion and transmission and other processes. These may include Oz, CO,, 
OH, NO, CO, NO,, NO, H,O, HO,, H, CH,, HNO, and others. The basic 
experimental data on the altitude distribution of these species are sketchy at 
best; theoretical understanding can hardly be any better. Their distribu- 
tions are believed to be governed by a combination of transport and chemical 


(8) 


for both their boundary conditions 


processes. The most elaborate calculations 


) 


consider only vertical trans- 
port, and have been sharply criticized’ 


and for their numerics. 


In view of these difficulties, it was decided to defer modeling of the 
minor species until development of the chemistry and IR models has pro- 
ceeded to the point where it is more clear just what is needed. Then a model 
will be provided using what data exist. Unfortunately, this model will un- 


doubtedly have a weaker basis in physics than the major species models. 
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5. SUMMARY AND CONCLUSIONS 


The preliminary ambient atmosphere model as described above has 
been programmed in FORTRAN for exercise and evaluation. The program 
has been used to calculate entire atmospheres for local times of 0, 4, 8, 12, 
16, 20 hours and solar 10.7-cm flux values of 75, 150, and 225 x 107°" watts/ 
m?” c/s (CIRA 65 Models 2, 5, and 8). Casual inspection shows that the re- 
sults are correct everywhere to within a few per cent, except that in regions 
where a species is unimportant it may be given less accurately. Unfor- 
tunately, a very thorough comparison with the data-base models would have 
to be automated, and would involve a vast amount of keypunching atmosphere 
tables, so this has been foregone. It is also difficult to illustrate compari- 
sons graphically when the quantities agree to a few per cent while varying 


over 17 orders of magnitude. 


No special care has been taken to optimize this version for speed, 
Since it is to be re-programmed at GRC for that purpose. Nevertheless, 
the current version gives results for a single altitude in about 1 msec on the 
CDC 7600. 


It is clear that the preliminary model is still incomplete; lacking are 
(a) latitudinal and seasonal variations, (b) minor species models, (c) multi- 
temperature models, and (d) consideration of excited states. Work on (a) 
has been deferred pending receipt of the newer CIRA atmosphere models, 
expected soon. Work on (b), (c), and (d) has been deferred until work on 
the chemistry and IR models has proceeded far enough to establish a better 
definition of model requirements. Meanwhile, the preliminary model has 
been delivered to permit early running of the overall ROSCOE Code. 


6. REFERENCES 
: L. G. Jacchia, Static Diffusion Models of the Upper Atmosphere 


with Empirical Temperature Profiles, Smithsonian Contributions to 


Astrophysics, Vol. 8, No. 9, 1965. 


10. 


a, 


12. 


L. G. Jacchia, Revised Static Models of the Thermosphere and Exo- 
sphere with Empirical Temperature Profiles, Smithsonian Astro- 
physical Observatory Special Report 332, 5 May 1971. 


G. V. Groves, Atmospheric Structure and Its Variations in the Region 


from 25 to 120 km, AFCRL-71-0410, 27 July 1971. ae pat 


CIRA 1965, COSPAR International Reference Atmosphere 1965, North- 
Holland Publishing Co., Amsterdam, 1965. 


U. S. Standard Atmosphere, 1962, Available from the U. S. Govern- 
ment Printing Office, Washington, D. C. 20402. 


U. S. Standard Atmosphere Supplements, 1966, Available from the 
U. S. Government Printing Office, Washington, D. C. 20402. 


D. F. Strobel, Minor Neutral Constituents in the Mesosphere and 
Lower Thermosphere, Radio Science 7, 1 (1972). 
T. Shimazaki and A. R. Laird, Seasonal Effects on Distributions of 


Minor Neutral Constituents in the Mesosphere and Lower Thermo- 
Sphere, Radio Science 7, 23 (1972). 


D. R. Bates, Some Problems Concerning the Terrestrial Atmosphere 
Above About the 100 km Level, Proc. Roy. Soc. London, A253, 451 
1959). 


J. A. Stein and J. C. G. Walker, Models of the Upper Atmosphere 
for a Wide Range of Boundary Conditions, J. Atmos. Sciences 22, 
965). 


J. C. G. Walker, Analytic Representation of Upper Atmosphere Den- 
sities Based on Jacchia's Static Diffusion Models, J. Atmos. 


Sciences 22, 462 (1965). 


J. S. Nisbet, On the Construction and Use of a Simple Ionospheric 
Model, Radio Science 6, 437 (1971). 


27 


| 
i 
| 
| 


4.2 COMMENTS, REVISIONS, AND EXTENSIONS TO LOWEN'S 
PAPER [Lo-73a] 


a. p. 15 The 120 cards in the original ATMOS (noted in the abstract) 
have increased to about 300 in Revision 09 of the Subroutine 
ATMOSU with an additional 300 substantive comment cards. 
The total number of cards in the deck (for ATMOSU, 
SPCMIN, IONOSU, their DRIVER, and their associated 


subroutines) is about 2318. 


b. p.17  CIRA-1972 [CI-72] has now been published. 


c. p. 19 Inthe first of Eqs. (3), M should appear as M(z’) in the 


integrand. The remaining equations are correct. 


. 20 The analytic fit to the profile in Fig. 1 is given by Eq. (2). 


. 
i) 
_ 


p 

e. p. 21 In Eq. (4), read p(z) instead of p(x). 
p In Eq. (10), read [Ar] instead of [A]. 
p 


. ol The coefficient 0.009 in Eq. (10) is close to the value, 
0.008615, required to reproduce the value given by CI-65 
for [Ar] at 120 km, and it is also close to the value, 
0.00934, given by US-62 for sea level composition. 


hb. p. 21 In Eq. (12), note that the factor |(M,/M) - 1] equals the 
quantity f appearing in Eqs. (16) and (17). 


i, pe 21 The coefficient 4.625 x 107° in Eq. (13) is required to re- 
produce the value given by CI-65 for [He] at 120 km, but 
it is quite different from the value, 5. 24 » 1078, given by 


US-62 for sea level composition. 


j. p. 22 Read Ar instead of A (twice). 


In Eq. (17), note that M,[O]/2Lp = f. 


The [O|-values for 40 < z(km) < 80 were taken from the 
noon, summer curve in Fig. 5 of Ref. 8. (Provision for 
day and night values of [O] is made in ATMOSU.) 


The Stein and Walker model assumes that the temperature 


profile is 


£ : e -T¢ 
f= T > &,= Tis?” 


ao 


(17a) 


Nisbet's Fourier-series fits to the local-time (HL) depend- 
ence of the exospheric temperature T, and to the tempera- 
ture gradient Tt at 120-km altitude, with coefficients that 
depend on the value of the 10.7-cm solar flux (F) associ- 
ated with the CIRA-1965 tables, are given by 


24 


(18a) 


(18b) 
Values of the coefficients are given in the computer program. 
p. 24 Equation (19) should read 


Zz (z - 120) (R, + 120) 


— gdz° = —————__—_- (19-Rev) 
8190 “120 — e 


geopotential altitude above 120 km. 


p. p. 24 Equation (20) should read 


te, B190/7 eT (20-Rev) 


q. p. 24 Delete footnote. 


r. Pp. 25 Equation (25) should read 


M = Lo/))n, (25-Rev) 
i 


Ss. Pp. 25 Add Section 4.4, Density Scale Height. 
By starting with the definition of density scale height, 


p(z)/(-de/dz) , (26) 


and using the analytic expressions given in the paper, one 


finds that for the low-altitude model 


5 16 | 
M,x10° ¢ 1 k-1 
BH = | eee g kz 
p R Ty g/Ty, d k 


- 2.66710 x 107” 2: 8383 _ __2 


(27) 


and for the high-altitude model 
2 n.(z) m; 
H = i 


p R. + 120-2 ce ; 
(aD n, (2) mp hi bit | 
2 


(28) 


To provide a continuous density scale height across the 
boundary between the low- and high-altitude models, the 


following parabolic transition function can be used, 


t. 


p. 25 


H(2) = alc. - 110) + bles 0) se, —410— a0) = 120 


(29) 


fo) 
il 


H (110) 


b = derivative of density scale height at 110-km altitude, 
taken to be given by [H_, (110. 5) - H,(109. S)i/1.0: 


a = [H, (120) - 10b - H,(110)|/(120 - 110)? . 


Add Section 4.5, An Evaluation of the Departure from 
Hydrostatic Equilibrium. 


The possible departure of the ROSCOE atmosphere 
from hydrostatic equilibrium is of interest, and we report 


our findings here. 
4.5.1 Low-Altitude Model in ATMOSU 
In the low-altitude portion of ATMOSU, the pressure, 


given by Eq. (3) of Lo-73a, is (initially) obtained by inte- 
grating the hydrostatic equation, 


dp a 
az pg. (30) 


However, owing to limitations in fitting ¢/T, a correction 


factor is used to obtain the final expression for the pressure, 


b 

A 

p(Z) = p3(Z) e ae (31) 
where Z is in kilometers and 
11 e F 

pg(Z) = p, expl- ¢ > m4 gt (32a) 

k=0 

32 


8 


A = 9.4144 x 107 (32b) 

b = 2.833 (32c) 
5 

c=M,107R . (32d) 


To numerically test the extent to which the ATMOSU 
atmosphere departs from hydrostatic equilibrium, we 
evaluated the quantity 


107° dp/az _ 


FEHSEQ = -— e i ., (33) 


which is the fractional amount by which the acceleration at 
a point due to the pressure gradient fails to balance the 
acceleration due to gravity. By differentiating Eqs. (31) 
and (32a) and using Eqs. (2) and (6) of Lo-73a, one can re- 
write Eq. (33) as 


tc Ane * & 


g/ Ty M, 


FEHSEQ = (34) 


The results of evaluating Eq. (34) are given in Section 4. 5.3. 


4.5.2 High-Altitude Models 


The high-altitude model in ATMOSU is closely re- 
lated to the CIRA-65 model [CI-65]. So, before examining 
the ATMOSU high-altitude model, we consider the CIRA-65 


model. 


4.5.2.1 CDIRA-65 Model. One can test the departure of 
the CIRA-65 atmosphere from hydrostatic equilibrium by 
first noting that the pressure scale height, a quantity tabu- 
lated in CIRA-65, is defined by the expression 


33 


= p 
Hy” TC ap7az) ° - 


so that, if hydrostatic equilibrium obtains, the relation 


Pp. .. 
at = pe (36) 
p 


should be satisfied. The quantity FEHSEQ in Eq. (31), 


FEHSEQ = q = =" 4 (cgs units) , (37) 
p 


has the same meaning as in Eq. (33). An evaluation of 
Eq. (37) for the CIRA-65 Model-5 8-hr atmosphere shows 


that FEHSEQ varies from -9. 96 x 1074 at 120-km altitude 


to -6. 53 x 107? at 800-km altitude. These errors are 


probably due to errors in numerical integration. 


4.5.2.2 High-Altitude Model in ATMOSU. The ATMOSU 
high-altitude model is the Stein- Walker model as adapted 
by Nisbet to CIRA-65. To evaluate Eq. (33), one needs an 
expression for dp/dZ, which one can show, after consider- 


able algebra, to be 


dp -7¢ 
—— ae n 
: &120 2 
mn (1+, +7,)ae77)) 
-— n.{|y. + : (38) 
a 1 ee | 


34 


soe agement ee 


+ th nn RE a CE A A AT a TR IRR . 


ete 


4.5.3 Numerical Results 


For the low-altitude model in ATMOSU, one finds 
that the pressure derivative is too small to maintain hydro- 
static equilibrium by a fractional amount that ranges from 
0 at 0 km to 0.020 at 120-km altitude. 


In the high-altitude model, for HL = 11.66 hrs and 
SBAR = 157.57, the pressure derivative is too small to 


maintain hydrostatic equilibrium by a fractional amount 
that increases from 90% 10°” at 120<ken altitude to 
1.0 x 107? at 320-km altitude and then decreases to 1. 8 


10°? at 800-km altitude. 


2.3 SYMBOLS IN LOWEN'S PAPER AND ITS EXTENSIONS AND 
THEIR FORTRAN NAMES IN ATMOSU 


Symbol Fortran Definition 
a SA | Eq. (21). 
[Ar] SNI(4) Argon number density. 
A; A(I) Fourier coefficient for T. 
Bi B(I) Fourier coefficient for T. 
C, C(I) Fourier coefficient for T_. 
f SF M,[O|/2Lp; see Eqs. (16) and (17). In 


SFDA F(ZH) ATMOSU, SFDAF(ZH) is used for the day- 
time profile of f. 


fy, DD(I) Coefficient in polynomial fit to f, Eq. (17). 
In ATMOSU, DD(I) is used for the daytime 
profile of f. 


F FF The 10. 7-cm solar flux associated with the 
SBAR CIRA~1965 tables. 


Symbol 


go 


Fortran 


GG 
GA F(ZH) 


GZ 
AA(I) 


GDTM 
GDTMAF(ZH) 


Z 
o(z’) ; 
g(2 ) a’ GDTMI 
hit )” GTMIAF(ZH) 


11 
ays kz 
1 


[He] 
HL 

Hy 

H, (110) 


? GKKZ 


GKKZAF(ZH) 


SNI(5) 
HL 
HRHO 


HRO110 


H, (110. 5) HR1105 


H, (109. 5) HR1095 


H_, (120) 


HRO120 


dH (110Vaz DDZ110 


k 


L 


m 


SK 
BIGA 


SMI(I) 


Definition 


Acceleration of gravity at altitude z. 


Sea-level gravitational acceleration. 
Coefficient in polynomial fit to g/Ty,- 


Ratio of g to Ty Eq. (2). 


Integral of g/T yy 


Derivative of g/Tyy- 


Helium number density. 
Local time, in hours. 
Density scale height. 


Density scale height at z = 110 km. 


Density scale height at z = 110.5 km. 
Density scale height at z = 109.5 km. 
Density scale height at z = 120 km. 
Derivative of density scale height at 110 km. 
Boltzmann's constant. 

Avogadro's number. 


Mass of species i. 


Mean molecular weight. 
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Symbol 


N 


Fortran 
BIGMS 


BMBMS 


SN 
SNS 
SNI(I) 


SNIZ(I) 


SNI(I) 
SNI(3) 
SNI(2) 
Pp 


PZ 


TIF 


Definition 
Mean molecular weight at sea level. 


Ratio of M to M, for day (noontime) 
conditions. 


Total number density. 
Total number density if no dissociation. 
Number density of species i at altitude z. 


Number density of species i at 120-km 
altitude. 


Molecular nitrogen number density. 
Atomic oxygen number density. 

Molecular oxygen number density. 
Pressure. 

Sea-level pressure. 

Basic factor in pressure, given by Eq. (3). 
Universal gas constant. 

Radius of a spherical earth. 

Fourier coefficient for T_. 


Kinetic temperature. 


Molecular-scale temperature, Eq. (1). 
Temperature at 120-km altitude. 
Exospheric temperature. 


Altitude. 
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Symbol 


Fortran 
ALP(I) 


GAM 
GAMT*SMI(I) 


ZZ 

PI 

RHO 
DRODZN 


TAU 


Definition 


Thermal-diffusion coefficient. 


Eq. (20). 


Geopotential altitude (above 120 km). 
3. 141592653590. 

Mass density. 

Spatial derivative of p. 


Variable controlling the temperature gradient 
at 120-km altitude, Eq. (22). 
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3. AUXILIARY SUBROUTINES FOR ATMOSU | 


The purpose of the five auxiliary subroutines ZTTOUT, 
JULIAN, SOLCYC, SOLORB, and SOLZEN is to convert inputs that are 
convenient for the user to the inputs required by ATMOSU, SPCMIN, 
and IONOSU. It is assumed the user will locate his coordinate system 
in space and time by stating the geographic north latitude and east 
longitude, the date, and zone time (based on 15-degree intervals of 
longitude) in a 24-hr system. These auxiliary routines determine the 
universal time, Julian day number, local (apparent) time, the solar 
zenith angle viewed from the origin, an index denoting day or night, 


and the 19. 7-cm solar flux. 


These subroutines (except ZTTOUT) had their origin in the 
AFWL WORRY code (where they were known as JULIAN, SOLCY, ORB, 
and ZSOL) and were revised when they were incorporated into the early- 
version ROSCOE code |LL-75|. These routines, to which ZTTOUT was 
added, were further revised and laden with comment cards under the 


contractual effort for the current ROSCOE code. 


3.1 SUBROUTINE ZTTOUT 


Subroutine ZTTOUT converts a Gregorian calendar date (speci- 
fied by stating the year in the 20th century ([YRS), the month (IMONS), 
and the day (IDAYS)) and zone time (ZT) at a given east longitude (PLON) 
to the Gregorian calendar date and mean (or universal) time (UT) at 


Greenwich. 


See Table 3 for a summary of inputs and outputs for ZTTOUT. 


Table 3. Summary of ZTTOUT Input/Output Variables. 


INPUT VARIABLES 


Argument List 
None 

TIME Common 
IYRS - 


IMONS_ - 


IDAYS - 
ZT - 


PLON - 


Number of the year in the 1900's (e.g., 1974 
becomes 74) at east longitude PLON 


Number of the month (e. g., February becomes 2) 
at east longitude PLON 


Day of the month at east longitude PLON 


Zone time for the 15-degree longitude interval 
containing PLON (decimal hours) 


East longitude of point P (radians) 


OUTPUT VARIABLES 


Argument List 
None 
TIME Common 

IYRS - 


IMONS - 
IDAYS - 


UT - 


A possibly revised value of the input parameter, 
corresponding to Greenwich 


A possibly revised value of the input parameter, 
corresponding to Greenwich 


A possibly revised value of the input parameter, 
corresponding to Greenwich 


Universal time corresponding to the zone time ZT 
(decimal hours) 
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3.2 SUBROUTINE JULIAN 


Subroutine JULIAN converts a Gregorian calendar date (speci- 
fied by stating the year in the 20th century (IYRS), the month (IMONS), 
and the day (IDAYS)) to Julian day number (DAYJ) for use by Subroutine 
SOLORB. 


See Table 4 for a summary of inputs and outputs for JULIAN. 


Table 4. Summary of JULIAN Input/Output Variables. 


INPUT VARIABLES 
Argument List 


IYRS - Number of the year in the 1900's (e.g., 1974 
becomes 74) at Greenwich 


IMONS - Number of the month (e.g., February becomes 2) 
at Greenwich 


IDAYS - Day of the month at Greenwich 
Common 


None 


OUTPUT VARIABLES 
Argument List 


YRFJ  - Julian day number (a half integer) at 0 hrs UT on 
January 1 of the year of interest 
VEQJ  - Julian date for vernal equinox 
DAYJ - Julian day number (a half integer) at 0 hrs UT on 
the day of interest 
Common 
None 
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3.3 SUBROUTINE SOLCYC 


Subroutine SOLCYC computes the 10. 7-cm solar flux (SBAR), 
an input to ATMOSU through ATMOUP Common, based on an assumed 
sinusoidal 11-year (or 4018-day) variation. The maximum value of 
250 for SBAR, associated with Model 9 of the CIRA-65 atmosphere has 
been assigned the date of 1 June 1958. The minimum value of 65 for 
SBAR is associated with Model 1 of the CIRA-65 atmosphere. 


See Table 5 for a summary of inputs and outputs for SOLCYC. 


Table 5. Summary of SOLCYC Input/Output Variables. 


INPUT VARIABLES 
Argument List 


DAYJ  - Julian day number (a half integer) at 0 hrs UT on 
the day of interest 


Common 


None 


OUTPUT VAKIABLES 
Argument List 


None 
ATMOUP Common 
SBAR- - Average 10.7-cm solar flux |1. 0OE-22 W/(m" Hz) | 
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3.4 SUBROUTINE SOLORB 


Subroutine SOLORB computes the north latitude (SOLLAT) and 
east longitude (SOLLON) of the apparent (actual motion) subsolar point, 
given the Julian day number at 0-hours UT on 1 January of the year of 
interest (YRFJ), the Julian date at which vernal equinox occurs (VEQJ), 
the Julian day number at 0-hours on the day of interest (DAYJ), and the 


universal time (UT). 


See Table 6 for a summary of inputs and outputs for SOLORB. 


Table 6. Summary of SOLORB Input/Output Variables. 


INPUT VARIABLES 
Argument List 


YRFJ  - Julian day number {a half integer) at 0 hrs UT on 
January 1 of the year of interest 

VEQJ - Julian date for vernal equinox 

DAYJ - Julian day number (a half integer) at 0 hrs UT on 


the day of interest 
TIME Common 


UT - Universal time corresponding to zone time ZT 
(decimal hours) 


OUTPUT VARIABLES 
Argument List 
SOLLAT- North latitude of subsolar point (radians) 
SOLLON - East longitude of subsolar point (radians) 
TIME Common 


GAT - Greenwich apparent time (decimal hours) 
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3.5 SUBROUTINE SOLZEN 


Subroutine SOLZEN computes COSCHI, the cosine of the solar 
zenith angle at a point P, given the geographic north latitude (PLAT) 
and east longitude (PLON) of the point P and the north latitude (SOLLAT) 
and east longitude (SOLLON) of the subsolar point. The day-or-night 
parameter IDORN is +1 for daytime, i.e., if COSCHI 2 0.0, andis -1 
for nighttime. The local apparent time (HL) is also computed from the 
Greenwich apparent time (GAT) and the east longitude of the point P 
(PLON). 


See Table 7 for a summary of inputs and outputs for SOLZEN. 


Table 7. Summary of SOLZEN Input/Output Variables. 


INPUT VARIABLES 
Argument List 
SOLLAT - North latitude of subsolar point (radians) 
SOLLON - East longitude of subsolar point (radians) 
TIME Common : 


PLAT - North latitude of point P (say, grid origin) 
(radians) 
PLON - East longitude of point P (radians) 


OUTPUT VARIABLES 
Argument List 


None 
ATMOUP Common 
IDORN - Parameter for day or night. If COSCHI is the 


cosine of the zenith angle of the sun at point P, 
IDORN is 1 for daytime, i.e., IF(COSCHI. GE. 0. 0), 
and is -1 for nighttime, i.e. , IF(COSCHI. LT. 0. 0) 


HL - Local apparent time (decimal hours, e.g., 2230 
hours becomes 22. 50 hours) 


4. MINOR NEUTRAL SPECIES 


4.1 SUBROUTINE SPCMIN 


The ROSCOE high-altitude chemistry module ‘Vol. 11] requires 
the minor neutral species O, CO,, N, and NO. Analytic-fit profiles 
for day and night at all altitudes are provided for O and CO, in Sub- 
routine ATMOSU. The profiles for N and NO are provided in Subroutine 


SPCMIN. 


The ROSCOE low-altitude chemistry module "Vol. 117 requires 
in addition to O, CO,; N, and NO, the minor neutral species H,0, 
0,('d,) Oz, and NO,, which are also provided by SPCMIN. 

The inputs and outputs for SPCMIN are summarized in Table 8. 
The nature of the functions used for fitting the adopted data-base values 


"Vol. 14b7 at noon or midnight in various altitude ranges is given in 


1 
Tables 9 through 16 for O, CO,, N, NO, H,O, O,( 4), Ox, and NO,- 


Table 8. Summary of SPCMIN Input/Output Variables. 


INPUT VARIABLES 
Argument List 


KK - Calculation flag 
=1, calculate initialization parameters 
= 2, calculate atmospheric properties 


ZH - Altitude of interest (km) 
ATMOUP Common 
IDORN - Index for day or night 
=+1, day 


=-1, night 


(cont'd) 


S1Z2N 


DATA 
ALTKM(47) 


NALTOD 


NALTND 


NDEGND 


NALTNN 


NDEGNN 


NALTNO 


NDEGNO 


NKMH20 


NDGH 20 


NALTO2 


NDGO2D 


NKMNO2 


ALTODN Common 


Table 8. (Continued). 


No density at 230-km altitude for use in 
N-density initialization. 


Altitudes at which minor species densities are 
specified as data 


Number of altitudes between 0 and 130 km used 
to establish the arithmetic function used for 
daytime O densities between 0- and 120-km 
altitude. 


Number of altitudes between 40 and 230 km 
used to fit the daytime N densities. 


Degree of the polynomial used to fit daytime N 
densities between 40- and 230-km altitude. 


Number of altitudes between 85 and 230 km 
used to fit the nighttime N densities. 


Degree of the polynomial used to fit nighttime 
N densities between 85- and 230-km altitude. 


Number of altitudes between 0 and 120 km used 
to fit the daytime NO densities. 


Degree of the polynomial used to fit daytime 
NO densities between 0- and 120-km altitude. 


Number of altitudes between 0 and 120 km used 
to fit H9O densities. 


Degree of the polynomial used to fit the Hg0 
densities between 0- and 120-km altitude. 


Number of altitudes between 0 and 50 km used 
to fit daytime Og(ta g) densities. 


Degree of the polynomial used to fit the daytime 
O(1Ag) densities between 0- and 50-km altitude. 


Number of altitudes between 0 and 160 km used 
to fit the daytime NO9 densities. 


(cont'd) 


NDGNO2 


ODA Y (27) 
ONITE(18) 
CO2(25) 
ANODAY (25) 
ANONIT(18) 
ANDAY(47) 
ANNITE (47) 
O2SDGD (47) 
O2SDGN(47) 
O3DA Y (27) 
O3NIT (27) 
H2ODN (25) 
SNO2D(33) 


SNO2N(33) 


* 
See Vol. 14b. 


Table 8. (Continued). 


Degree of the polynomial used to fit the day- 
time NOo9 densities between 0- and 160-km 
altitude. 


Noontime data-base values of [O| at altitudes 
0(5)130 km* 


Midnight data-base values of [O| at altitudes 
0(5)85 km * 


Data-base values of [CO9] at altitudes 0(5)120 
km * 


Noontime data-base values of [NO] at altitudes 
0(5)120 km * 


Midnight data-base values of [NO] at altitudes 
0(5)85 km* 


Noontime data-base values of [N]| at altitudes 
0(5)230 km * 


Midnight data-base values of {[N] at altitudes 
0(5)230 km * 


Noontime data-base values of (O2(1Ag)| at 
altitudes 0(5)230 km * 


Midnight data-base values of [O9(1A,)| at alti- 
tudes 0(5)230 km * 


Noontime data-base values of [03] at altitudes 
0(5)130 km* 


Midnight data-base values of |O3] at altitudes 
0(5)130 km * 


Data-base values of [HO] at altitudes 0(5)120 
km * 


Noontime data-base values of |NOg| at altitudes 
0(5)160 kni* 


Midnight data-base values of |NO9g| at altitudes 
0(5)160 km * 


(cont'd) 
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Table 8. (Continued). 


OUTPUT VARIABLES 
Argument List j 


None : 
ATMOUP Common | 

SNI( 7) -N density, 1/em® | 

SNI( 8) - NO density, Pom 

SNI(13) ‘ 0,(4,) density, 1/em® 

SNI(14) - 0, density, Lew” 

SNI(15) - NO, density, tiem 

SNI(16) - H,O density, i /em” 
ALTODN Common 

NALTOD - See input 

ALTKM(47) - See input 

ODA Y(27) - See input 

ONITE(18) - See input 

CO2(25) - See input (Note that the CO9g densities from 0- 


to 100-km altitude are reset in Subroutine 
ATMOSU by using a constant mixing- ratio of 
3.2% 10°4,) 


Table 9. Fit Functions for O Profiles.” 


Altitude Range, 
km Description 


Day 
0-120 Not O but f is fitted by a 12th-degree polynomial 
> 120 ATMOSU high-altitude model 
Night 
0-60 Constant at data-point value 
60-75 Exponential, with slope determined by data points 
at 60 and 75 km 
75-90 Sth-degree polynomial, to match data points at 
75(5)85 km and daytime fit-function at 90 km and 
derivatives of 60-to-75-km fit-function at 75 km 
and daytime fit-function at 90 km 
90-120 Daytime fit-function 
> 120 ATMOSU high-altitude model 


* Fits are made not in SPCMIN but in ATMOSU. 


Table 10. Fit Functions for CO, Profile. 


Altitude Range, 
km Description 


Day or Night 


0-100 Constant mixing ratio of 0.00032 in ATMOSU low- 
altitude model 
100-120 6th-degree polynomial, to match ATMOSU low- 


altitude-model value at 100 km and data points at 
105(5)120 km and derivatives of low-altitude-model 
function at 100 km and ATMOSU high-altitude- 
model function at 120 km 


> 120 ATMOSU high-altitude model 
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Table 11. Fit Functions for N Profiles. 


Altitude Range, 


km Description 
Day 
0-40 Constant at data-point value 
40-230 8th-degree polynomial, determined by least 
squares for data points at 40(5)230 km 
= 230 Proportional to No," 
=. t 
IN] = {ENV ] 9] 
Night 
0-85 Constant at data-point value 
85-230 6th-degree polynomial, determined by least 
squares for data points at 85(5)230 km 
= 230 Proportional to No," 
IN] = {INVIN,]} — LNg] 
230 


“This procedure makes [N] dependent on the time to the extent that 
[No] is dependent on the time. 
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Table 12. Fit Functions for NO Profiles. 


Altitude Range, 


| km Description 

| Day 

0-120 12th-degree polynomial, determined by least 
squares for data points at 0(5)120 km 

120-125 Parabolic transition function, determined by the 


120-km data point, the derivative at 120 km (esti- 
mated by uSing the daytime fit-function values at 
115 and 120 km) and a prescribed value at 125 km. 
(The 125-km value is determined by the require- 
ment that the slope of the function be continuous at 
125 km. See Vol. 14b.) 


> 125 Exponential, determined by the prescribed value | 
at 125 km and a solar-flux dependent value at | 
215 km. 
Night | 
0-50 Constant at data-point value | 
50-55 Exponential, determined by data points at 50 and 
55 km | 
55-85 8th-degree polynomial, to match data points at 


55(5)80 km, daytime fit-function value at 85 km, 
and derivatives of 50-to-55-km fit function at 
55 km and daytime fit-function at 85 km 


85-100 Daytime fit-function 


> 100 A prescribed altitude-dependent fraction of the 
daytime fit function 


Table 13. Fit Functions fcr H,O Profile. 


Altitude Range, 


km 


Description 


Day or Night 


0-120 


= 120 


Altitude Range, 


= 100 


Day 


12th-degree polynomial, determined by least 
squares for data points at 0(5)120 km 


Exponential, 


[HO] ss [50] a9 exp[-0. 166(h - 120)] , 


where [Hg0]};99 is determined from the fit func- 
tion from 0 to 120 km. 


Table 14. Fit Functions for 0,('A,) Profiles. 


Description 


10th-degree polynomial, to match data points at 
0(5)50 km 

Exponential, determined by data points at 50 and 
75 km 


oth-degree polynomial, determined by data points 
at 75(5)90 km and derivatives of 50-to-75 km fit- 
function at 75 km and 2 90-km fit-function at 90 km 


Exponential, determined by data points at 90 and 
105 km 


Night 


Constant at data-point value 


Exponential, determined by data points at 70 and 
80 km 


dth-degree polynomial, determined by data points 
at 80(5)95 km and values of daytime fit-function and 
its derivative at 100 km 


Daytime fit-function 
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Table 15. Fit Functions for O, Profiles. | 


Altitude Range, 
km 


0-40 


40-75 


75-90 


= 90 


Night 


0-55 
55-70 


70-75 


75-90 


= 90 


Day 


Description | 


9th-degree polynomial, to match data points at 
0(5)40 km and derivative of 40-to-75-km fit- 
function at 40 km 


Exponential, determined by data points at 40 and 
75 km 


5th-degree polynomial, to match data points at 
75(5)90 km and derivatives of 40-to-75-km fit- 
function at 75 km and 2 90-km fit-function at 
90 km 


Exponential, determined by data points at 90 and 
105 km 


Daytime fit function 


Sth-degree polynomial, to match daytime fit- 
function at 55 km, data points at 60(5)70 km, and 
derivatives of 0-to-55-km fit-function at 55 km and 
70-to-75-km fit-function at 70 km 


Exponential, determined by data points at 70 and 
75 km 


5th-degree polynomial, to match data points at 
75(5)90 km and derivatives of 40-to-75-km fit- 
function at 75 km and = 90-km fit-function at 90 km 


Exponential, determined by data points at 90 and 
105 km 
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Table 16. Fit Functions for NO, Profiles. 


Altituae Range, 
km Description 


0-160 12th-degree polynomial, determined by least 
squares for data points at 0(5)160 km 


> 160 Exponential, with slope determined by fit-function 
values at 140 and 160 km, and passing through fit- 
function value at 160 km 


Night 
oe [NO2 | ight % INO} say 8 [NOg lay : [NO] ight 
55-65 Exponential, with slope determined by fit function 
at 55 km, and passing through data point at 65 km 
65-82 Exponential, with slope determined by data point 
at 65 km and by daytime fit-function value at 
82-km altitude 
> 82 Daytime fit function 


4.2 AUXILIARY SUBROUTINES 


A brief description of the operation of Subroutines FITTER 
and SOLVE is given in Section 1. 


Ca | Subroutine FITTER 


A summary of inputs and outputs for Subroutine FITTER is 


given in Table 17. 


4.2.2 Subroutine SOLVE 


A summary of inputs and outputs for Subroutine SOLVE is 


given in Table 18. 


Table 17. Summary of FITTER Input/Output Variables. 


INPUT VARIABLES 


Argument List 


NPTS - Number of data points 

X(I) - Values of the independent variable, e.g., 
altitude (km) 

Y(I) - Values of the dependent variable, e.g., Species 
concentration (em73) 

NO - Degree of polynomial to be fitted. 

IKIND - Index for kind of equation to be fitted. 


NO 
=1 if equation is In(¥) = }> A, =< 
n=0 


NO 

; et ~ n 

= 2 if equation is Y= An D4 
n= 


ISIGN - Index for sign of exponents 
= 1 for negative exponents 
= 2 for positive exponents 

Common 


None 


OUTPUT VARIABLES 
Argument List 


Z(J) - The least-squares fit coefficients. Z(1) corre- 
sponds to A, Z(2) to Aj, ete. 
Common 
None 
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Table 18. Summary of SOLVE Input/Output Variables. 


INPUT VARIABLES 
Argument List 


A(I, J) - Element (I,J) of matrix of constant coefficients 
for NO simultaneous linear algebraic equations 
NO - The number of equations 
Common 
None 


OUTPUT VARIABLES 
Argument List 


X(K) - The least-squares fit coefficients. These are the 
same as the output Z(K) from FITTER and the 
same as DD(K) in ATMOSU. 


4.3 PLOTS OF MINOR NEUTRAL SPECIES PROFILES 


Comparisons of the fit-function values with the data-base 
values [ Vol. 14b7 of minor species densities are given in Figs. 3 through 
10. Broken lines (solid for noon, dashed for midnight) connect data- 
base values at 5-km intervals. Circles represent fit-function values at 


5-km intervals. 


Plots of day and night values for five of the minor species den- 
sities of interest in the D region appear in the very recently published 
handbook by Knapp and Schwartz |KS-75, Fig. 8-1], reproduced here 
as Fig. 11. To aid in comparing our results with the handbook results, 


we have used their scales to replot our fit functions for the same spe- 


cies, shown in Fig. 12. 
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Fig. 11. D-Region Minor Neutral-Species Profiles 
from KS-75 (Fig. 8-1). 
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Fig. 12. Selected Minor Neutral-Species Profiles from 
Subroutines SPCMIN and ATMOSU. These 
profiles may be compared with those in Fig. 11. 
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5. AMBIENT IONOSPHERE (SUBROUTINE IONOSU) 


See Table 19 for a summary of inputs and outputs for I[ONOS| 


Dot E- AND F-REGION IONOSPHERIC PROPERTIES 
The E- and F-region chemistry module requires the following 
quantities: 


a. q, the effective total ion production rate that reproduces 
the ambient ionosphere when used with the chemistry 


model (em? eee”) 


ive ous , -3 
b. O', the positive atomic-ion density (cm —) 
es -3 
c. M _, the positive molecular-ion density (cm —) 
d. TL» the electron and (No vibration) temperature ( K) 


| RANC IV uses unsatisfactory steady-state formulas |GE-70, 
Eqs. (2-276), (2-274), and (2-275)| for q, O , and M.. The derivation 
|GE-70| of the RANC IV equations is incomplete because there are a 
number of tacit, (unnecessary) simplifying assumptions. The exact 


solutions may be derived as follows. 


Equations (2-270) and (2-271) of GE-70, for steady-state con- 


ditions, become 


[O"] = 0 = q, - BO"|- a, [0°] jlo" | + (M"] (1) 
[M*] = 0 = dy + B[O"]- @,[M~} {[o"] +[M"]} , (2) 
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Table 19. Summary of IONOSU Input/Output Variables. | 


INPUT VARIABLES 
Argument List 


Jd - Calculation flag 
yd JJ=1: calculate initialization parameters 
| JJ=2: calculate ionospheric properties 
ZH - Altitude of interest (km) 
ATMOUP Common 
IDORN - Parameter for day or night. If COSCHI is the 


cosine of the zenith angle of the sun at point P, 
IDORN is 1 for daytime, i.e. , IF(COSCHI. GE. 
0.0), and is -1 for nighttime, i.e., 
IF(COSCHI. LT. 0. 0) 


SNI(1) - Ny concentration (1/em*) 
SNI(2) - O, concentration (1/em*) 
SNI(3) - O concentration t/a’) 
Jidh - Heavy-particle temperature ( K) 


ALTODN Common 


ALTKM(47) - The array of altitudes at which minor species 
are specified as data in SPCMIN 


RATE Function Routine 


Reaction rate coefficients for chemical reac- 
tions (This routine prepared by Knapp and 
Jordano £ Vol. 111.) 


DATA 
HE BOTD - Altitude below which the daytime electron den- 
sity decreases exponentially and above which 
the logarithm of the daytime electron density 
increases parabolically (km) 
EBOTD - Daytime electron density at altitude HE BOTD 


(1/em3) 


(cont 'd) 


68 


HF2MXD 


EF2MXD 


EDDSCH 


F2DSCH 


HE BOTN 


EBOTN 


HF2MXN 


EF2MXN 


EDNSCH 


F2NSCH 


TXT120 


TXT200 


TXT800 


DQDAY(18) 


Table 19. (Continued). 


Altitude at which the maximum daytime elec- 
tron density occurs (km) 


Daytime electron density at altitude HF2MXD 
(1/em3) 


Seale height with which the daytime electron 
density decreases below altitude HE BOTD (km) 


Scale height with which the daytime electron 
density decreases above altitude HF2MXD 


Altitude below which the nighttime electron 
density decreases exponentially and above 
which the logarithm of the nighttime electron 
density increases sinusoidally (km) 


Nighttime electron density at altitude HE BOTN 
(1/em3) 


Altitude at which the maximum nighttime elec- 
tron density occurs (km) 


Nighttime electron density at altitude HF2MXN 
(1/em3) 


Seale height with which the nighttime electron 
density decreases below altitude HEBOTN (km) 


Scale height with which the nighttime electron 
density decreases above altitude HF2MXN. 


The difference between the electron tempera- 
ture and the gas temperature at 120-km alti- 
tude in the ambient daytime ionosphere ( K) 


The difference between the electron tempera- 
ture and the gas temperature at 200-km alti- 
tude in the ambient daytime ionosphere ( K) 


The difference between the electron tempera- 
ture and the gas temperature at 800-km alti- 
tude in the ambient daytime ionosphere ( K) 


The effective total ion production rate at alti- 
tudes 0(5)85 km that reproduces the ambient 

daytime D-region ionosphere when used with 
the chemistry model (ion pairs {cm sec }) 


(cont'd) 
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DQNIT (18) 


OUTPUT VARIABLES 


ATMOUP Common 


SNI(9) 
SNI(10) 


SNI(11) 


SNI(12) 


IONOUP Common 


EFE 

E FOP 
EFMOLP 
TX 
QDEF 


Table 19. (Continued). 


The effective total ion production rate at alti- 
tudes 0(5)85 km that reproduces the ambient 
nighttime D-region ionosphere when used with 
the chemistry model (ion pairs/|em” sec }) 


Electron concentration for ZH = 90 km (1, ea) 


Atomic oxygen ion concentration for ZH = 90 km 
(1/em”) 


Molecular ion concentration for ZH = 90 km 
(1/em3) 


Electron (and Ng vibration) temperature ( K) 


See SNI(9) above. 
See SNI(10) above. 
See SNI(11 

See SNI(12) above. 


) above. 


The effective total ion production rate that re- 
produces the ambient ionosphere when used 
with the chemistry model 
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where 
[O*] = positive atomic-ion density (em7°) 
ven = positive molecular-ion density (em™*) 
-1 


q, = atomic-ion production rate (em7? sec 


do = molecular-ion production rate (em ~ sec ~) 


B = KIN, | + KO, | 


ky = jon-molecule interchange rate coefficient lene eee) 
Ky = ion-molecule charge-exchange coefficient fom /aec} 
oS effective two-body collisional-radiative recombination rate 
coefficient for atomic ions (em’/sec) |KJ-74b| 
2 tae! fait eee 
= ky (TY) + ky 9(T,) [e]+1.5x10 © [e] 1p 


Ky, (7) = yadiative recombination rate coefficient for the reaction 
O° +e>O 4 hy 


k,,(T_) = collisional-radiative recombination rate coefficient for the 
reaction O'+e+e7O+e 


an dissociative recombination rate coefficient for the reaction 
M' + e> products (em?/sec). 


By assuming charge conservation, 


[e] = [O°] +[M'] , (3) 


Let 


% 1 > (6) 
where 
Ea em 
and 
a= Sy + Hy + (8) 


By adding Eqs. (4) and (5) and using Eq. (8), we have 
q = {aIM*]+a,[07}t[e] . (9) 


We have five equations |(3), (4), (5), (6), and (8)| and six vari- 
ables; we can get analytic solutions for the five variables if we regard 


the electron density as given. 


Equations (4) and (6) give 


or 


q = {8 +a felt to*] (11) 
Use Eqs. (9) and (11) to eliminate q: 
JaglM*]+o10"]} le] = {e+e fel}r or) . (12) 
Eliminate [M*] from Eq. (12) by using Eq. (3) to get 
; af el? 
lO} = (13) 


B + le|[agf + @(1 = f)] 
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Eliminate rom from Eq. (11) by using Eq. (13) to get 


i tet yar lt 
a, le|~ 51+ [e|B 
PES: ali ae en ea (14) 


1 + |e] gi {af + a (1 - ft 


Thus we have expressed q as a function of |e]. This equation differs 
from the corresponding equation in RANC IV |GE-70, Eq. (2-276)] by 


containing terms involving a. 


The reader who likes Eq. (13) for [oO] can use it, of course, 


; == 
but we have used another expression for |O' |, 


iS 
+ 
| 
> 
2 
_ 
_ 
or 
~_— 


obtained by solving Eq. (11) for (o" |, since q is known from Eq. (14). 


From Eqs. (5), (8), (6), and (3) we get 


It may be shown that the sum of [o" | and [M’], as given by Eqs. (15) 
and (16), satisfies the requirement of charge conservation expressed by 
Eq. (3); this is not true for the corresponding RANC IV equations 
|GE-70, Eqs. (2-274) and (2-275) |. 


In Subroutine IONOSU we use Eqs. (14), (15), and (16) to com- 


pute q, [O |, and [M’] after prescribing analytic fits to nominal profiles 


ee Te ee ee ee ee 


of E- and F-region electron density |Ri-73, Fig. 1] and electron tem- 


perature |Ev-73]. 


The prescribed electron-density profiles in the E- and F-region 
for noon and midnight conditions are shown in Fig. 13. The fit functions 


used to obtain these profiles are described in Table 20. 


The prescribed electron temperature profile and the heavy- 
particle temperature profile in the E- and F-region for noon and mid- 
night conditions are shown in Fig. 14. The fit function used to obtain 
the electron temperature profile is described in Table 21. 

For approximately mean solar-flux conditions, SBAR~ S 
158 x 1 Ww ne Hat, profiles of q are shown for noon and midnight 
conditions in Fig. 15 and the corresponding values of O* and M’ are 


shown in Fig. 13. 


0.2 D-REGION IONOSPHERIC PROPERTIES 


The D-region chemistry module requires the following 
quantity: 

q, the effective total ion production rate that adequately 

reproduces the ambient ionosphere when used with the 

chemistry model. 

The modeling of q in the D region (and lower) is offered with 
reservations; it may need to be improved if experience shows that this 


topic is more important than it is presently assumed to be for radar. 


For the D region, q is determined by specifying data points at 
30- and 60-km altitude and by requiring the fit function to be continuous 
with the value of q derived from the E- and F-region model at 90-km 
altitude. The fit function is extrapolated below 30-km altitude for 


modeling convenience and not on a physical basis. 
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Table 20. Fit Functions for E- and F- Region 


Altitude Range, 


km 


90-100 


100-300 


> 300 


90-100 


100-360 


> 360 


a P ; ee 
Based on Fig. 1 in Ri-73. 


Day 


Night 


Electron Density Profiles. 


Description 


Exponential, determined by data-point value 
(EBOTD) at 100-km altitude (HEBOTD) and scale 
height EDDSCH 


Parabola, determined by data-point values EBOTD 
and EF2MXD at altitudes HEBOTD and HF2MXD 
and vertical slope at altitude HF2MXD 


Exponential, determined by data-point value 
(EF2MXD) at 300-km altitude (HF2MXD) and 
scale height F2DSCH 


Exponential, determined by data-point value 
(EBOTN) at 100-km altitude (HEBOTN) and scale 
height EDNSCH 


Sinusoid, determined by data-point values EBOTN 
and EF2MXN at altitudes HEBOTN and HF2MXN 
and vertical slope at the same altitudes 


Exponential, determined by data-point value 
(E F2MXN) at 360-km altitude (HF2MXN) and scale 
heigt F2NSCH 
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Table 21. Fit Function for Electron Temperature Profile. 


Altitude Range, 


km Description 
Day 
=< 120 Same as heavy-particle temperature 
= 120 The difference between the electron temperature 


(Tx) and the gas temperature (T) is prescribed to 
be zero at 120-km altitude and 500 K at 200-km 
altitude. The parabola 


z 
< 


Te - T = 500 [(ZH - 120)/80| 
is then used. 
Night 


= Q Same as heavy-particle temperature 


The data adopted are based on the calculations of Webber 
|We-62| for the ion-production rate due to galactic cosmic rays. Webber 
|We-62, Fig. 2] presents results in the altitude range from 30 to 90 km 
for two geomagnetic latitudes (50 and 70°) and for sunspot-minimum 
and sunspot-maximum conditions. For the geomagnetic latitude of 50 , 
Webber |We-62] finds Vax 0°04 and q.j, = 0.08 at 60-km altitude 


mir 
and q =2.1l andq_. =4.5 at 30-km altitude. We adopt solar- 
max min 
cycle mean values of 0.06 and 3.3 at 60- and 30-km altitude, respec- 
tively. The interested reader may also wish to consult Ra-72 (Fig. 2.3) 


and Po-73a (Figs. 2 and 3). 


The profiles of q in the D and adjacent regions for noon and 
midnight conditions are shown in Fig. 16. The fit functions used to ob- 


tain these profiles are described in Table 22. 
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Table 22. Fit Functions for Effective Ion Production Rate 
in D and Lower Regions. 


Altitude Range, 
km Description 


0-60 Exponential, determined by data-point values at 
30- and 60-km altitude 


60-90 Exponential, determined by data-point values at 
60-km altitude and daytime value of q from E- 
and F-region model at 90-km altitude 


Night 
0-60 Same as daytime 


60-90 Exponential, determined by data-point value at / 
60-km altitude and nighttime value of q from E- 
and F-region model at 90-km altitude 


81 


—— = *- -- . 


6. DRIVER, LISTING OF COMPUTER PROGRAM, 
AND SAMPLE PROBLEM RESULTS 


A short driver routine is provided to exercise the ATMOSU, 
SPCMIN, IONOSU, and associated routines. The required input con- 
sists of the year, month, day, zone time, geographic colatitude and 
longitude of the point of interest, a set of test altitudes, and the number 
of such altitudes. Input quantities are more specifically described be- 
low in Table 23. The driver, after reading and writing the input, first 
initializes the ATMOSU routine by the call ATMOSU(1,120.). The 
driver next loops over the test-altitude array, exercises the ATMOSU, 
SPCMIN, and IONOSU routines for each altitude, and prints out the re- 


sultant atmospheric and ionospheric property values. 


A listing of the driver, ATMOSU, SPCMIN, IONOSU, their 
associated subroutines, and outputs from two sample problems are 


provided. 


The quantities in the output block at each altitude are labeled 
in the headings. The last four entries (E, oO", M , and N_) in row-two 
of the output block at each altitude (= 90 km) are computed by Subroutine 
CHEMQ and are included for comparison with the quantities E, O', and 
M’ in row-one computed by IONOSU. Subroutine CHEMQ, prepared by 
Knapp and Jordano see Vol. 11) for use with the NRL Simple- 
Chemistry module, computes steady-state ionization for the E- and 


F-region. 


Table 23. Input Quantities to DRIVER 


a. NALTS - Number of test altitude values 
(FORMAT-I5) 
b. ALTA(I) - Test altitude values (km), eight values per card; 


NALTS values in total 
(FORMAT-8F10. 2) 


' c. F¥RS - Number of the year in the 1900's at east longi- 
' 4 tude GLO (e. g., 1974 becomes 74) 
: IMONS - Number of the month at east longitude GLO 
(e.g., February becomes 2) 
IDAYS - Day of the month at east longitude GLO. 
ZT - Zone time for the 15-degree longitude interval 


containing east longitude GLO. 


GCO - Geographic colatitude of grid origin or whatever 
reference point is desired (degrees) 


GLO - Geographic east longitude of grid origin or 
whatever reference point is desired (degrees) 


(FORMAT 315. 3E10. 4) 
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53 


71 


173 


203 


aq? 
ei! 

213 
ais 
A? 
Pe} 

22? 
Pe? 
22h 
Per 
ae? 
230 
24? 
ass 
235 
P%A 
ey? 
aun 
Put 
Pur 


Oro 


WRITE (56,2002) 
Pune FORMAT (5x,4( 
KPITE (6,2904) 
PUN0G FORMAT (//,AH 
& Ra 7Y sfj 
e UA DPFEG) 
CONVERT G 
G00 = GCOerR AD 
GIO = GLOI®RAD 
INENTIFY 
PLAY & PINPeG 
PIUN & 6LN 


* INTTYVAl T7E 


WRITE (6,8029) 
R020 FORMAT(//20H 


CALI ATMOSUC] 


wEITE(H,2°9046) 
20%% Fi~MAT €/7,84 
* AM 77 SF} 
# Ue RAD) 
AHRITE(6,2905) 


CT,ALISCT),»L32,NAL TS) 
TR, 2X,F19,2)) 
TY2®S,TMINS, IDAYS, 27,601, 6L9 


TYRS S150) Orn IMONS #275,19h [NAYS sf5/ 
2,4,14H HRS GCM st12,4,14H DFG GIO sF1A,.4, 
Ch) aND BLO FRUM SEGREES TU KADIANS, 

DEG 
NEG 
THE GRIT) ORIGIN AS THE POINT P, 
Ci) 


THE ATMNSU ROUTINE 


TNTTIALIZATION CALLS S/S) 
01?0.) 

TYRSeIMONSe INAYSsZT¢509,6L0 

TYRS 315,)0H IMUNS 315,104 IDAYS s15/ 
2,4,)4HK HRS GCN zE12,4,14H RAD GLO) SE Lets 


TOORN,HTsGATSPLAT,PLON 


DRIVER, 58 
DRIVER,59 
DRIVER, 60 
DRIVER A} 
ORIVER be 
URIVER AS 
DRIVER, 44 
OR IVER, 6S 
ORIVER 66 
DRIVER, A7 
DRIVER AB 
UR IVER, 49 
URIVER,70 
OFIVER,71 
DRIVER,72 
DRIVER,73 
URI VER, 74 
DRIVER, IS 
VRIVER, 76 
DRIVER, 77 
ORIVER, 7A 
DRIVER, 79 
DRIVER, ARO 
UPIVER,AY 
ORI VER Re 


2ONS FORMAT (//,8H TRURN 
* PLAT SE17,4,10» 


efS,10n UT 8F1{2,4,104 
PLON sE12,.4) 


GAT SbE1eede1 On 


UORIVER 


BS 


WRITE CH,29GBYHL , SAAR 


PONS FORMATO //,5H HL s,F10,3e5Xe7H SAAR &,F10,3) 


t *& |(NP QVER TEST ALTITUDFS 

wRITE(6,8002) 

BONE FORMAT (1H0,17P9H ALY N2 O02 nN 
@ AR He C2 E n+ rae: 

QNEF 71O0%09(S¥, GHIsECs S¥),ek, 10H (CC SEC)) 

WRITE (6,8093) 

RONS FORMAT F1HO,9¥%, 1204 N NI) NUD 
* 3 Han F Ne Me 
* Ne A19xX,1905x,444/C0, $%)) 
wGITE (AH, ANG) 

ROAM FORMAT (440,9%,72H PRESSIIRE FERSEQ DENSITY 
a TEMP F TEMP /10¥,724 NYNES/CM¥82 
eC ra PEG * DEG « ) 


ON SO TEteNALTS 

74 s AaLTser) 

CALI ATM ISUCA,7H) 

CALL SPCMIN( 2,74) 

Call TONUSUCA+ 7) 

ENEQN 2 0,0 

Pls 2,0 

ENPN 2 0,f 

{Fe 2H bt oh 3} GO TH 8s 
CNe2 & SNIC}) 

the 3 § 

Ci 2 SNY(3) 
C1) s SNIfR) 
CNas s SNT(7) 
CNeA = 1,8 
CNP = 0,0 
Cre & O06 
CENF 2@ 0,6 
¥¥ 2 TX 

re Ss FR 


DRIVER,AS 
DRIVER, RO 
DRIVER AT 
URIVER, RB 
VRIVER, AO 
DRIVER, 90 
DRIVER, 9] 
ORIVER, 92 
UDRIVER, 93 
OPI VER, 94 


e 

. 
ORIVER, AU 

. 


N2CSNCIVRIVER, 95 


DRIVER, 96 
DRIVER,97 
ORIVER, 98 


NEN SC HTDRIVER, 99 
GRAMS/CURIVER, 100 


VRIVER, 101 
DRIVER, 102 
ORTVER,103 
VETVER, 104 
UST VER, 105 
DRIVER, 104 
URIVER, 107 
URI VER, 108 
ORIVER, 109 
ORIVER, 119 
ORIVER, IAS 
URIVER, 112 
ORI VER, 1413 
ORIVER, 144 
DRIVER,NAS 
DRIVER,ILS 
DRIVER, 117 
DRIVER, ILA 
DRIVER,119 
URIVER, 120 
ORI VER, 12! 


Pu? TG. Tt ORIVER, 122 


2u4 CALL CHEMQ(QDFFE,ENPO,PQ,FNFQ) DRIVER, 123 
Pu? 45 EMPAQ & FNFEQeUPQWENPY URIVER,.1e4 
25? WPITE(C6,8005) PHACSNT( J) eJelebde CSNIC I), JE9, 11) ,90VFR,SNIC7), DRIVER, 12S 
* SNI(A),SNICIS) oSNICL3),SN1014),8NI 016) ,ENEGQ, MPO, ORIVER,126 

* FMPQ,ENPQ, PR, FEHSEN, PHD, HHH, TT, SNI( 12) ORIVER, 127 

ROOS FORMAT CIXFU.2,1NEL2,5/(10X%,10E1203)) DRIVER, 128 

B35 SO CONTINUE ORIVER,129 
C ORIVER,130 

San wQITE (6.9050) ORIVER, 131 
9050 FORMATC//,204 END UF TEST PRORLFEM) DRIVER, 1 $2 

{ buy GN TM yoyo DRIVER,133 
tug END URIVER 134 


ATMOSU 


SUAROLTINE ATMOSUCJJ, 74) ORIVER,155 

C ORIVER,V36 
C ATMOS!) POMPIITES THE PRIPERTTES MF THE UNUTSTURHKFL) ATMOSPHERE, DRIVE®8,137 
€ GIVEN THE OL TITIOF 7K. AFTER ASSUCTATED SIIBRULITINES ChiMPUTE URIVER 158 
Cc THE LOCAL APPARENT TIME HL, SILA FLUX SB4k, AND MAVeNKeNTIGHIT DRI VER, 1 $9 
is PARAMFTFR [TAQRN, URI VER, 140 
c ATMOSI! TS REVISTON 99 (96/02/75) OF ATMS DEVELOPED BY PB, wn, URIVER, 144 
C LOWFN (SEF, AN AMAITPNT ATMOSPHERE MIDEL FOR ROSCHEs FP, 187, ORT VER, 142 
C VOL, 5S OF PROC, UNA 497% ATMOSPHERIC EFFECTS SYMPNSTUM, DNA ORIVER, 148 
C 31359Pe5, S TUNE 1973,) URIVER, 144 
(3 REVISTIIN 92 (NOs07/74) PROVIDES URI VER, US 
( 1, IN HIGHeALTITUDE MUNEL, FOR USE UF GAaF(129,) INSTEAD OF DRIVER, 146 
f GaFkCO,) © G2 IN COMPUTING GAM AND 27, URIVER, 147 
c >, NFENSITY SCALE HFIGHT FUR RLITH LOW@ AND HIGH@AL TY TUDE VRIVER,IGA 
4 C MODELS, wTTH AN AN HOC PAaRARUI IC TRANSITION FROM 110° TODRIVER, 149 
C 120eKM ALTITUNE TN PROVIDE & CONTINUUNS DENSITY SCALE ORI VER,150 
( HEIGHT ACRIISS THE ROUNDARY RETWEEN THE Tw MUNELS, ORTveR, 151 
t 3B, ALTFRED FURMULA FOR 0 DENSITY UN FIRST CALL AND aT LiW URIVEP,152 
C ALTITUDE S11 aS TU USE SFeFUNCTINN DTRECTLY, DRIVER, 153 
C OU, COMMENT CARDS, DRIvFR, 154 
c REVISTUN 083 (10/25/74) PROVIDES DRIVER,155 
Cc S. PRIVTSION ENR DAY UR NIGHT VALUFS (UF ATOMIC OXYGEN ORIVER, ISA 
C (NBTATNED FROM THE MINOR SPECTES SURRIIUTINE SPCMIN) URI VER 157 
C FOR ALTITIIDES BELOW 120 KM, ORIVER,158 
3 6, AUTOMATED PRICENURE FOR EVALUATING CONSTANTS TN DENSITY ORIVER, 159 
C SCALFeHFIGHT FORMIILA USFU IN THE 11ne Til {enekm™ DRIVER, 1609 
c TRANSTTTUN PEGION, DRIVER, Jo} 
C 7, PRICENURE FOR LETTING SMLAK FLUX SHAR, AN INPUT Tf ORIVER, 162 
C ATMOS, BE NETERMTNED BY THE AUXILIARY RULITINE SUICYC, ORIVER, 163 
C A, PRICENURE FOR LFTTING THE LOCAL (APPARENT) TIME wl, URI VER, 164 
c aN TNPUT Th aATMOSII, BE AETERMINED BY THE AUXILIARY ORIVE®,.165 
C SUARMUTTINE SOLORK, ORIVER, 166 
€ 9, PRICENURE FIR LETTING THE DAY JR NIGHT PARAMETER TUDRN URIVER, 167 
C KE NR TERMINED BY THE AUXILIARY SURROUTINE SOLTES, DRIVER, 168 
C RFVISTON OW (12/0R/7U) PROVILFS DRIVER, 169 
C 19, CAHRUN DTD TDF AS THE SIXTH SPEFIFS IN ATMUSU, KITH DRIVER, 170 
C PROFILE SPECIFIED BY R, F, MYFRS [IN J2/07/7U, DRIVER, I71 
( 11, EVALUATTUN OF DEPARTURE FRUM HYDROSTATIC FQUILIRRTUM, ORIVER, 172 
t 12, & FLAG, ZHELAG, TM INSURE THatT SUBROUTINES FUNUSU AND ORIVER,173 
t SPCMIN ARE CALLED aT THE SAME ALTITUDE AT WHICH aTMOSt! ORIVER,174 
G a4§ {4ST CalLLFD, URIVER,175 
C 13—e DAY AND NIGHT PROFILES AF ATOMIC OXYGEN SPECIFIFD BY DRIVER 17H 
c Ke FL MYERS UN 11/09/74 AND 11/25/74, RESPECTIVELY, DRIVER, 177 
G 14,, CORRECTED PSLICEDURE FOR EVALUATING CONSTANTS TN DENSITY DRIVER, I7A 
C SCALFeHFIGHT FURMILA IISFD IN THE Iie TN {eneKm OR] VER,179 
C TRANSTTION REGION, DRIVER, 180 
¢ 15, CORRECTED CONSTANT IN LOweALTITIDE FURMULA FUR AENSTTY ORI VER, 14! 
t SCALE HEIGHT, ORLVER 182 
C RFVISTON 95 (Ne/yuus7ys) PROVIDES DRIVER, 1BY 
t 16, INTERFACE wYT SPCMIN WHICH Ntw COMPUTES NENSTTIES OF DRIVER, 184 
C H2tle Ne NO, NM Qe TOCSINGLET DELTA G),y AND OF, DRIVER, 18S 
C 17, INTERFACE WITH TUNOSU WHICH NOW COMPUTES THE FREECTIVE DPIVER,1BA 
C TON PRLUNUCTTON Rate at all aL TITUDES, DRIVER,187 
C REVISTON Nm (O4/0R/75) PROVIDES DRIVER,1HA 
t tRe REVISED NIGHT PRUFILE UF ATOMIC UXYGEN SPECTFITEN ey URIVER, 169 
C BeFese MYEFS UN 02/27/75 (MINOR CHANGE HELOw 69 KM), ORT VER,199 


NQAOKNKCOAOOHANMVAAMAAAAAGAANMNOCANARAANONAANMAE CH 


NAHAONOCCANMKOAKOAAANIAHAAAAAAAM 


en, 


REVISED) DAY AND NIGHT PRUFIL ES MF NY 
SPECTFIFD KY RLF, MYERS UN Nu /05/75, 
REVISED DAY AND NIGHT PRUFILES NF at 
SPECIFIED HY AF, MYERS UN fus117/7S, 


REVISTUN OF (04/24/75) PRIIVIDFS 


ele 


REVISED PRUFEDURE FOR SPECIFYING AND 
VFRNalL FUUINOX (PER RY We. LOWFN (927 


REVTSTUN MAR (95/23/75) PROVIDES 


2? 


REVISFH PROFILE UF WATER VaPOR SPECT 
ON OS/10/75, 


REVISTON 99 (06/02/75) PROVIDES 


23, 


INPIIT 


CORRECTED FORMULA IN HIGH@PALTITUDE M 
OF NEPARTURE FROM HYDRPUSTATYC ERUILY 
PARAMETERS 


ARGUMENT LIST 


JJ = CALCULATION FLAG 
Ss ie CALCULATE INTFITALT ZATION 


TRIC OXTOF 
UMIC NITROGEN 
USING NATE UF 
26/750), 


FIED BY B, F, 


UDEL FOR FVALUATION 


BRIUM, 


PARAMETERS 


= 2e CALCULATE ATMUSPHERIC PROPERTIES 


TR © ALTITIIDOF OF INTEREST (K™) 


ATMONUP COMM(IN 


He, SRA, ITNURN 


TIME COMMON 


IYRS, IFMUNS, TOAYS, ZT, PLAT, PL 


AL TORN COMMUN 


NALTUO, ALTKM(47), NDAY( 27), ONT 


(MITPUT PARAMETERS 
ATMOIIE COMMUN 


PP, RHO, TT, SNI (6), HRHO, FER 


AL TORN COMMUN 


COMMUNSATMIUPS 
COMMIINSTIMES 
COMMUNZALTORNA 


$1ZAN 


HL» SHARKS TDIIRN,PP,RMI,TT,SNI( 16 
TYRS,IMONS, ITDAYS,2T,PLAT,PLONUT 


VARTARLFS JN ATMQUP 


HE oe LOCAL TIME, HRS : 
SRAw = AVER, 10,7efm SMLAR FLUX, Jebe 
INQDRN = INDEX FOR DAY UR NIGHT, FOR U 
DAYTIME PRUF{LE JFCTOORN,GE.0) 
PRUFILE IF(TDORN,LT,9) 
PR = PERSSIIRE, DYNES/CMRRO 
RRO = DENSITY, G/T MaeZy 
TT = TEMPERATURE, DEGREES KELVIN 
SNI(1) #& NA, {/CMeey CFROM ATMNSU) 
SNY(A) = 02, {/Emeey (FROM ATMS) 
SNT(S) & OF {/CM*eZ CFROAM ATMOSU) 
SNT(4) = aR, t/CMPet (FROM ATMNSU) 
SVI(5) = WE, '/CMHRR (CF RIIM ATMINSI!) 
SNT(6) & CM2s 17CM**3 (FROM ATMNSU) 
SNICV) & Noe t/Omeay (FROM SPCMIN) 
SNT (A) = ND, 1/Cmeet (FROM SPCMIN) 
SNI(9) & Es {/CmMeey (FROM JONNISU) 
SNIC}I{9N) = Ue, t/Cmewt CFROM LONNSU) 
SNI(11) s Me, {/O meaty CF ROM TANS!) 


SNI (12) 


TX, AEG XK CEROM TONOSU) 


88 


ON 


TRE 


MYERS 


TEC18), CM2025) 


SEW 


)pHRHN, PERSE 
GAT 


22 W/( hae? KZ) 
BELOW 120 KM, 
AND NIGHTTIME 


NALTOM,ALTKM (U7) ,NDAY (27) ,UNITE CIA), S1Z2N,CN2025) 


USE 


DRIVER 194 
Ds 19? 
ORIVER,193 
DRIVER, 194 
ORIVER, $95 
ORIVER, 196 
DRIVER,197 
ORIVER, 198 
DRIVER,199 
ORIVER, 200 
ORIVER,201 
DRIVER, 202 
DRIVER,203 
ORTVER, 204 
DRIVER, 2US 
DRIVER, 206 
VDRIVER, 207 
URI VER, 20R 
DRIVER, 209 
DRIVER, 210 
DRIVER, 211 
DRIVER, 242 
ORIVER, 2143 
URIVER,2i44 
DRIVER, 215 
DRIVER, 216 
ORIVER 217 
DRIVER, 21 
DRIVER, 249 
DRIVER, 220 
DRIVER, 225 
DRIVER, 222 
DRIVER, 223 
URIVER,22u 
DRIVER, 225 
DRIVER, 226 
DRIVER, 227 
DRIVER, 2eR 
DRIVER, 229 
DRIVER, 239 
DRIVER, 231 
URLVER,232 
URIVER, 233 
URIVER 234 
DRIVER, 235 
URIVER, 23h 
DRIVER, 237 
DRIVER, 23R 
ORIVER, 239 
URIVER, 240 
DRIVER, Au} 
ORIVER, P42 
VRIVER. 2445 
DRIVER, ?u4 
ORIVER,245 
DRIVER, 2u6 


| 
| 
| 
| 


C SNIC13) s UACLIDR), I/0Meed (FROM SPCMIN) VR] VER, Au7 
c SNICY4) = (13, {/CMeey CR ROM SPCMIN) DRIVER, AUR 
Cc SNIC15) = NNM2, 1/CM*e d (FROM SPCMIN) VRIVER, 249 
Cc SNIC16) = Het, {yrmeaS (FROM SPCMIN) DRIVER, 250 
Cc HRH = DENSITY SCALE HEIGHT, KM DRIVER, 251 
Cc FEHSFy s FRACTYONAL FRROR IN HYDRUSTATIC EQUILIRRIIM, DRIVER, A252 
C ORIVER, 253 
C URIVER, 254 
COMMONS7THCHE XS ZHFLAG DRIVER, 255 

C URIVER, 256 
DIMENSION 4(6),805),0(4)65(5),8H05),A4A012).09(13) ORIVER, 257 
DIMENSION SNT7(63, SMI Cod, ALP fod DRIVER, 25" 
NIMENSION ACCO,21) Fe XCOIsXCC7),Z2ZI4A2ON(S) »-ZIIN(5S),ONZI (9) DRIVER, 259 
DIMENSION ZIMYC(S)6Z2ICMS(S),CIAzZT(5) ORI VER, 260 
NIMENSION FoAav(e?) DRIVER, 261 

€ DRIVER, 262 
C VEFINTTIONS UF DATA QUANTITIES DRIVER, 263 
C STGMS = SFAPLFVEL MEAN MULECULAR WEIGHT, G/ML E ORIVER, 264 
L P2 = SFAMLFVEL PRESSIIRF, DYNFS/CMa8e URIVER, 26S 
ig RIGA = AVUGANRRO NUMRFR, PARTYICLES/SOLE DRIVER, 266 
C RR os UNIVERSAL GAS CONSTANT, ERG/(MULE DEGew) ORIVER, 267 
C (SET TN SUBROUTINE, RReSKe5IGa) DRIVER, 26Rf 
C DRIVER, 2609 
DATA BIGMS,P7,R1GA / 29,96,1,01325FHNe,be022169EF23 / ORIVER,270 

Cc DRIVER, 274 
C Sk z ROLTZMANN CONSTANT? ERG/(DEG eK) DRIVER, 272 
¢ NOFG «© OFGREF MF POLYNOMTAL TH BE FITTED FOR THE DRIVER, 273 
C OAYTIME () PROFILE DRIVER,274 
C DRIVER, P75 
C CAUTION ewww NDEG MUST NOT EXCEED 12 wITHOUT MAKING DRIVER, 276 
C APPROPRIATE CHANGES IN PROGRAM, ORIVER, 277 
DATA PT,SkK £4 3,.141592653590,1,38Nb22E ein 4 , NDEG / 12 ORI] VER, 278 

C ORIVER, 279 
c GZ = SFA=LFVEL GRAVITATIONAL ACCELERATION, CM/SEC 82 DRIVER, 289 
C Qt = RADTUS OF SPHERTCAL EARTH, KM ORIVER, 281 
Cc DRIVER, 282 
j DATA G7, KE £ DRO,K24, be SHTESF ONS / DRIVER, 283 
Cc DRIVER, 284 
C TS = TOTAL NIMBER OF SPECIES VRIVER, 265 
C SMI(CT) = MASS NF Ney 2, fp AR, HE,» AND CU2, GRAMS URIVER, 286 
C ORITVER, 267 
OATA 18, CSMTCT eLEt,6) 4 be ByOSITF ord, SB 85SFe2S, P,HSHTEMSS,URIVER, PBA 

* Be ASSSF ers, FeHUAUF OD, 7, 50R0F a2 S/URIVER, 289 

C DRIVER, 290 
C ALPCT) = THERMAL HIFFUSION COEFFICTENT ORIVER, 291 
C DRIVER, 292 
DATA (ALPCT),Tet¢6) 4 UROede 90,49, 0,0 / DRIVER, 293 

C ORI VER, 294 
Cc DRIVER, 295 
Cc CNEFFICTENTS GeSUR@k FUR G/TM DRIVER, 2946 
t DRIVER,297 
OATA (AACT), 121,12) 7 be SISUSARPEFON, PeOSUSLUDSF HO, DRIVER, 298 

* 2eSULSH4 {HE OP, wr TUUTUU TOF eS, DRIVER 299 

* PeUSBATAAKHF aNd, Se S4OVMUIBSEROK, Ayal TAANGATF aN, DRIVER, 300 

+ 09,59552721 SF aly, Bek ASSSSHTE CLI, eA ~OAN0STAUF ej, URIVER, 301 

* UATIGDNIAPFOI7, Oe S/ ORT VER, 302 


172 
174 


17h 
177 
Pol 


Cr we 


C8 


6 # 
& a 


im) 
ao 


aon CeAaAAASAAMAAa aA 
an 
er 


a 


* 
* 


* 


* 
* 


* 
* 
* 


® ARITHMETIC STATEMENT? FUNCTIONS Tu CALCULATE 
* G/TM, INTEGRAL OF G/TM, AND Ge 
GPIMAF( a@ Y = CCECOOOOOE AdCL2)*eAG & AACL1))#AD © AA(1N)) AO 

+ BATG))#AQ © AA(8))#AQ & AACT)) MAD © BAL HI) RAD 

# AAC5))#AN & AA(G)) FAQ © AACE) MAD © AAC 2)) HAN © AATY) 
GIMTAF( AQ ) & CC CCOECEECOOCAACI2AI/1I2*AQ © AACITI/11,) #49 

@ AACLOI/IO, eA & AAC9D/9,) RAD + AALB)D/A,) WAU 

¢ BACTY/7,) #AQ @ AAC KH) /6,) #490 + AA(5)/5,) eA 

+ ARCUY/4O,) AD & AA(3)75,) AD & AAC2I/2,) AQ © AACS) I*AU 
GaFC RW ) & GZ/ 01, 0FRASRE) #42 


SFLAFC Au) 


x Aw 
ew oF 


GkK7AF( AQ) = 


6M FT 


RQ s 
cri 


GG os 
GAT 


GAT™ 


& x ARITHMETIC STATEMENT FUNCTION USED Tr) CALCULATE M/MSTAR DAY, 


FYPC COCCECOCOECOCDN (13) eR & DOC 12)) #bO ¢ AD(11) eR 
NNC10))#BQ + AHEC9)) #469 © HYVCB))#BD © NN(7)) eR 

DN (6))#RQ © UN(S5)) #8 ¢ ON(G) aR + DN(3))eRG 
NNC2))I#AW # ON(1) ) 


“+e 


ITHMETIC STATEMENT FUNCTION USED CALCULATE DENSITY SCALE 


IGHT (KM), 


TO 


COOCOCOEE AACY2I#11 8A © AACLI IRL, HAL 
+ AATLO) HQ V RAD © AACDIER,IEAG * AACR) HT, ) HAND 


+ AAC 7)%6,)wAD © AAlH)#E,)RAQ + AA(S)#U,)eAN 

+ AAC AYeZ, HAD © AAC EIRP, HAD © AAP) 
STATEMENTS 100 TO 20001 ARE DANE JUST UNCF, UN A CALL TC 
ATMNSU(T, 129), TO SET UP NEEDED PARAMETERS AND TU EVAL UATF 


SNL AR@FLUXeNEPENDENT FUIIRTER CUEFFICIENTS USED IN COMPUTING 


THE TIME wOFPENDENT VALUFS OF TAll, THE VARTARLFE 
TEMPER ATIIRE GRANITFNT AT THE LOWER BIUNDARY, TIF, THE 
EXUISPHERTC TEMPERATIIRE (SEE J, Se NISRET, RADTU SCIFENTE VOL, 


be P, 4t7 (1971)), AND THE CUFPFICIENTS IN THE PaRARULIC 
TRANSTYTON FUNCTION FOR THE DENSITY SCALE CHEIGHT RE TWEEN 
THE LOWe AND HIGHeALTITUDE MONELS, 
SUBSEQUENT CALLS, TO aTMOSUC2, 24), 
WHERE AFYER A LiiWeaL TT TuNE MOOFL IS 
LESS THAN 420 KM AND & HIGHOALTITUDE MODEL 
OU (109,209), JJ 


GO 1M STATEMENT 200 
USFD FOR aLTITIIDFS Z 
TS USER NTHERWISE, 


INTTIALTZATTON 


SKeHTGA 

= 1,9F ¢95*RIGMS/RR 
COMPUTE GRAV, ACCEL. 
INTEGRAL ‘IF G/TM AT 
GaF( 7H ) 

2 GOTMAFC 7H ) 

I 2 GYMYAF( ZMH ) 
COMPUTE PRESSURE, 
ATCORAING TO THE 
THE ANUNDARY CUNDTTIOUNS at 


™ 


G, G@ DIVIOFD BY MUL, 
120 KM, 


SCALE TEMP, » AND 


OND TEMPERATURE AT {20 KM 
MODEL, THESE VALIIEFS PROVTUF 
KM FOR THE HIGH@AaLTITUE 


DENSITY» 
LIweAL TITUNE 
12N 


90 


MODEL 


DRIVER, 303 
OPIVER, 304 
DRIVER, 305 
DRIVER, 3046 
ORIVER, 307 
DRIVER, 30R 
UVRIVER, 309 
URIVER, 319 
DRIVER. 3145 
DRIVER, 312 
DPIVER, 313 
ORIVER, Bu 
DRIVER, 315 
URIVER, 316 
URIVER, 317 
URIVER, 311A 
VRIVER, 319 
ORIVER, 320 
URIVER, 323 
DRIVER, 320 
URIVER, 323 
DRIVER, 324 
DRIVER, 325 
DRIVER, 326 
DRIVER, 427 
URIVER,32A 
DAIVER, 329 
DRIVER, 3309 
VRIVER, 331 
DRIVER, 3452 
DRIVER, 343 
DRIVER, 354 
DRIVER, $35 


CONTROLLING THEURIVER, 356 


OR LVER, 337 
DRIVER, 338 
DRIVER, 330 
DRIVER, 340 
DRIVER, Sut 
DRIVER, 342 
ODRIVER, 34% 
DRIVER, tag 
DRIVER, 34S 
ORI VER, tun 
DRIVER, 3u7 
URIVER, yur 
DRIVER, 3u9 
DRIVER, 350 
ORIVER, 35} 
UKIVER, 352 
DRIVER, 35% 
URI VER, 354 
ORIVER, 355 
ORIVER, 356 
DRIVER, 357 
DRIVER, 358 


203 
ais 


2e1 
P22 
20k 
sn 
34 
236 


Pul 
pus 
put 
259 
25? 
254 
271 
oTs 
300 
303 
wo? 
3409 
34 
417 
323 
325 


330 
333 


344 
335 
yan 
Bue 
sat 
was 


tus 
459 
45? 
355 


357 
vol 
tou 
te? 
Lea 
y74 
316 
4ol 
uot 


nag 


1a 


16 


PP = PLEEXP(el CI RGDIMT & FUL dur eNBa Liaw, h33) 
RHU cs HIGMS&GNTIMsRRRPPSGL, 
NOw CALL TRE & AUYILTASRY ROUTINES, 
CALL ZIT jUT 
COLL JULTANFTYRS, TMONS, THAYS,YREJ,VEuS,Davu) 
Cat Sil CYC(DAYA) 
CALL SOLURBCYRE J, VEQI,NAYJ,S0LL4T,SOLI UN) 
CAL! SUL ZENOSOLL AT, SUI LON) 
CALL SPCMIN( 1574) 


COMPUTE M/MSTAR (NEPENUS ON DAY US NIGHT PROFILE OF 
APTER CUMPLITING 


PRINCTPLE, RUT WE COMPUTE YT MANLY FOR LAY) 
COEFFICTENTS DD(I), 

RAPRR = P,eRATGASKR 

D9 108 Nzy,NALIND 

w7Z = al TKM(N) 

wG6G s GAFC w2ZZ ) 

wGDTM Bw GATMAF( w2Z ) 


wPP = P7REXP (CCL ®GTMTAF( w7Z ) & GUI UUEwURaWwZ7 eee KY3)} 


MMARHM = WORST KACRKMWPRRRGUTM) 

FONAY(N) Ss AMARHI FOND AYV(N) 

CONTINUE 

DMN 106 131,13 

DN(T) s 0,0 

CONTINUE 

CALL FITTERCNAL TON, ALTKY,FUAYeNDEGy 1 » 2 » DN) 
USING DAY () PRUFILE 

SF 2 SFNAF( ZH ) 

OMBMS = 1,07¢ 1, © SF ) 

$Z zs RMAMSe&GG/GDTM 


COMPUTE THe SPECIFS NUMRER NENSITTES aT 120 aM, 


COMPUTE TOTAL NUMRAER DENSTIY,N(S/OMe83) 
SN = RIGA/SRIGMS#@RHUSRYRMS 


VEIT VER, $59 
DRIVER, 360 
DRIVER, 3o! 
DRIVER, 302 
DRIVER, Se3 
ORIVER, You 
ORIVER, 36S 
VRIVER, $66 
URIVER, 367 
URIVER, 368 
DRIVER, 69 
URIVER, 370 
URIVER, 371 
DRIVER, 372 
DRIVER, 37% 
ORIVER, S74 
ORIVER, 378 
ORI VER, 376 
URI VER, 377 
DRIVER, 37A 
DRIVER, 379 
ORT VER, 380 
DRIVER, At 
URIVER, 382 
DRIVER, SB 
DRIVER, 3K4 
ORIVER, BRS 
ORI VER, 3B6 
OPIVER, 3a 
URIVER, 3BA 
ORIVER,3K9 
ORI] VER, 399 


COMPUTE TOTAL NIIMBER DENSITY TF NO DISSOCTATION,NSTARC{/CM@eSIORIVER, 394 


SNS = KIGARRHO/RIGHS 


COMPUTE DENSITIES (1/CMae3) UF N&, M2, ND, ARy 


SNI7(1) = 0,7828N8 
SNIZ(2) 2 4,211*SNS @ SN 
SNIZ(3) = 2, @SNS*SF 
SNIZ7(4) = 0,0098#SNS 
S®17(5) 8 4,625Fe05*SNS 
SNIZ(4) = C€N2025) 


RE{P0 2 RFO120, 
GESK = GGsSK 
CC ws OLeHl 442, 


FF = SBAR 
COMPUTE FOURIFR CNEFFICTENTS USED FOR TAU AT 

AC({) & #2,2101SAEPU2 © 1,97N080F ENS * FF 

AC2) = 6, 71235RF OOS © 1,1B11N7TFeNS * FF 

AC3) B 2, 7UB1RNECOU © 3,39NS22F en? * FF 

AfUY B@ @S,6b5S477E ROU + A, 66901 HE en? * FF 

A(5) © #U,65275RAE S08 + 2,327930FeNT7T *& FF 

A(6) © FB, GRUS5SUF ODS © I,1CAIS7TE SN? *& FF 

BLY 2 e$, UNT39RAE SOS + 1,90N959F eos # FF 

RO!) 2 eS, UPKS97E COU © ULL0F51 SF eno # FF 

HOS) = e2,SIRIRGECOK w@ S,3uspigken? # FF 


91 


URI VER, 392 
ORIVER, 39% 
DRIVER, 394 
ORI VER, 395 
DRIVER, 396 
ORIVER, 397 
DRIVER, 39R 
VRIVER, 399 
DRIVER, GUN 
ORIVER,401 
DRIVER, Gu 
DRIVER, YOY 
ORI VER, 4OUu 
VDRIVER 40S 
DRIVER, W046 
DRIVER,UOT 
DRIVER, UOR 
UDRIVER, 409 
DRIVER,4IN 
VRIVER, 414 
DRIVER,412 
ORI VFR, UL3 
DRIVER, aja 


406 
440 


ays 
418 
420 
ep 
42s 
4e7 
&%3p 
4ya 
437 
Guy 
4“gu 


Gu6 
ago 
us? 
45% 
4s4 
as? 
Gas 
472 


AHORA AAaAaeananranoannan 


fon) 


110 


HOU) FS 1, 3ROKUSE COU «6 O,97TS3AGE NT & FF 

BOS) =F +], 55K99UFe OU + 3,951R11F wo? * OF 
COMPUTE FOURPTER COEFFICTENTS l1SED FOR TIF, 

COL) F 45,4USSZRE SOP? © 4, S2RBOTFENU *& FF 

C2) = Hp ,179RY9F 402 @ H,U9S SHOR en *& FE 

COS) S HF,115091E 401 @ 4, 7OHRIBE WN & FF 

COG) = +4, OAOSQEE HON & U,15N6APFene * FF 

COSY © s6e3ROQOKIF 400 © 1 41S THOF Rene * FF 

COO) & +1, 0454R7E400 © 1,995652F ene * FF 

SOL) & wtel386OtEtul @ T,29R7UGF ent & FF 

SC2) & 41. $59KHRE4O1 & 2,815 729F NS & FF 

SO3) 3 #9, R5915AREMO1 & RLISRBARLEwNe & FF 

STUY = +7,OALI3PF ONL © T,1S17TOREMNe & FF 

SOS) & 2.92551 SEMUL © 4, 6252 tor ene & FF 
COMPUTE TAL (1/4KM) AND TIF (DEGREES KELVIN) 

TAU s A(1) 

TTF S.4089 

Ph EtG M145 

BY = -Y 

SFI = SIN(CCRFY) 

GFT = COSCOCaFY) 

TAU & TAS © CFT Halley) ¢ SkTRACT) 

ThE 2 FTF « CPC Ley ©& SF resitth) 


SA 2 


TO PROVIDE & CONTTNUDUS DENSITY SCALE HEIGHT ACROSS THE 


BOUNDARY RETWFEN THE LUWe AND HIGHeALTITUNE MODFLS, 


wk 


PARAL TC TRANSITTON FUNCTION, 


HRA) = FHRI20 ® ZHM1{0n02 + SA & ZHMI10 * HROTIO 
WHERE 
HRV{10 = NENSTTY SCALF HETGHT AY {10 KM 
ZH¥M110 = 7Helln, 
84 = APPROXTMATE MERTVATIVE LiF DENSTTY SCALE HETGHT 
AT 1100ekKM ALTI TUNE 
= HRI1L05enR 1095 
MR11N5 = MENSTTY SCALE HETGHT AT 110,5 KM, 
HP1N95 & NMENSTTY SCALE HEIGHT AY 109,85 KM, 
FHRI20 2 (HRUIL20 © 10,8SHK = NRULIN)S(120,H110,) ea? 
IN THES INTTTALTIZATIUN CALL WE NEFLD) TO COMPUTE THE RENSTTY 
SCALE HEIGHT AT 120 KM, HRUIT20¢ ACCORDING TO THE HIGHeALTITIIDE 


MODFL, WHICH DEPENDS ON HL AND SAR, AND ALSO THE DFENSTTY 
SCALE HEIGHTS ACCORDING TO THE LUwWeALTITUNE MOOFL At 
110,5 «4, aNp {09,5 wm, 


COMPUTE SMALL Ay 


COMPUTE COEFFICIENT 


¢ 


GSM’ zs 
3 0,9 
NPUNIN £ AYN 


RHUL 


bn 


TEP =e TLYAT KE 


QF MaSiiRetl IN GAMMA®PSIIBS] 


1, OF e0SahGSK/ (TIF a TAt) 


120 181,15 


SNZSMT @& GNIZ(1)#S5MT(7) 
GAM 3 GAMTaeSmT(T) 
Al GAM 


wh) 


NRORIN 
CONTINUE 
HRLIY 20 


s 4LP(T) 
RH) # SNZSMT 
e& PRADIN + 


# GAM 


+ 120 


SNZSmMTel(GAmM ¢ ALGAM{H#SA/(1,0S4)) 


= RHO/SNRANIN/TAU 


92 


USF A 


qo KM, 


DRIVER, 43S 
DRIVER ,UiA 
URI VER 417 
DRIVER,U1R 
DRIVER, 419 
URI VER, 420 
URI ver, U2) 
DRIVER, 42? 
URI VER, 423 
URIVER,deu 
VRIVER, 425 
URI VER, 426 
ORIVER,427 
ORIVER,4OR 
DRIVER, 429 
DRIVER, 430 
ORIVER, 431 
DRIVER, 432 
ORIVER, 433 
DRIVER, 434 
DRIVER, 435 
ORIVER US 
DRIVER, 437 
DRIVER, 435A 
DRIVER, 459 
URI VER, Gun 
DRIVER, ut 
DRIVER, 442 
DRIVER, 443 
DRIVER, Guu 
DRIVER, das 
DRIVER, 4u6 
URIVER, Gu? 
DRIVER, GUA 
DRIVER, 449 
DRIVER, 459 
URIVER, 451 
URTIVER, 452 
DRIVER, 453 
OP IVER, 454 
ORL VER, USS 
DRIVER, USS 
DRIVER, 457 
ORIVER,USR 
DPI] VER, 459 
DRIVER, 460 
ORIVER, 461 
ORI VER, 46? 
DRIVER, G63 
ORI VER, Ubu 
ORIVER 46S 
DRIVER, 466 
ORI VER, 467 
DRIVER, GOR 
DRIVER, 469 
URIVER,U70 


S34 
Su? 


556 
557 


S74 
S77 


A165 
hen 
eu 


534 
AVA 
640 
AU 
6u7 
659 
45! 


AS? 
454 
&SA 
hat 
47? 


ATA 


CAARANANNMALCMAGANDANMANANAO 


Rony 


1%0 


* 


* 


* 


COMPUTE DENSITY SCALE HEIGHT AT 140 K™, 
GNTM s GNYMAFC 419, ) 
MRIIVEN = 1 ,OSCCOYREDTM & 2, HH7TOIISCE MNT HL Ig OHR1 ASS 
© A2,0/(RE+1190,0) » GKKZAFC 110,09 )/469TM) 
COMPIITE DENSITY SCALE HFIGHT aT 110,5 KM, 
GNIM s GOTMAFC 149,5 ) 
HPL1OS Ss L,A/(CCL#GDTM © AP ~hETO9IS2F ANT HEL10,5 eet R33 
@ 2, 0/(RFO110,5) = GKKZAF( 110,5 )/GDTM) 
COMPUTE NFNSITY SCALE HEIGHT af 109,5 KM, 
GAIM = GOTMAFE 109.5 4 
MRINIS S L,O/(CCL#GDIM = A hO7TU995SZEMOTHINGI, Sees, 835 
= 2,0/7RF 4109.5) = GeKZAF( 199,5 )/GNTM) 
SR = HRIYNSeHRINGS 
FHRIQN S AOL *(HRALAQ @ YO,eSR © HKINIO) 
WRITE CH,HOO1) TIF, TAU 
FOARMATC/S/,7H TIF & sEi3y607H TAll = ,E13,6;/7) 


AT NIGHTTIME, 0 DIFFERS FROM DAYTIME OM ONLY BELOW ALTITUDE 
ZTUN(S) © ON KMy IFC ZH LT eZTUN(S)), WHERE ZTON(1) = 60 KM, 
SNECS) = TNZIC1) = ONITECAS) 32 4,1 « 

IFC PH GEZTONC YS) GAND, ZHeLT,ZION(2)), WHERE ZION(C2) & 75 KM, 
SNIC3) = ON7T CQ) aFXP(ZMAINSONSCHL) WHERE 

INZT(2) = UNITE(16) s 4,9NE+0A 

ZM20N & ZHeZTIONN( 2) 

NONSCHT = ALUGCONZIC2I/0N7ICLIIAC7IONC2@71ONG1)) 
IF CZH.GT,ZIANC2) ZANDe PHal Te ZINC 5)) WHERE ZINNCS) & 90 KM, 
SNICSY & tageeC COOK CL) eZM20N © KC2YIHZM2ON © K03)) eZMAEN 

# X(4))eZMPUN © K(S)IHZMPUN © ¥(K)) 

WHERE THE CONSTANTS X(1),Te1,6 SRE NETERMINED SO THAT THE 
SLUPE UF SLOGIOCSNE(3)) af ZINNC!Y) = 7S KMp OLUZe7, AND AT 
ZTONCS) = 90 KM, OLOZ5Z, TS5 CONTINUOUS AND ALOGIO(SNIC3)) 
EQUALS THE NIGHTTIME VALUES FORK U aT Z7I0NC2) = 75 , 
ZTON(3) & AN , ANA ZIONC4) & BS KM AND EQualS THE DAYTIME 
VALUE FOR 1 AT Z7IANC5) & 90 KM, ALOGIO(IDAYZS), 
THE NIGHTTIME () CONSTANTS ARE NOW SET, 

ZTON(1) & ALTKM(13) 

INZTCT) @ ONTTE CY 3) 

DN 1350 Ts2,5 

ZTUNCT) & ALTKM( TOYA) 

HNZTCT) s&s ONI TEC] 444) 

CONTINGE 

Zee = ZIIN(S§) 

TN RESET ON7105) TU 17S PROPER VALUF WE NEEM TO FIRST 
CALCULATE Ma ES ey 


COMPUTE GRAV, ACCFLe Ge G OT VIDEL KY MUL, SCALE TEMP, TM, 
AND INTEGRAL JF GTM AT ALTITIVE ZHA, 


6G « GAFY 7a ) 
GATM & GITMARE FHA J 
GNT4] & GYMTAF( ZHe ) 

COMPUTE PRESSURF ant DENSTTY AT ALTYTTUDE ZH2 
PP = PJMEYP (al CIMGDIMT © D,Up, UF enbeZHeee?, R33) 
RH 3 RT GMS@eGNTMsRREPD GG 

COMPUTE M/MSTAR DAY at ALTITUDE 2He 
SF 2 SFNAF( ZHo ) 


93 


DRIvea 
ORI VFR 
DRIVER 
DRI vEeR 
URIVER 
DRIVER 
URI vE® 
ORIVER 
ORIVER 
ORIVES 
DRIVER 
ORI VER 
DRIVER 
URIVER 
DRIVER 
DRIVER 
DRIVER 
DRiveR 
URI VER 
DRIVER 
DRiveR 
ORI vER 
VRI vie 
DRI EPR 
ORIVER 
URIVER 
DRIVER 
DRIVER 
URI VER 
ORI VER 
DRIVER 
ORI] VER 
OR[ ver 
ORI VER 
QRIVER 
ORIVER 
DRIVER 
DRIVER 
ORI VER 
ORIVER 
URIVER 
DRIVER 
DRIVER 
DRIVER 
URIVER 
DRIVER 
URI VER 
URI VER 
VRIVER 
OR] VER 
ORTVER 
DRIVER 
DRIVER 
ORI vER 
URIVER 
URL VEs 


703 BMAMS 3 1,0/¢1, + SF) ORI VER,& 
C COMPUTE TOTAL NUMRES DENSITY, NOV/CMe#S) AT ALTITIIOF 7H? ORIVER,* 

705 Sof RIGA/BTGYS*RHU/SBMHMS ORT VER,¢& 
C COMPUTE TOTAL NUMRER DENSTTY TR NO NISSOCTATION, DRIVER,¢ 

C NSTAR (1/CM*e3) ORIVER,¢ 

749 SNS = BIGARRHIYAIGMS DRIVER,§! 
741 (INAYZS = 2, *SNQRSF DRIver,* 
trys WUNZTCS) = ONAYZS DRIVER,S 
14 ANSCHT s ALNGCONZTC2)/0NZT (1) 9/0ZTUN( 2) ©2TUN(4)) DRIVER,§* 
723 OLU727 = ONSCHYMALUGIO( FXP(1,0) 3 DRIVER,* 
730 x(5) £ DOLNZAzZ DRIVER,¢ 
731 Kfod = ALNGIOCONZT(2)) DRIVER,* 
735 oo 435 733,85 URIVER,§® 
748 ZIMAUN(T) = ZINNCTI@ZTUNC2) ORIVER,S 
24S 13S CONTINUE ORIVER,S 
Pub nm 4an Tet,3 DRIVER,® 
woh 2112 = ZIM2NNCT +2) DRIVER,§ 
152 OCI,4) = Z21%2"7TT2 DRIVER,' 
75% on 140 J=4,% DRIVER,* 
757 OCl,GeJ) =f 7ITARNCL Sed) DRIVER,§ 
765 140 CANTINUF DRIVER,§ 
770 2718 & ZIMANN(S) URIVER,® 
74 OfOed) = P,*ZTIS DRIVER, § 
773 oN 145 tat,3 URIVER,§ 
101 Pbk al Set DRIVER,§ 
10% NOG,GeJ) = FITSa( CRI OL, SF IG eNO (4, 5eJ) DRiVER,¢ 
1945 145 CONTINUE ORIvVER,§ 
1044 mM 4150 Yat,3 DRIVER,S 
1024 O(T,5) & ALNGIOCONZT C1429) @ KCS)*ZIMAUN(T #2) @& x6) DRIVER,® 
1031 150 CONTINUF URTVER,§ 
C TO SET N(N,5) WE NEED THE DERTVATITVFE OF ALUGIOCSNT(3)) AT ORIVER,§* 

‘a ALTITUDE ZINN(S) #& 9H KM, DLUI57, GIVEN BY URI] VER,* 

C PLOZ5Z2 2 ALOGIOE ExPC1,9) JECNCSFIZNZ © 1, 0/HRHN), ORIVER,€ 

C EVALUATED aY 759N = ZIUN(S) & 90 KM, DRIVER, € 

C COMPUTE DENSITY SALE MFIGHT AT 90 KM, DRIVER,S 

103% GOT = GITMAFC 90, ) ORI veER,S 
134 HPII9OD S LeNZITCL*AGOTM w A,OHTNIVHAE MO TRIV, AH] HSS URIVER,S 
* © 2,07(RE+90,) = GKKZAFL 90, I/GOTM) @ DRIvER,§® 

105% 250" = ZINNES) ORIVER,§ 
1954 MN7ST & ALNGIO® FKPECT OD DECOCOEECOCE CHS, tOO C13) eZ50N ORI VER,§ 
* * UPON CLA) OFSIN © LO HONCT 12) Z5ON & 9,8DN(10))4e750N DRIVER,S 

* + B,*DN69)Y*Z5NN + 7LeNN( KI) *Z5UN * 6, #0907) 18759N URIVER,§® 

+ ¢ S,#uUN (6) V*Z75ON © U,eDDE5)I*Z50N © $4 81D0(04) 2 RZ5NN DRIVER,$ 

* $ 2.eNN(3)Y*#ZSON © ONC2)) = Ler sHRNDD DRIVER,* 

Li by Ut4e5) 3 DLIZSZ=xt5) DRIVER, ! 
1120 NO = 4 DRIVER,' 
112? CaLk Sul vECN,©,Nt) VRIVER,§ 
c TA PRAVTDE & CUNTINIIGUS TRANSTTIUN IN THE Clie DENSITY BETWEEN DRIVER,' 

C THE ALTYTIINE OF 190 KM, RELIIw aHIC 4 A CONSTANT MIYING RATT OP IVER,! 

C TS ASSUMED, AND THE ALTTTIUDF MF 130 KM, AT WRICH THE aTMoSU URIVvVER,! 

C HTGHeALTITUNE MODEL (RASED ON OTFFUSIVE EQUILTBRIIIM) REGINS, ORI] vER,! 

C wE USE THF POL YNOMY AL VRIVER,! 

t LOGIOFSNI fod) = SUME ¥OCT)eZMICNerel(7Fel)), Ter,? DR1IvER,! 

C WHERE THE CONSTANTS xC(T),1e1,7 » ARE DETERMINED SU SHAT THF ORT VER,! 

C SLUPE HF ALMGIOCSNI(6)) at ZICU2(1) & 100 KM, CLE 712, AND DRIVER,! 

C At 7ICOPCS) = 126 KM, OLEZSZ, IS CONTINUOIIS AND ALIIRLOCSNT (46) ORIVER,! 


D-A0%e 013 GENERAL RESEARCH CORP SANTA BARBARA CALIF F/6 4/1 
THE ROSCOE MANUAL+ VOLUME 14A. AMBIENT ATMOSPHERE (MAJOR AND MI==ETC(U) 
JUN 75 D A HAMLINe M R SCHOONOVER DNA001~74=C=0162 


UNCLASSIFIED DNA=3964F=14A 


. FEB 
a le | 


1124 
1133 
1134 
1135 


1136 
1149 
1142 


L144 
1157 


1162 
11aa 
1166 


1170 


1210 
1216 
1229 
1?26 
1230 
131 
1734 
1235 
1236 
124? 
1259 
1253 
1254 
1254 
1264 
1264 
1309 
1301 
1304 
1344 
1314 
1332 
1333 
1335 


OOM 


no 


NANaD 


an 


140 


lag 


1465 


170 


175 


ENUALS THE VALUES FOR CNM2 AT FICUA(T) © 19V,195,110,115s A 
$20 KM FOR T21,5 , 
THE CN2 CANSTANTS ARE NOw SET yee 

ON too T31,5 

ZTCNACTY & ALTKM(T #20) 

CNP7ICI) & COACT+20) 

CONTINUE 
RESET COe@71(1) TO THE VaLii€é ORTAINEN 
MODEL AT ALTITUNE ZTIcHet1) = 100 KM, 


FROM THE LOweal TITUNDE 
TU DD THIS wE MUST FI 


COMPUTE GRAV, ACCFL, Ge G DIVTOFD KY MUL, SCALE TEMP, TM, 
INTEGRAL MF G/TM AT 100 KM, 
COMPUTE GRAV, ACCFLe Gs G DIVIDED BY MUL, SCALE TEMP, 1H, 
INTEGRAL OF G/TM aT 100 KM 

GG = GAF( 190, ) 

GNTM 3s GUTMAFC 100, ) 

GNTM] s GTMTAF( 100, ) 
COMPUTE PRESSURF AND DENSITY at 190 KM 

PR = PZREXP(eCCIRGDIMT & 9,4, duFen8#100, 82,833) 

RHO 3s RIGMS&GDIM/RREPP/GG 
COMPUTE TOTAL NUMRER DENSTTY TR NT NDISSNCTATION, 
NSTAR, AT 190 KM, 

SNS = BIGASPHI/RIGMS 

CN27TC1) & 3,70FeNG & SNS 

xP OF) S ALUGIO(CHAP7T(1)) 
THE SLOPE VF ALOGTOCSNI(o)) AT ALTITUDE ZICNeC1) = 100 KM, 
OLC717, 18 GIVEN Ry DLE71Z = ALOGIO(EXKP(1,0))e( 01, /RAN) 


REDERHGIZNZI) 3B ALIGINCE XP (19d) 8 wi, /HRHD), 
COMPUTE DENSITY SCALE HFIGMT AT 10y KM, 

RRUL09 3 TLASECOIRGDTM © 2, h670995eF OOP e190, ee] 635 

@ A,OsERES100,) © GKKZAFC 100, )/UNTM) 

DLC71IZ & Cas AsHRALNOYRALNGIO( FxP(1,9) ) 

xe (A) 3 OLC712Z 

OM {ou Ts?,5 

ZIMICCI) 3 7ITHACTIEZICNE1) 

CONTINUE 

On 168 21,4 

ZT12 sw ZIMIC(T 41) 

OCI1,5) = 721 %2e7T12 

ON 165 Jst,4 

OCL,SeJ) 2 7] T2eO fT e6eJ) 

CONTINUE 

ZTIS #& 7IM{C(S) 

(5,5) 2 2,*7T1S 

On 170 Js1,4 

FJl = Jel 

0(5,5eJ)) 

CONTINUE 

oN 175 Ts1,4 

OCle6) & ALNGIOICNSZIC1O1)) © xO (Ce)eZIMICfI +1) © ¥CC7) 

CONTINGE 

DLC757 = ALNGIO‘ FxP (1,0) 

005,60) = DLE Z7S7exC (4s) 

NO 8s § 

C4LL SOI VEC(N,¥C,N) 
COMPUTE N2 AENSTTY AT 2800eK™ ALTITUNE FOR USE IN NeNENSITY 
INITIALTZATITON Th SURRUNITINE SPCMIN, 


= 71TS# (CPI t ote d/FI1 #005, 603) 


DETAY#(SASMICOVEGAMT IS (SA01, 0) 


Nn 


RST 
AND 


ANY) 


OR] VER,5&3 
DRIVER,SA4 
ORI VER, 585 
DRIVER,5S86 
DRIVER,587 
DRIVER ,S6A 
DRIVER, S69 
ORIVER,S90 
ORIVER,S91 
ORI VER,592 
ORIVER,S93% 
ORI VER, 594 
ORI] VER,S9S 
DRIVER, S96 
ORIVER,S97 
VRIVER,59R 
DRIVER,S99 
DRIVER, 600 
URTVER,601 
DRIVER, 602 
URIVER,6043 
DRIVER, 604 
DRIVER, 60S 
URIVER, 406 
DRIVER, 607 
DRIVER,608 
DRIVER, 4609 
DRIVER, 610 
ORIVER, ALT 
VRIVER,612 
DRIVER A143 
DRIVER, 14 
DRIVER, 415 
ORIVER, 616 
DRIVER, 417 
ORIVER 618 
ORI VER,619 
ORIVER,620 
DRIVER, 621 
DRIVER, 622 
ORIVER, 623 
DRIVE® 624 
DRIVER, 625 
ORIVER, #26 
ORIVER 627 
ORIVER, 462A 
ORIVER,629 
DRIVER, 630 
ORIVER,A31 
DRIVER, 432 
DRIVER, 633 
ORIVER, 634 
ORIVER, 435 
URTVER 636 
ORIVER,A37 
DRIVER, 635A 


1337 
1343 
1351 
1355 
1357 
136! 


Cert 
1372 
1374 
1375 
1375 


1376 
1402 


140? 
1404 


1404 
1404 


1407 
1440 
1440 
1442 


1445 
1416 
1420 
1422 
1426 
1434 


1442 
1446 


1460 
14627 
1463 


On 


£77 


Cccee 
Cc 
€ 
Cc 


c 


2t0 


Z7 = RELA9ACALTKM( UT) e120) SCRE FAL TKMCUT7)) 

EYTZ = EYPCeTA'I#77Z) 

TTOTZ & ( TIF SCTIFeTZ)eF TZ J/TZ 

GAM = GAMT#SMI(1) 

ALGAM1 B ALP(1)4GAM41,0 

S1ZAPN = SNIZ(N)*ETZERGAM/TINTZ*&#AL GAM] 
EVALUATF ATMUSPHERTC PROPERTIES AY 90eKM ALTITUNE PRINR 
TO INITTALIZING IANNSU, 

ZHSAVE = 7H 


ZH s 90, 
JiMP sn 
GN Tr) 210 
Jime = 2 
INITIAL TZE TONOStU ROUTINE, 
CALL ITONQSUCL,Z4) 
24 = JHSAVE 
SFT ZHFLAG (ARRITRARY NEGATIVE VALUE) 
ZHFLAG = e2n, 
RETURN 
CONTINUE 
IFC ZH,EQ,ZHFLAG 2 RETURN 
AN FRRONENUS CONDTTTON WILL UCCUR TF IONOSU UR SPCMIN JIS 
CALLEN WITH JJ=2 aNnD A GIVEN VALUE NF ZH YF ATMNSU HAS NOT 
BFEN CALLED FIRST wITH Jyz2 AND FOR THE SAMF VALUE OF ZH, 
THE VARTARLF 7HFLAG IS UISFbD TO DETECT THIS CUNDITION AND 
TO MAKE THE REQUIRED CALL TY AlMQSu, 
ZHFL AG TS INITIALIZED TN AN ARBITRARY NEGATIVE VALUF IN 
THE INET TIALTZATION Caty TO ATMUSU, 
ZHFLAG = 7H 
CONTINUE 
REZHI] & 1,0/( REOTH ) 
IFC ZH ,GE, 120. ) GA TO 250 
(s LOWeALTTTIUDFE MODEL (ZH ALT, 1200) 


COMPUTE GRAV, ACCEL. AT ALTITINVE 72H, 
GaAF( ZW} 
COMPUTE GRAV, ACCFL, 
GONTM 3 GOTMAFC 7H ) 
COMPUTE INTEGRAL MF 
GNTMT = GYMTAF( ZH ) 
COMPUTF FUNCTION NEEDED FOR MENSTTY SCALE HETGHT 
GKKZ & GKK7ZAF( 7H ) 
COMPUTE PRESSURE (DYNES/CMRH#O) 
PP = PZREXP(eCC1EGNTMT + 9,4, UUF wehbe ZHeee 633) 
COMPUTE DENSITY (G/EM*eet) 
RHI) = BIGMS#GDTM/RREPP/GG 
COMPUTE DFNSITY STALE HEIGHT (KM), 
IF (24 GF, 110.) GY TH 230 
HPHO = Y,O/(COLHGDTM © 2, 6670995 CF E07 HZHeeL 83S = 2,0 HRKF ZHI 
* © GKKZ/GNT™) 
69 TO 24) 
ZHH110 © 7H © 110, 
HOH £& (FHRION&7HMII0 


GG(CM/SEC eRe), 
b6 = 


DIVIDED AY MOLECULAR@SCALE TEMPERATURE, 


G/M, 


+ S8)eZ7HMI10 * HROLI0 


DRIVER,&39 
ORIVER Aun 
DRIVER, 641 
ORIVER 642 
DRIVER, AUS 
URIVER, 644 
DRIVER, 64S 
DRIVER, 646 
DRIVER ,447 
ORIVER, 648 
DRIVER, 649 
ORIVER,659 
DRIVER, ASS 
DRIVER, 652 
URIVER, 653 
ORIVER, 654 
VRIVER,655 
DPIVER, 456 
ORIVER,657 
URIVER, 658 
ORIVER,659 
DRIVER, 660 
ORIVER, 661 
DRIVER, 662 
ORIVER, 663 
URIVER, 664 
DRIVER, 665 
DRIVER, 666 
ORIVER, 667 
ORIVER, 668 
DRIVER, 66° 
DRIVER, 670 
ORIVER, 471 
URIVER,472 
DRIVER, 673 
DRIVER, 674 
ORIVER 675 
DRIVER 674 
ORIVER,677 
ORIVER,678 
ORIVER, 479 
ORIVER, 689 
URIVER, 681 
ORIVER, 682 
VRIVER, 4683 
DRIVER, 684 
ORIVER, 665 
ORIVER, 486 
DRIVER,AB7 
ORIVER &BR 
ORIVER, 489 
DRIVER 699 
DRIVER, 491 
ORIVER, 692 
DRIVER, 693 
URIVER 694 


1467 
1472 


1474 


1477 
1501 
1802 
1504 
1506 
1519 


151? 
1814 
1517 
1829 


15934 
1535 
1540 
1542 
1547 
1551 
155% 
1555 
1556 
1562 
1563 


1601 
140? 


1604 
1613 
1615 


1615 
1615 


1420 
1625 


16359 
1639 
1431 
143? 
1433 


Nan 


Ann 


C 


ceccce 


C 
C 


us 

eug SF = 8 
BMBMS 

cn 

EY <8 


COMPUTE NUMBER DENSTTIES OF SPECIES, 


DA 


ING DAY PROFILE 

FNAF( 7H ) 

= 1.07(". + SF) 

MPUTE TEMPERATURE (DEG K) 

MR49#6G/G0IM 

wE PRESCRIBE THE 
YeNIGHT NFRPENNENCE OF OG AND USE THE LOweALTITUDE MODEL TU 


COMPUTE THE ASSOCTATED SLIGHT DAYeNIGHT DEPENDENCE OF 0? , 


SNS = 
SN e& § 
SNI(4) 
SNI (2) 
SNI(3) 
IF( IN 
en 
IF( éH 
IFC ZH 
Z™20N 
SNI(3) 


GO TU 


ZM20N 
SNI¢3) 
GN TO 
SNIC3) 
SNI¢4) 
SNIC5) 
IF( Zw 
zMiCoe 
SNI (6) 
* 

Gf Tf 
e246 SN] fb) 
qn 
FE 


we 

247 FFHSEQ 
TFC Ju 
RETURN 


HT 


COMPUTE THE GEGPUOTENTIAL ALTITUDE ARUVE 120 KM, 


250 CONTIN 
77 2 R 


COMPUTE 


ETZ = 
yY¥ & ¥ 


COMPUTE RaTT) OF TEMPERATURE TU 
2 ETS 

A) 

069 

29,0 

2 4,9 


TTOTZ 
PP z 0 
WH) os 
HRPUDZN 
DPPDTH 


IF( fHeZION(1) 9 


HIGA®RHN/BIGMS 
NS/RMAMS 

= 0,7&*#SNS 

2 1,.211*SNS e SN 

= 2,*SNS*SF 
URN,GF,O0 2 GO TO 245 
MPUTE NIGHTTIME VALUE OF O 
eGEeZS5UN Y GY TO 245 
s7IONf2) 2 242,242,244 
zs JReZIONC2) 

= 10,8 (0( CCX (9)*#Z429N © X(2))#7M20N & X03)) #72M20N 

¢ ¥(U)Ve7MAUN & X(5))4ZM20N # X(A)) 
24s 
244,245,243 
as ZHe7ION(2) 

© INF TOP RE KPC IMAONRNNSCHI) 
24s 

= ONZI01) 

= 0,009*SNS 

S U.62°5Fe0S*#Sn§ 
sles (00. } GO ToO72e4 

= ZHezICne(t) 

= 10,ee (CC OCOxE (1) ecyiCo2e + xCf2))e2MiCue + ¥CC3))#Z41C02 
* XC (ud #241002 & XCF5))#7241C02 + ¥CC6))*ZMICU2 * KCC7)) 
aur 

= 3,20Fen4u * SNS 
MPUTE FRACTTUNAlL FRRUR FROM HYQORUSTATIC EQUILIBRIUM, g. 
HSEQ = 91, 0F e05@DPPNZH/(RHIHGG) © 1,0 

EB e2,6A709952PF OL? & RR & ZHew},833 / (SIGHS * GDOTM) 

ERE 2ebATOIISPE@{S 2 1,0EeVS & 2,853 © 9,4, UF eB 

2 e2,667099S2F eA * RR w ZH**1,8335 / (BIGMS * GDTM) 
MP,EQ.0 3 GO TA 177 


GHeALTITIINE MODEL (26 ,GE, 120,) 

Z7Z(KM), 
UE 

EL 20e(ZHel 20, aREZHI 


EXP (CeTAlia22) 
IF wm (TIF STZ) *ETZ 


TEMPERATURE AT 120 KM, 


ORIVER, 695 
DRIVER, 696 
DRIVER, 697 
DRIVER, 4698 
OCRIVER.699 
ORIVER,700 
ORIVER, 701 
DRIVER,702 
ORIVER,703 
ORI VER, 704 
DRIVER,705 
ORIVER,706 
DRIVER,707 
ORIVER,708 
ORIVER, 709 
ORIVER,719 
DRIVER,T11 
DRIVER, 712 
DRIVER, 713 
ORI VER, 7144 
DRIVER, 715 
ORIVER, 716 
ORI VER, 717 
URIVER, 718A 
ORIVER, 719 
ORIVER,720 
ORIVER, 721 
ORIVER, 722 
DRIVER, 723 
ORTIVER, 724 
ORIVER,725 
DRIVER, 726 
DRIVER,727 
DRIVER,72A 
ORIVER,729 
DRIVER, 730 
ORIVER,731 
DRIVER, 732 
ORIVER, 733 
DRIVER, 7346 
ORIVER,735 
ORIVER,736 
DRIVER, 737 
DRIVER, 748 
ORIVER,739 
ORIVER, 740 
DRIVER,741 


THE TEMPERATURE AT THF GENPOTENTIAL ALTITUDE, TICOEG KIORIVER, 742 


ORIVER, 743 
DRIVER, 744 
ORIVER,745 
DRIVER, 746 
URIVER, 747 
ORT VER, 74R 
URIVER, 749 
URIVER,750 


1634 


14406 
lou? 


1651 
1461 
1603 


146604 
1671 


1673 
1675 
1707 
1744 
1716 
1722 


1724 
1724 


260 


NN 260 T2#t,TS 
COMPUTE GAMMAeSURe], 
GAM = GAMTaSMT(T) 
ALGOM{ 3 ALP(T) # GAM @ 1,0 
COMPUTE DENSITIES (1/CMen3) OF Ne, M2, O, AR, HF, AND C2, 
SNICLY © SNIZCT)®ETZa#GAM 4 TIDTZ*e®AL GAM 
COMPUTE TOTAL NUMRER DENSITY (1 /UMaed), 
PP s PP * SNICT) 
COMPUTE TOTAL MASS DENSTTY (G/CM883), 
Qui) = RH) @ SMTCTY®SNI(T) 
COMPUTE A PORTION OF THE SPATTAL DERIVATIVE UF THE NENSTTY, 
SGAFT = SNICT) a (GAM @ ALGAMURFTZe (TTF eTZ2)/TT) 
DRUNZN B&B ARNHIN + SGAFT#SMI CT) 
COMPUTE A PORTION UF THE SPATTAL DERIVATIVE OF THE PRESSURE, 
DPPNZH & NPPOZH « SGAFT 
COMPUTE SPATIAL DERTVATTVE "1F PRESSURE, 
DPPNZH = € GAF( ZH Y/GAMT ) #(SARPPHETZ © TT®OPPOZH/TIF) 
COMPUTE FRACTIONAL FRRUR FROM HYDROSTATIC EQUILTBRIIIM, 
FERSEQ s @(NPPENTH/(RHMAGAF( ZH )) * 1,0) 
COMPUTE PRESSURE (DYNES/CMR82), 
PP s PRPeTTeSk 
CIMPUTE DENSITY SCALE HFIGHT (KM), 
DRUNZIN = DRIDIN®TAUM(RE1 20077) eRETHI 
HARD = RHN/ARONIN 
RFETIRN 
eNov 


98 


ORIVER,751 
ORIVER,752 
ORIVER,753 
ORIVER, 754 
DRIVER,755 
ORIVER, 756 
DRIVER, 757 
DRIVER, 758A 
ORIVER,759 
URIVER, 760 
DRIVER, 761 
ORIVER, 762 
ORIVER, 763 
ORIVER, 764 
DRIVER, 765 
URIVER, 766 
ORIVER, 767 
ORIVER,76R 
DRIVER, 769 
URIVER,770 
ORIVER,771 
DRIVER,772 
ORIVER,773 
ORIVER,774 
ORIVER,775 
DRIVER,7I4 


BLKCHM 


aLUCK DATA RULKI HM 


COMMUNSOHE MRS 
FoR RFACTTUN I, 


AR(CO6), BR( 66), 


CK (66) 


RATE = ARCIV(T/30N) ##ARC IRE KPCOCR(T)/T) 


DATA AR/5S,0F e110, deUF el OeleOFol2n4 VE 1,4 ,VEal0,1,0F wie, 4 UE mle, 
Le2FasGee OF atin eMFoesard UFereslycEol9,1, OF m7, 3, 0E MTs 
PelFeNy, ce e4F ot selecFol7s1 VEIT, 5, 3E 12,6, 9F H12,2,7F ols, 
7, 0Fel te), 3b 09,5, 06 10,2, SEMIN, 1, 4E 29,2, TEM, 3, OE M19, 
14d, 3, SEI 2 LEMU UG, OF O10, U4 UEe1{,2,5F oi, 1,bF oR, 
LeOF elt SOF atin te OFoljs6,cF msi ,le5Fo11,8,5F 013, 4, OF M10, 
Pe dF eI, 4 NE O11,1, 5EV,4, NE wl12 7, PEOL1, AE e7,5,bF #26, 
be OF eA, 1 OF O71 CE woe 1 SF O41 NE wlO, 1, E15, 7, 0EeI10, 
2e0F HO 2ecF at Ores 0Fat7,5,0F9729,7,5Fa9,5,9E 07,9, 0F mle, 


be dF el 2,1 OF eles beech ahs 
DATA RRS 09% 
“5.4%, 
0.5, 
0.0 
1,98, 
on n2, 
1.5, 
71,9, 
2075. 


O00 
+2a8 


0,9 
C0 RSs 
40.50 
0,0. 
1,5, w1.71, 
0,150 0,06 
O,Ue. 1,45, 
1,05 0,08 
N6%, 3,00 
Q,%5 015s 


9.0% 


DATA CRS 0206 
"9000.5 
00% 


Oe%De 


0a%e 
02% 
31900,, 
90% 


0,06 

0.06 
U30Mes 
P600Nss 
T194., 0005 623%,, 
0,0. 277M 1e.s 0,9, 
237H9,, 0,9, ATTA,, 
900, S42004, 2u60n,, 
V—%, O,0e 460%, 
2370%,, 3900,/ 


orwrryowv- 


OeOs 


#0,35, 


9% 
2A ,99, 
0,0, 
#30, 
1,6, 


Zhe B16 
0,06 
23590,, 


0,0. 
5,0, 
lee 
e1,0, 
0,9, 
40,4, 
70,5- 
0,0, 
o1,%, 


0,0, 
oles 
40,5, 
40,5, 
70,97, 
0,%, 
71.5, 
OUs 
0,%, 


0,0. 
e075, 
1,0, 
0,0 » 
°0,¢e3%, 
0,0, 
en u, 
S25 58 
°,1S, 


20060,, 
Q,0, 
0.0, 515005 

379004, 8906, 

C5423S ee Oe Me 

0.0% C1079 
6113,, 0,0, 0,0. 
GAG, Onde, Ne tis 
S1429,, 0,0. 0. 


0.0, 
0,0. 


0.0, 
0,V0, 
00, 
5590,., 
PBUSR,. 
0,0. 


690,, 


20,75, 
71.0 


0,06 


0,0, 


0.0. 
0,0> 
0,0, 


0,0, 


0+ 


ORIVER,777 
DRIVERP,77A 
ORIvErR,779 
DRIVER, 789 
ORIVER, 784 
DRIVER, 78? 
ORI VER, 783 
DRIvER, 784 
DRIVER, 785 
VRIVER,786 
VRIVER, 787 
DRIVER, 76R 
DRIVER,789 
ORI VER, 790 
URI VER,791 
ORTVER,792 
ORIVER,793 
VRiVER, 794 
ORT VER, 795 
VP (VER, 796 
VRIVER,797 
DRIVER, 79R 
DRIVER, 799 
DRIVER, RON 
ORIVER,RO1 
VRIVER,802 
DRIVER, AD 
DRIVER, ROU 
DRIVER,&805 
VRIVER,AO6 
ORIVER,ROT 
DRIVER,ROR 
ORI VER, 809 
DRIVER,AIN 
VDRIVER,AI1 
DRIVER,AL2 
VRIVER,ALS 
URIVER, Alu 
DRIVER,ALS 
ORIVER,AL6 


[SP 


FITTER 


SUBROUTINE FITTERONPTS,%,Y, Ne IKINO,LSIGN,Z) 


ce) 
an 


n 
C 


NPTS 
xCT) 


Yy(1) 


NQ) 
IKING 


TSIGN 


7 
ion] 


705) 


oO 
oO 


DIMPNSTON 
NOL s NMe4 


INPI'T PARAMFTFRS 


NUMBER OF DATA POUYNTS 
VALUES OF THE INDEPENDENT VARTARLE, E,Gye 
ALTITUDE, KM 


SUBROUTINE FITIFR USES THE METHOD OF LEAST SQUARES TO COMPUTE 
THE COEFFICTENTS, Z2(J)eJ21,NO IN & PULYNOMIAL UF DEGREE NO 
WEPRESENTING THE NEPENDFNT VARTABLE YCI) COR, OPTYONALLY, 
NATURAL LOGARTTHM) SPECTFTEN (AND GIVEN EQUAL WEIGHTS) AT 
NPTS VALUFS UF THE INDEPENDENT VARTABLE XCI), 


17S 


VALUES OF THE DEPENDENT VARTARLF, E.G, SPECIFS 


CONCENTRATION, 1,/CM*es 
DFGREF OF POLYNOMTAL TO BE FITTED 
INDEX FOR KIND OF EQUATION TO BE FITITFD 
1 IF FQUATION JTS 

LNCY) £ AN @ AqaxX © APRKHHD & 5,4 + ANS KORN 
2 TF FQUATION [8 

yY 2 AN + ALwxX © APRKMH2 © 6,4 + ANE REN 

INDEX FOR SIGN OF EXPONENTS 
1 FUR NEGATIVE EXPONENTS 
2 FOR POSITIVE EXPONENTS 


OUTPUT PARAMETERS 


THE LFAST*SQUARFS FIT COEFFICIENTS, 
ZC1) CuURRFSPONDS TU An, Z(2) TO Al, ETC, 


AC20,e1)6 KC109), YC10%)» 26209 


12 NM2@ = NMe2 

{4 On 9 fat,NOt 

15 09 9 Jzt,NUP 

4 AfI,J) = 0, 

25 9 CONTINUE 

$4 OM 20 TetsNPTs 

3h Rs yy) 

BA AC(1,1)9 = A(1,1) + 140 
ui Gf TO (10¢12). TKIND 

5? 10 R = ALUG(R) 

69 12 § = XT) 

6? GO Thi (94916)5 TSTGN 

70 14 § = 1,0/S 

1h 16 0 ® t,v 

77 AC1L,sNN2) 2 AC(H,NiI2) + k 
10? NN 18 Je2,NOq 

114 Y 2 Wes 

15 Af(y,J) 8 A(1,/) ¢ @ 

114 18 AfJ,NNI2) & ACJ,NII2) © DER 
121 DN 20 Ke2,NN]4 

134 Q@ 2 W#8 

135 PO ACAANOL) & ACK,NO1) # @ 
1u5 Oo 30 [seen 

1ur DA 30 Je1,NN 

155 ACT,J5) & ACTal, Jeo!) 

154 BO CONTINUE 

16s CALL SL vFCA,Z,NI1) 

{70 WE TILA 

tnt END 


DRIVER AIT 
URIVER,AIR 
URIVER,AI9 
ORI VER, R820 
VRIVER,AR2) 
ORIVER,R22 
URIVER, R23 
ORIVER,ROU 
ORIVER,R 25 
ORIVER, R26 
OR] VER, B27 
ORIVER,B2R 
ORIVFR,& 29 
ORIVER,A30 
VRIVER,ASI 
ORIVER, 832 
ORIVER, R34 
ORI VER, R34 
DRIVER ASS 
ORIVER,R 36 
DRIVER, 837 
ORIVER,A3R 
DRIVER, R39 
URIVER,RUN 
URIVER, Rul 
URIVER, Ru 
URIVER,RuU3 
URIVER, Aud 
ORT VER, AUS 
DRIVER, Rub 
DRIVER, RUT 
DRIVER, AUR 
URIVER, AUD 
DRIVER, ASO 
DRIVER,ASI 
ORI VER AS? 
DRIVER, ASS 
ORIVER,RS4 
ORIVER,RASS 
ORIVER,RS6 
DRIVER, AST 
DRIVER, ASA 
DRIVER,ASO 
ORIVER, B60 
URIVER, Fol 
DRIVER, Ao2 
ORIVER, ROS 
ORTVER, ROG 
ORIVER, ROS 
URIVER R66 
UR] VER RO? 
DRIVER, ROR 
ORT VER, 869 
VRIVER,ATO 
ORIVER,A71 
URIVER,&72 
DRIVER,R74 
ORIVER,A74 
DRIVER,RIS 


IONOSU 


SUBRUNUTINE TONDSUCII,74) 


[3 SUBROUTINE TONGSY PROVIDES THE PRNPERTIFS UF THE AMRIFNT 
[2 IONOSPHERE REQUIRED BY ALL THE CHEMISTRY MUNULES, 


t HOWF VER, THIS THIRD VERSION OF TUNOSU IS LIMITEN IN THAT THE 
C PROFILES OF JTONOSPHERIC PRIIPERTTES ARF REPRESENTATIVE BUT 


(¥ NOT 


oO 
_ 
a 


THE 
(4) 


(2) 
(3) 
(4) 


> 
Oo 


THE 
(4) 


iam 
nr 


> 
c 


THE 


(1) 


NECESSARILY THE FINAL SFLECTINNS, 


Ex AND FeREGINN CHEMISTRY MODULE REQUIRES. 6, 

YOysCCMawd SEC) SF FFQ, THE CFRFECTIVE TOTAL INN 
PRIIDUCTION RATE THAT REPRODUCFES THE AMBTENT IONOSPHERE 
WHEN IISF() WITH THE CHEMISTRY MODEL, 

Ve(1sCmaed) = EFOP, THE POSTTIVE ATOMIC INN OFNSITY, 
Me(1/CMaed) 2 EFMOLP, THE POSTTIVE MOLECULAR TON NENSITY, 
TX(DEG KY, THE FLECTRON AND NO VIRRATTUNAL TEMPERATURE, 


DeREGTON CHEMTSTRY MUDULE REQUIRES so, 
Wey sCCrews SEF)D) & HQ, THE EFFECTIVE TUTAL TUN PRODUCTION 


DRIVER, R76 
ORIVER,AT7 
DRIVER, ATA 
DRIVER,A79 
DRIVER, ABD 
DRIVER, ABI 
DRIVER, Rb? 
DRIVER, BBS 
DRIVER, AB4 
DRIVER, ABS 
DRIVER, Abb 
URIVER, ABT 
DRIVER, RBA 
ORI VER, R89 
ORT VER, ADO 
ORI VER, AOL 
DRIVER,A92 
ORIVER, 893 
DRIVER, AGU 


QATE THAT REPRODICES THE AMBIFNY TUONUSPHERE WHEN I'SE() WITHORIVER,&9S 


THE CHEMISTRY MODFL, 


INPUT PARAMETFRS 
ARGUMENT LIST 


JJ = CALCULATION FLAG 
= ty CALCULATE INITIALIZATION PARAMETERS 
= 25 CALCULATE TONUSPHERIC PROPERTIFS 

7H @ ALTITIIDE OF INTFREST CF KM) 


ATMQUIP COMMUN 


IMURN, SNT(1), SNTC2), SNT(3), TT 


AL TORN COMMON 


ALTKM(47) 


RATE FIINCTTON 


RATE 


OUTPUT PARAMETERS 
ATMISIIP COMMON 


SNIOC 9), SNI(10), SNIC11). SNI(12) 


TONNE COMMUN 


FFE, EFUP,FEFMIIP, TX, UDEF 


VARTABLES TN TOANOUP 


EFESSNI( 9) © ELECTRON DENSTTY IN Ew AND 
FeRFGTON, 1/CM*@5 
FFEOPSSNIC{0) @ ATOMIC OXYGEN ION DENSITY IN Fe AND 
FeKEGTON, L/CM883 
FEMOLPSSNIC{1) © MOLECIILAR ION DVENSTTY IN Fe AND 
FeREGIUN, 1/CMe8S 
TxESNI(C 12) © ELECTRUN AND NO VIBRATIONAL 
TEMPERATURE, DEG K 
QNFR «© EFFECTIVE TOTAL ION PRUBUCTITON RATE, 
1/(CMaext SEC) 


RFUITREN QUANTITIES FOR THE Ew AND FeREGION CHEMISTRY ARE 


HATAINEN 48 FOLLOWS ees 

EFQ TS CUMPLITFO FRUM & MUNIFIFD FORM Kh EQ, (28276) OF 
DASA 2u97e1 (GE#7%),, 

EF = FFALPDaFRF aE FR a1, © EFALPRREPF/FRFRETAI/(1,9 © 


DRIVER, R946 
ORIVER, AO 
URIVER,R9A 
ORL VER,R99 
ORIVER,900 
ORI VER, 901 
URIVER,902 
URIVER,903 
URIVER, 904 
ORIVER,905 
VRIVER,994 
DRIVER, 907 
ORIVER,90AR 
ORIVER,909 
DRIVER, 910 
URIVER,911 
DRIVER,912 
ORI VER, 913 
DRIVER, 914 
ORIVER,915 
DRIVER, 9146 
ORITVER 917 
DRIVER, 918 
ORIVER,919 
ORIVER, 920 
VRIVER, 921 
DRIVER, 92? 
DRIVER, 923 
DRIVER ,%2u 
VRIVER, 925 
DPIVER,926 
ORIVER,927 
DRIVER, 928 
URIVER, 929 
VDRIVER, 950 
DRIVER, 9351 


CFF ALPOMEFSF + FEAL PR#( 1 PE FSF) mEFE/SEFHE TA) VRIVER, 932 


WHERE VRIVER, 933 
EFE = ELECTRON DENSTTY PROVIDED aS INPUT DATA TH ORIVER,934 
TONOSU (1/0Max3) DRIVER, 935 
FFSE = SNIC3)/0SNI 03) & 2,e(SNT(1) SSI (2))) URIVER,9$6 
SNTC1) = w2 CONCENTRATION ORIVER,937 
SNTC2) = 02 CONCENTRATION VRIVER,93R 
SNT(3) = 9 CONCENTRATION ORIVER,949 
FFALPND = DTSSOCTATIVE RECUMBTNATIUN RATE CNEFFICTENT FORK URIVER, 940 
THE REACTICIN Me © —& & PRIMUCTS (CMee3/SEC) ORIVER, 9a 
= RATF(13,TX) # RATF(14,TX¥) WhERF RATE CINOY, TEM) [S ORIVER, Su? 
THE FIINCTTON ROUTINE FOR Fo AND FeREGION RATE DRIVER, Ou 
CORFFICTENTS, (CMe 3/SFC) DRIVER, Suu 
FFRETA s RATFECIO,TX)*SNI CL) © RATECO,TT)#SNI (2) (1/SEC) DRIVER,94S5 
RATFE(10,TX) = RFACTTUN RATE COEFFICTENT FOR THE REACTTIIN DRIVER, 946 
()}¢ @ NO s Ne @¢ N DRIVER, 947 
RATFC 9,TT) = KFACTIUN RATE COEFFICTENT FOR THE REACTION DRIVER, 9UA 
ne ¢ 02 2 bet + 0 DRIVER, 949 
(2) EFOP Ts COMPUTED FROM a MIUIFTED FURM UF FUu,(2e274u) UF ORIVER,950 
DASA Pu97"{] (GE"70),, VRIVER, 954 
EFUP = FRSFREFU/S(FFRETA + EFALPRREFF ) DRIVER,952 
wHE PE DRIVER, 983 
FRALPR = EFFECTIVE THURBIDY CULLISTUNAL eRADITATIVE RECOMBINA® DRIVER, OSU 4 
TTUN RATE COEFFICIENT FUR ATUMIC TUNS DRIVER 95S 
= RATE(11,TX) # RATF(1SeTX)REFE © 1,5FS07#SORT(FFF)/ DORIVER, 95H 
T¥e*3 (CMARS/SEC) DRIVER,9S7 
(3) EFMALP TS COMPUTED FROM A MODTFTED FORM UF FU, (29275) UF ORIVER,9SA 
DASA 24u97@} (GE#70), DRIVER, 959 
EFMOLP = ((1,eEFSF)BEFM + FFRETAREFE)/S(CEFALPDRERF «© FFRETAY ORIVER, 960 
(4) TX(MEG K) ITS FOMPUTFD FROM AN TTERIM PRESCRIPTION, OR] VER,961 
cc DRIVER ,%6? 


FLECTRUN DENSTTY PROFILES FOR NOMINAL MIDLATITUNDE DavytyTMe AND DRIVER, 963 
NIGHTTIME CONDITIONS TN THE Ew AVN FeREGIONS ARF PROVIDFL AS URIVER, 964 
APPRIXIMATE FITS TU CURVES IN FIG, 1 TF Hy RISHRETH, PHYSICS URIVER, 965 
AND CHEMISTRY UF THF TUNOSPHERE, CUNTEMP, PHYSICS, VOL, 14, ORIVER,%66 


P, 229(1973) (RIe75), URIVER 967 
(qa ORIVER,96R 
FOR DAYTIME FLECTRUN DENSTTY eee ORIVER, 969 
cc VRIVER,970 
ASSII4ME PARARILIC TNCREASE IN LUG DF ELECTRON DENSTTY FROM VRIVER,O7$ 
ALUGLOCFROTN) & 5,0 47 ALTITUDE HEBNTO = 100.0 KM TN ORIVER,972 
ALOGLOCEFAMYDY = ALNGIO(C7,~SF495) AT ALTITUDE HFAMX¥D = 590, KMsDRIVER,973 
FOLL IWED AT HTGHER ALTITUDE KY EXPONENTIAL DECKF ASE WITH URIVER,974 
SCALE HEIGHT FENSTH = 20, KM, BFLOW ALTTTIDE HERUTD, ASSUME DRIVER 97S 
EXPONENTIAL DFCREASF WITH Salt MFIGHT FUNSTH = 5,0 KM, DRIVER, 976 
URIVER,977 

IF (7H{GT HF AIMYD) EFF = EF OMXNME RKP (CHF OMXN@ZH)/FONSTH) DRIVER,9O7A 
IF C7HGEFROITN g AND, 7HALE HR AMYIY ORTVER,979 
FEE = FEOMyN ef Ugee (FRE AR (HE AYKU M(H) ae) URI VER, ca0 
WHERE THE COEFFICTENT EFEA TS DETERMINED SU THAT FFF & FBOTN ORIVER, 981 
AT ALTITUDE HEROIN, ORIVER, 982 
TeEoe FFFA & ALIGIOCERUID/EF 249) /(HFAMXDCHEBIITN) a22 DRIVER, 983 
ayTH URIVER,9R4 
HF OmMyN = ALTTTIHDF OF FAMAx IN DAYTIME, KM DRIVER, 98S 
FFAmMyYy s ELECTOAUN DENSTTY AT FamMAx IN NAVTIMEF, 1/0 Naat DRIVER, ORS 
FROTO s FLECTRON DENSTTY AT HFROTD, T/C mMaed URI VER ,967 


102 


CACAO MCHA ALOAAOSSOANAOS AO OCAM A AMATO KC AAMAMAHOAODMOANMOCMeOR Al CMAs oe 


> 
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is 
a 


> 
oO 
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Lae] 


+ 
o 


o 


im} 
ao 


TECZALTHEBUTD) EFE = ERITDSE KPO (ZHeHEBITNO) ZERDSCH) 


FOR NIGHTTIME ELECTRUN DENSTIY ee, 


ASSIIME SINUSITD INCREASF IN LAG JE FLECTRON DENSITY FRUM 
ALUGIOCEBOIN) = 8,0 AT ALTITUDE MERDIN B 100, K™ TL 


FOLLIWEN AT HIGHER ALTITUDE BY EXPUNENTTAL DECREASE 
HEIGHT FONSTH = PND, KMA BELIIM HEKNIN, ASSUME EXPOUNERTIAL 
DFCREASF WITH SCALE HEIGHT FONSCH = 5,U KM, 
TE(7H,GT,HEOMEN) FFE = EFOMXN®FXP( (HF OMXN@2H) /FONSCH) 
TF C7H,GE,FROTN AND, ZH,LE,HRAMKN) 
ALUGTOCEFFE) = ALIGIOCE RNIN) ¢ O,5N#ALNGIOCEFOMXN/EBDTN) 

& (CL MOSIN(PING (A, eLHeHERUTNGHFOMXN) / 

CHE OMXN@HEBOTN) )) 

IFCZ7H,LTHERUTN) EFE = ERDTN®EXP((ZHRHEANTIN) /EDNSCH) 


EL ECTRIUN TEMPERATURE PROFILES IN THF Fe AND FeREGTUN ARE 

DATATNED, FOR (NGAN)Y PAYTIME CUNDITIUNS, RY PRESCRIRING THE 
DIFFERENCE RETWEEN THE FLFECTRON TEMPERATURE TX AND THE GAS 
TEMPER ATI2G TT AT TwO ALTITIIVES AND USING A FARABNLIC FIT 
TN THIS DTEFFRENCE, FUR NIGHTTIME CUNDITIUNS, we 


FOR DAYTIME FLECTRIIN TEMPFRATIIRE,,, 


ALTITIIDE, KM TXeTT, NEG Kk TT(CTkKAe6S, MINEL@S, BeHR) 
1en iD 


200 500 


TXT120 335 
TxTe00 933 


THESE VALIIES NF 
REPIIRTEN Fy JeVv,y FVANS 
FIR 196A AND 1967, PLANET, SPACF SCT, VUL~y 21, PP, 7639792 
C1973), (Fve73)) AND THE CTRA®I9KS MODEL ©S BeHR ATMOSPRKERE 
(fl*65), 


T¥eTT ARE CONSISTENT WITH THE VALiIES MF TX 


PECZTH ELT SU) TK = TT 
IFC 7HGE.120,) TT¥T = SARTC ZHMI2eNsA Y 
WHERE 
7HM1?P0 = 7Helen, 
A = BO, 4 SN0e #82 


THE QEQIITREN QUANTITY FOR 
AS Ful LAWS, oe 
DA TS FORTEN Tu EQUAL THE VALIIE UF FFQ AT 


THE DeREGTUN CHEMISTRY TS UYRTAINEN 


THe oATTiM OF THE 


GRIP (90ekKmM) AND 15 DETERMINED BY INPIT DATA AT LOWER 
AL TITIIDFS, 
NOTE ge, GOFF 3 PO OK GNEF £ FRO NEPENDING TIN THE 


ALTITUBE Zhy 


FIR DAYTIME yey 


TF C75 .LF AO.) 
PH = MVNAY(Y) & ONY INT ae (2H 72077421307) 
ILIO7F SF ANVWNAVCVSI/NGNAY (7) 
ZHMI07 = TRHOALTKM(7) 
Z13“07 = 


AL TKM( 13) e4LTKM(7) 


DRIVER, IBA 
URIVER,9K9 
URI VER, 990 
VRIVER,991 
URI VER, 992 
VRIVER, 993 
ALUIGLOCFRAMX¥NY © ALIGIOC, OF 405) AY ALTITUDE HRAOMXN = 360, KM,ORIVER, 994 


aTTeK SCALFURI VER, 995 


URIVER,996 
DORIVER,997 
URIVER,99R 
ORI VER,999 
URIVER, 1000 
ORIVER, 1094 
DRIVER, 1002 
ORIVER,1003 
ORIVER,1UNG 
URIVER, 1005 
DRIVER, 1096 
ORIveER, 1007 
ORIVER,10%8 


ASSUME TxeTTURIVER, 1009 


URITVER, 1010 
DRIVER,1011 
DRIVER, 1012 
DRIVER,1013 
URIVER, 1034 
DRIVER,1015 
ORIVER,1016 
ORLVER,1017 


(MILI STUNE HILL THOMSUN SCATTER RESULTSDRIVER, 101A 


DRIVER,1019 
URIVER, 1020 
ORIVER, 1021 
ORIVER, 1022 
VRIVER,1023 
DRIVER 1024 
URIVER,1025 
URIVER, 1026 
URIVER,1027 
URIVER,1026 
URIVER,1029 
VRIVER,1030 
ORIVER 1031 
DRIVER, 103e 
VDRIVER,1035 
DRIVER, 1034 
DRIVER,1035 
ORIVER, 10386 
ORIVER, 1037 
ORIVER, 10368 
VRIVER,1039 
DRIVER,1040 
DRIVER, 1044 
ORIVER INU 
URIVER, 1043 


C TF COUg el TeZh pAND, ZHelLT, 90) DRIVER, 1044 

C OQ = NUNAY(13) * ANG See ZHM71372459M1 5) VRIVER,1UUS q 

ce WO1913 = FFAZIO/DOLAY(4 $) URIVER, 1046 4 
€ ZHM71%3 = 7HeALTKM( 43) DRIVER, 1047 

C ZIIMLS = AL TKM(19) AL TKN (135) ORIVER,104R 4 
crc DRIVER,1049 
[% FOR NIGHTTIME see URI VFR.1050 
C URTVER,1051 
5 IFC ZHLF 60.) DRIVER,1082 
C OF = OWNITO7) & GNI30O7TH#C7HMZN7/Z1 507) DRIVER, 1053 
€ QAN{3O7F = NANITCH3ZI/NGNIT(7) ORIVER, 1054 
c UPIVER,1055 
€ LF CAO. Ve2e eANO, CHET 906) DRIVER,1056 
C ON =F OGNIT(LS) *& ONIGI3*#*#OZHM71 3771913) ORIVER, 1057 
c QN191% Ss FFOZI9O/0OONIT(43) DRIVER,1056 
CCC ORI VER, 1089 
DTMENSINN AYMAY(C18),90NIT(18) URIVER, 1060 
COMMONSATMIUPS HL SBARSTDOURN, PP,RRO,TT,SNI (C16) sHRHO,FEHSF OQ URIVER, 1061 
COMMONSTONUIIPS EFE,EFUP,FFMULP,1¥, O0FF ORIVER, 1042 
COMMUN ZALTUONS NALTON,ALTRM (U7), NDAY C27) ,UNITE(IB),S1Z2N,CO2(25) DRIVER, 1063 
COMMIIN/ZTHOHEX/Z ZHELAG DRIVER, 1064 
cee DRIVER, 1065 
DATA ERITO,HERDTD EF AM¥D sHFOMKN,FENSCH,ENDSCH 4 Le0F 4O5e1 2 0F HNC r DRIVER, 1066 
* TDF HOS S.0F HNP, 2 e0F 4020560 / ORIVER,1067 
DNATA ERUTN, HEROIN, FF SMXN, HE OMXN, FONSCH,EDNSCH 4 LOE MNSe Le VF eNerDRIVER, 1068 

tO UF HI 3, OF 492A IE HAS, O S URIVER,1069 3 
DATA THXTV20,TXTPOOGTXTAOD 4 0,0,5,0E402,1,86 403 / ORIVER,1070 
DATA PT 4 3,141592658590 / DRIVER, 1071 
C INTERIM VALUES 06/10/75 DRIVER, 1072 
DATA COADAYOT I, TSLpIBISAHI go dy de D0 ee MMO, S*0g/ DRIVER, 1073 
C INTERIM VALIIES 06/10/75 ORIVER,I074 
PATA (OANTTOTI,TEle,1BV/A#0, 6b, sek ge I Nbr SHS ORIVER,1075 
crc ORIVER, 1076 
GO TO (109,200), JJ ORIVER,1077 
C INITIALIZATION, CALLED FROM SIiBRUIITINE ATMUSU LIIRTNG TTS DRIVER 1078 
C INITIALTZATTIIN, URIVER,1079 
11 190 CONTINUE DRIVER,10R0 
11 PTD? s Psa, ORIVER 101 
13 HPPRIDP = 0,508 (HE AMXNOHERIIN) DRIVER, 10R2 
{4 HAMAR? = 0,50 # (HF AMKN@HFE OTN) UVRIVER,1083 
en AlLG2D1 & O,5N*ALUGLOCER EMXN/EBOTN) DRIVER, S0AR4 
au FFEA & ALOGIOCERNTD/SEFOMKD) J CHF OMXDeHERNTD) #82 URIVER, 1085 
34 A = 89, / (509,#590¢) ORIVER, 1086 
G INITIAL TZ4TTIN FIP DwREGION eee DRIVER, 10R7 
C COMPUTE ELECT TRUN TEMPERATIIRE AT 9QekM ALTITUDE ORIVER,10A8 
36 TY 2 TF ORIVER, 1089 
37 TF(TONRN LTO) GY TO 150 ORIVER,1090 
‘ COMPUTE DAYTIME FLECTRON DENSITY AT 90 KM URIVER, 109} 
u2 EFe = ES JTD # EXPEC (ON, @HEBITIN) JENDSCH) ORIVER,NU9e 
51 GO Ti) LAO ORIVER, 1093 
c COMPUTE NIGHTTIME ELECTRON NENSTTY AT 9%eKM AL TITIIDE URIVER, 1094 
53 180 BFE 3 ERITN # EXP((90, © HEBUTN) /SFUNSCH) ORIVER,1095 
6% 1A EFALPN = RATE (13,7) @ RATE €44,7XK) ORIVER, 1096 
1? EFAL PR = RATE CLI, TXIFRATE (CIA, TKIRFFF SI SE@OTRSURT(EFE) STK ORS ODRIVER,1097 
110 EFSF BE SNT(3YS(SNT(3) © 2,e(SNIC1) © SNT(2))) VRIVER, 1098 
4446 EFRETA = RATE CYAN, THIYRSNTES) & RATECI,TTIBSNI( 2) DPIVER, 1099 


125 


140 
{44 
1458 
1u7 
150 
1S? 
154 
154 
1546 
164 
fol 
Yo! 


165 


FFIV710 s FRALPebPE REFER (1.0 ¢ FRALPROeEFE/EFRFTA) / 
(1,0 © (FRFALPD®EFSFE + FFALPRe(}  mEFSF))RFFE/SFF RETA) 

IFCIDORN, LT.0) Gu TO 990 

Q%1913 2 FFAZIQG/NIDAY()3) 

INLSO7F = OVNAYCHSISNYNAYC7) 

6 Ty 195 

QN1{9O13 2 FFNZIQ/NANTT (13) 

QNL3O7 & MONITCYSI/0QNTT(7) 

CONTINUE 

Z19M13 B= ALTKM(19V PAL TKM( 15) 

Z13MO7 = ALTKM(13IVPALTKM(7) 

WF TIIWN 

CONTINUE 

IFC Z7HNESZHFLAG 2 CALL ATMI)SUCE, ZH) 


AN FRRONENMIS CONDITION WILL OCCUR TF TONQSU IS CAaLLFD WITH 
Jus2 AND A GIVEN VALUF NF 24 TF ATMINSU HAS NOT REFN CALLEN 
FIRST wITH TJ22 AND FOR THE SAME VAlLUE DF 24, 

THE VARTARLF 7HFLAG IS JISED TN DETECT THIS CUNDITIUN AND 
TO MAKE THE REQUIRED Catt FO ATMUSU, 

ZHFLAG TS INITIALTZED TO AN ARBITRARY NEGATIVE VALUE YN 
THE INTYJALTZATTON CALL TO ATMUSU, 


TFC7H,GF 290.) GO TU ans 

SFT ELECTRIIN TEMPFRATUKFE FUR 7H,LT, 90, 
TX = TT 

ZFRO FRE, €FORP, AND £FMNLP FOR ZHALT.9U, 
EFE = €FUP = FFMJLP = 0,0 

PROCEED wItTH MO CALCULATION FOR ZH,LT,90, 
IFC TONRN,LT,0) GU TN 350 

COMPUTE DAYTIME DA 
TFC7H LF 260.) GM TU 325 

COMPUTE DAYTIME NA FOR Hoel T, 7H LT,90, 
ZHM714% # 7HeAl TKMO13) 
ON = AYNAY(13) #& ADI913#eOZ4M743/719M13) 
69 Ti) $85 
CONT [NIE 

COMPUTE DAYTIME OA FOR ZHyLF,60, 
Z4M707 = 7HeAal TKM(7) 
NA s NUNAY(7) w UNLSOFR#(ZHM7Z07/21 3M07) 
GN Ti) 385 
CONTINUE 

COMPUTE NIGHTTIME DO 
FREZW LE eO0'c)} GN TO S7S 

COMPUTE NIGHTTIME 10 FUR 60 ¢LTeZ7H LT, 90. 
ZHM714 & 7HeALTKM(13) 
NA B MGNITONS) & ONI9D1 See 02HM7137719M4 3) 
GM Ti) aS 
CONTINUE 

COMPUTE NIGHTTIME 09 FUR 7H,LFe60, 
Z4M7U7 & 7HeALTKM(7) 
OA = NINITCT) # UNL3SOTER(ZHMZN7/21 3M07) 
OOEF ss Ay 
SNI¢ 9) 


ia) 
2 
_ 
ae 
> 
nun 
— 


105 


URIVER,1100 
URIVER, 1101 
URIVER, 1102 
DRIVER 1103 
ORIVER, {104 
DRIVER, 11095 
DRIVER, 1106 
ORIVER, 1197 
DRIVFR,11068 
ORIVER,1109 
URIVER, 1110 
DRIVER, 1114 
DRIVER, (112 
URIVER, 1113 
VRIVER, 1414 
DRIVER, 1115 
DRIVER, 1116 
URIVER, 1157 
DRIVER,1118 
DRIVER, 1119 
DRIVER ,1120 
ORIVER, 1121 
URIVER,112e2 
DRIVER, 1123 
VRIVER, 1124 
DRIVER,1125 
ORIVER, 1126 
ORIVER,1127 
ORIVER,1528 
ORIVER 1129 
ORT VER, 1130 
ORIVER, A144 
VRIVER, 1132 
DRIVER, 1133 
URIVER, 1134 
VDRIVER, 1145 
URI VER, 1136 
ORIVER, 1137 
VPIVER, 1138 
ORIVER, 1139 
URIVER, 1140 
DRIVER, 11435 
DRIVER, 1142 
NRIVER, 1143 
VRIVER, 1144 
DRIVER,1145 
DRIVER, 1146 
DRIVER,1147 
VRIVER,11408 
DRIVER, 1149 
ORIVER,1150 
VRPIVER, 1151 
ORIVER 1152 
ORIVER,1153 
URIVER, 1154 
ORIVER,1155 


259 
?5? 


253 


255 
as? 
Pas 
Pot 
272 
to? 
toe 


sat 
B14 
B14 
16 
320 
325 


427 
332 
tan 
sue 
yas 
Yo! 
363 


372 


374 
40% 
net 
4es 
nut 


4g? 
45? 


uss 
46s 
Una 
465 
ue? 
a7! 
ur! 


SNIC12) = TX 
RETIIRN 

IFC INANE T,9 9 GO TN 259 
COMPUTE DAYTIME ELEC TRUN 
be AND FeRFEGIONS, 


DENSTTY AND TEMPERATURE OF 


ELECTRON DENSTTY 
IFC ZMHeHERNTD ) 219,212,212 
BFE = EB)TI w EXPO (ZHeHFHOTN) SENDSCH) 
GN TO 220 
LF (  ZHweHF AMX!) ) 
EF = EFQMxN * 
GN TO 220 
EFE = EFOMXN 
ELECTRIIN 
IF( ZHel2an, ) 
Tx = Ry 
GN TO 2A0 
ZHM120 = 7Hel2n, 
T¥ = TF € SGRT¢ 
G9 TO 2AO 


214,714,216 
Le k* CF RE RE (HE OMX De LH) a) 


* EXP (CHE OMXD=7HI/FANSCH) 
TEMPEWATIIRE 
22202244724 


ZHM4tean7a ) 


COMPUTE NTGHTTIME ELECTRIUN DENSITY AND TEMPFRATIIRE OF 
Ee AND FeREGINONS, 


ELECTRON 
IF( ZHeHERUTN 1 269,262,262 
tFe = ER ITN @ EYP CO (7HeHFROTN) ZENNSCH) 
6 TU 270 
TF ( ZHeHFAMXN ) 
EFE = ERVIN & 
69 TD 270 
EFRE s EFOMKN *# FP (CHF AMNXN@ ZH) /FANSCH) 
ELECTRIIN TEMPFRATIIRE 
t¥ | TF 


DENSTTY 


264,764,266 
19,88 ALGODL# (1 OeSIN(PINS a ( ZHeHAPRD2) /HPMADC))) 


COMPUTE EFW, FFOP, AND FFMILP 
EFG 
EFA, PN = RATE(43,7TX) © PATECI4,TX) 
EFAL PR = QATEYI1, TX) © RATECLA,TXIBEFF # LeSEPUTRSURT(EFEV S/T XaH§ 


EFSF = SNTCB)S(SNT(3) © Set (SNICLY*SNI(2))) 


FFRETA = PATE(YO,TX) # SNI(1) @ RATE(9,TT) * SNI(A) 

EFU = FFALPDSFRE RERF a (1, OFEFALPROEEFESFERETADS 

* CL, O+TEFALPNSERSFE + EFALPRO( 1, oF FSF) ®EFE/EFRFTA)D 

NER 3 FFA 
FEI\e 

BENS s FRESFREFOS(FRHETA © EF ALPRYFFE) 
EEMOLP 

FEMOLE = (CEL PEESF)*EFQ + EFRETASFFE)/S(FFALPDSEFE + EF BETA) 

SNIC9) 3 FRE 

SNIC10) = EFIP 

SNICI1) = EF MYO 

SYI(12) 2 T%¥ 

WE TIIWA 

ENT 


VFIVER, 1156 
VRIVER,1157 
URIVER, 1158 
ORIVER,1159 
DRIVER,1160 
UDRIVER, 1161 
VRIVER, 11462 
DRIVER 1163 
ORIVFR,1164 
URPIVER,11465 
URIVER,1166 
ORIVE® 11467 
DRIVER, 11468 
DRIVER 1169 
URIVER, 1170 
ORIVER, 1171 
DRIVER, 1172 
DRIVER,1175 
VRIVER, 1374 
VRIVER, 1175 
URIVER, 1176 
ORIVER 1177 
DRIVER,I178 
URIVER 1179 
DRIVER,14A0 
DRIVER, 1184 
DRIVER, 14R2 
DRIVER, 11A3 
ORIVER, 1184 
OPIVER,S1AS 
DRIVER,13R6 
UDRIVER,11AR7 
ORIVER,1188 
DRIVER, 1489 
DRT VER,1190 
DRIVER, 1191 
ORT VER,1192 
DRIVER 1193 
VRIVER, 1194 
URIVER, 1195S 
ORIVER,119%6 
ORIVER,1197 
ORIVER, 1198 
ORIVER, 1199 
DRIVER, 1290 
DRIVER, 1201 
DPIVER, 1202 
ORIVER 1203 
URI VER, 1204 
ORIVER, 1205S 
DRIVER, 1206 
ORIVER, 1297 
URIVER, 1208 
VPTVER. 1209 
UVRIVER, 1210 
ORIVER, 1211 
ORIVER, 1210 
DRIVER, 1213 


1? 
20 


20 


e3 
25 
59 
3? 
3h 
37 


CCC 


ca 
an 


oO 
c 


AGA AAOAArMAa NOOO RrAOnmm 
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oO 
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SUMRUUTINE TULTANCTYRG, ITMONG,IDAYSs VRE J, VF US,DAYI) 
SUBROUTINE JULIAN CONVERTS A GREGORTAN CALENDAR DATF TH 
JUL TAN DAY NUMBER DaydS FUR SURROUTINE SULIRB, 
SUBROUTINE JULLAN TS VALIP FUR YEARS 41994 19 1999 INCLUSIVE, 
INPI'T PARAMETERS 
TYRS @ NUMRER OF THE 
RFCOIMFS 74), 
NIIMRER OF THE MONTR 
Day DJF TRE MINTH 


YFAR IN THE 1900 S$ CEybye 1974 
IMUNS « 


INAYS « 


CF .Gee FEBRUARY KFCOMFS 2), 


UHTPUT PARAMETERS 
VRE J © JIL TAN PAY NUMRER (A HALF INTEGFE We) 
UN JANUARY 1 OF THE YEAR (1b 
VEQJ @© JUL TAN DATE FOR VERNAL FULINUK, 
NAYS «© JULTAN PaY NYyMBFR (CA Halk INTEGER) 
(IN THE DAY OF INTEREST, 


4Y 0 RRS UT 
INTEREST, 


AT O HPS UT 


DFFINITTON OF HATA 
DAYM(T) © THE CHMULATTVE NUMMER QF DAYS FRUM THE REGINNING 
OF THE VE Aw I) THE FN OF THE (Teej) TH MONTH, IN 
A NONwL FAP YEAR, 


OTMENSTON Mayne y?) 

DATA (MAYACTY,T3L 6129 FS Oop Stor SF er Mer I OOe ellen Ibleseieer 
* O43 ye PT Bee SONG, S54, / 

UVAYS = TODAYS 

VPS & [YRS 


THE FIRST TERM FOR OAaYS IS THE JULTAN DAY NUMBER AT 0 48S UT 
19090 JANUARY 1, THE THIRD TERM FOR Day t TS THE NIIMAER (IF 
EXTRA (LEAP wYFAR) LYS SINCE 4900 TM THE START OF THE YEAR 
OF TNYEREST, 


MAYS = P4I5SN20,5 * 8HS,*VRS © AINT( (v¥RS91,)/4, ) 

YRE XY = Nays 
VERNAL FullT%ax OCCURS wTTHIN ABUT 7 SECONDS IF TIME AT 
ON HAIRS ON 21 MARCH 1974. AT WHITH TIMF THE JULIAN Day 


NIIMRER 15 24UP}27,5 , FUR NEARBY YEARS THE JULTAN DATE FOR 
VERNAL FOUINOX WITL RE GIVEN AY VEU), 
VEGI = 24U2127,5 + 365,25*(YRSe74,) 


LFAP TS AN INDEX THAT ENUALS 9 
ENVALS 1, 2, MR § 


FOR A LEAP YEAR AND (THERWISE 


LFAP = MON(TYRS,4) 


IFC IMUNS.LT.3 9 GI TO 1 

IFC 1 EAP,EQN,O «6» 6DAYI s&s DAYJO1,9 
OAYJ @ PAY © DAYMCIMONS) © (MAYSe1,0) 
WETIIWN 

END 


ORIVER, 1214 
DIVER, 1215 
ORIVER,. 1216 
D@IVER,1217 
ORIVER,1218 
URIVER,1219 
URIVER, 1220 
URIVER, 1221 
URIVER, 1222 
DRIVER, $223 
URIVER, 1224 
ORIVER (1225 
VRIVER, 1226 
URI VER, 1227 
ORIVER, 1228 
ORIVER, 1229 
ORI VER, 12430 
DRIVER, 1231 
URIVER, 12%2 
DRIVER, 1233 
DRIVER, L236 
DRIVER,1235 
ORIVER 1246 
DRIVER,S237 
DRIVER, 123A 
ORIVER, 1239 
DRIVER, 1240 
URI] VER, 1244 
VRIVER,124e2 
VPIVER, 1243 
ORIVER, 1244 
DRIVER,1245 
DRIVER, 12U6 
DRIVER,1247 
DRIVER, 1248 
URIVER,1249 
URI VER, 1250 
ORIVER, 1251 
ORIVER, 1252 
DRIVER, 12535 
ORIVER, 1254 
ORIVER, 1255 
ORIVER, 1256 
ORIVER,1257 
ORIVER,1258 
ORIVER,1259 
DRIVER, 1260 
URIVER, 124} 
ORIVER, S262 
VRIVER 123 
DRIVER, 1e6u 
VRIVER, 1265 


RATE 


15 
16 


AAIAMHAONONAAAA AA OALCAOOAHAAOAHAAHAANNAANAOAHAHOAAOLCNHOHOAHS OAK OAC LN 


co 


FUNCTION RATE CINOK, TEM) 


MFACTTON RATE COEFFICTENT FOR A SPECIFIFD REACTION 


INPUTS 


INDY s REACTIO®R INDEX 
TEM & TEMPFRATIIVE (DFG K) 


(SEE HELUW) 


VRIVER, 12eho 
ORIVER, 1267 
DRIVER, 1268 
DRIVER, S269 
URIVER,1270 
ORIVER, 1271 
DRIVER, 1272 


VIRRATTONAL/ZELFCTRON TEMPERATURE FIK KEACTIUNS 7,4,10014,2SURIVER, 1273 


AND PY 
UUTPUT 


RATE = REACTICin RATE CUFFFICIENT 


RFACTTONS INCLUDED 


1 Ne ¢ J2 s Ne © OD 
BNe © OP fF Ne © ND 
S N@# @ NO FNS ON 
7 N@ @ F&F SN @ HNU 
9 Me + 2 s VWa4e 6 
1! Oe ¢ &€ s 0 + ANU 


13 NUe + EF 

15 “(4aS) ¢ 

17 N ¢O s NO @ HNU 
19 “(usd ¢ 

ei (4S) + 

23 N(4S) ¢ —E s N(2D) 
es. @¢ Ne s NT 6 N 
27 NQ0e4M2 = E*s}P4N2? 
29 Nek = UPeeh 

31 NdweM B NQweliQemM 
3B NbeeN? = Mess 

35 NONDQ#4*N2 3s NUetND 
$7 N2Q0«e4O2 = Uses 

39 N2Q4NN2e FS (Zeon) 


«#303 


G1! Ni2e#NN2 = JONNeONN 


U3 Ni)Pwe 2 = NDesNi? 
US NOPetHN = NICE 


U7 NOBZINN S&S NiIPeINNP 
U9 Xeeveoe4 = FRODIUCTS 


$1 1102eeE— = Wien? 
53 WOPesE = NED 
55 Nek & NethHNy 
S7 NeeNl fF NIPOF 
59 Net? = NEE 


61 ALO Axe & ALND OEM 


63 NeN2 2 NOS? 
65 NeetN) = NO Deel 


COMMON/SC HF MA, SYfbb), 


TMxXSTEM 


TMYSTFM/ $00, 
FY1S1.,0 

IF (RARCINOX),NE,O,9) 
Ex221,0 

TE (CROINDW) FNAL DO) 


GA TN 


BR ho), 


Ne + U2 2 024 ¢N 
N+ @¢ NQ = NQ@ # N 
N@# ¢ 0 3 Ut © N 

Ne @ E ¢ E ENE 
Ne * NP zs NU ¢ N 
OF 4 € * € 2 0 4 € 
NU #@ F = N(2D) #9 
N(20) ¢ 0 = NOt ¢ F 
N @ N & N@ @ HNU 
N(?0) #¢ O02 s NN ¢ 
N(2U) @¢ NU = NP #¢ 0 


N(20) ¢ &€ & N(4US) 
FeN2e02 = NQ84¢2 
Neer) = OFF 
NeerI24M = ieee 
Net = NSeeNe 
NGerNi) = VONDO#4+0? 
N$eeN = N2e4Ne 
NSeeNG = NN Qe4N2 
NONQetNG & NO2e4NO? 
NesetnNJ2 = NOPe4O2 
NUPOF & NUC ORNL 
NOP SAN = NIG eeNU 
¥t*Y¥e = PRODUCTS 
ALS+FeM = AL OM 
WUeer & 14D 

(+E oéM = YOM 

Net s Yer 

Neen = NOE 
ALUCOEOM = ALUPeOM 
ALO*#ee SsAL40 

NU¢U? s&s ONO? 
NU2wsh = U2esV 


+ 


CK (bb) 


ErtsTMyee (BRC TNDXK)) 


DRIVER, 1274 
DRIVER, 1275 
DRIVER, 1276 
URIVER, 1277 
ORIVER,1278 
ORIVER 1279 
ORIVER,12R0 
DRIVER, 1284 
CRIVER, 1282 
ORIVER, 1283 
DRIVER 124 
DRIVER, 12RS 
URIVER, 1286 
ORIVER,12R7 
ORIVER, 128 
DRIVER, 4289 
ORIVER 1294 
DRIVER 1291 
VRIVER, 1292 
DRIVER, 1293 
ORIVER, 1294 
DRIVER.I295 
VRIVER,12%6 
DRIVER, 1297 
ORIVER, 1298 
ORIVER,1299 
DRIVER,1300 
ORIVER, 1304 
DRIVER, 1592 
DRIVER, 1395 
DRIVER, 1304 
ORI VER, 1359S 
ORTVER 1306 
ORI VER, 1507 
DRIVER, 13598 
DRIVER, 1399 
DRIVER, 1310 
DRIVER, 1391 
ORIVER, 1312 
DRIVER, 1313 
DRIVER, I314G 
ODPIVER, 1315 
ORIVER, 1316 
ORI VER, 1317 
DRIVER, 1314 
DRIVER, 1319 
ORI VER, 1520 
DRIVER, 13?) 


RATE 


25 
%0 


190 


EXPNSeCRCINIKI/TMX 
IFCEXPU,GT,=79,9) Gil TU 2S 
RATE = 0,0 

GN TO 100 

EX2=SEXP(EXPN) 
RATFSARCINDY) REX {RE XO 


IFC TNOX,EQ,10) RATE = AMAXI(],5F oie, 


RETURN 
END 


109 


AMINE (C1, E010, 


RATE)) 


VRIVER,1522 
ORIVER, 1523 
ORI VER, {324 
ORIVER, 1525 
ORIVER, 1526 
ORI VER, 1527 
DRIVER,1328 
DRIVER, 1529 
ORI vER, 1330 


SOLCYC 


> 
ao 


Cig) MANANANSOAAANCANANAND 
a a in) (a) 
a a C1 £3 


an 


oO 
cae 
a 


SURR 


INPU 


gure 


cn 
OFFT 


DATA 
DATA 


pT2 
SAAR 
RET 
END 


OUTTNE SOLECYCCDAYI) 


SUBROUTINE SOLCYC COMPUTES THF SULAR FLUX SRAR, AN INPUT TH) 
ATMNSU THRELIGH COMMON ATMMUP, BASFD UN AN ASSUMED SINUSMI MAL 
11e¥R (NKR GO{ReDAY) VARKTATION, WITH THE MAXIMUM VALE OF 256 
FOR SRAR, ASSOCTATED WITH CIRA®HS MODEL 9, OCCURRING ON 

195R JUNE 1, THE MINTMUIM VALUE OF 65 FOR SRAR TS ASSNCTATEN 
WITH CITRA#ES MUDEL 1, 


T PARAMFTER 
NAYS @ JULTAN DAY NUMBER (A HALF INTEGER) AT 0 HRS UT 
UN THE DAY OF INTFREST, 


UT PARAMETER 
SHAR ew AVERAGE 10,70eCM SOLAR FLUX, 
SRAR TS AN TNPUT TU ATMNMSII 


1,0Fe2e w/(Maxe H7), 
THRUUGH CUMMIIN ATMOUP, 
UNJATMQUP/S HL, SBAR,IDURN,PP,RAD,TT,SNI (16) pHRHN,FEHSED 
NITION OF DATA 
NJ5R806 6 JILLIAN DAY NUMBER ON 
NISAN& 4 2436555,5 / 
PT / 3,141592653590 / 


1958 JUNE 1 = 2436355,5 


= eet] 
= 157.5 * 92,58*CNS( 
RN 


(DAavyJePJSB06)ePT 2/4018, ) 
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ORIVER, 1331 
DRIVER,1 $32 
ORIVER, 1333 
DRIVER, 1554 
ORIVER, 1335 
URIVER,13306 
ORT vEe®,1337 
DRIVER,1338 
ORI VER,1339 
ORIVER,1340 
DRIVER, 1 $41 
URIVER,1 $42 
DRIVFR,1343 
ORTVER,1 344 
DRIVER, 1345 
DRIVER, 1346 
VDRIVER, 1347 
ORIVER, 1348 
VRIVER,1349 
ORI vER,1550 
DRIVER, 1551 
ORI VER,13552 
DRIVER,1353 
URIVER,1354 
DRIVER,1355 
URIVER, 1356 
ORT VER, 1357 
VDRIVER,135H 


SOLORB 


SUBROUTINE SULURKOCYREY,VEOI,DaYI,SULLAT,SOLL UN) 


c 
on) 
OO 


SUBS JLAR PHTNT, 
JaniiaRy 4 OF TRE YEAR QF INTEREST 


aHI CH VERNAL 
O mRS UN THE 
TYME (UT), 


FQUINOXY OCCURS (VERS), 


NAY OF INTEREST 


cc 


VEQ) « 
payg 


ut 


INPIIT PARAMETERS 


VRE J © JIL TAN TAY NUMBER CA HALF INTEGFR) AT 0 RRS UT ON 


(LAY), 


SUARJUTINE S)LORR COMPUTES THE NUPTNH LATITUDE SLL AT ANN 
EAST LONGTTUDE SOL LON UF THE APPARENT 
GIVFN THE JUL TAN DAY NUMBFH AT 0 HRS UT ON 
(YRFJ), THE 


CACTUAL MOTTUN) 


JULTAN PATE aT 


THE JULTAN DAY NUMBER aT 


AND THE 


UNTVFRSal 


JANUARY 1 UF THE YEA OF INTEREST, 


JULIAN DATE FOR VERNOL FUUINUK, 


© JULTAN DAY NUMBER (A HALF INTEGER) af 0 WRG UT 
UN THE DAY NF INTFREST, 


e UNTVERSal TIME 


OUTPUT PARAMETERS 
GAT 


(CNECI Mac 


© GREENWICH APPARENT TIME 
GAT 1S PLACED JN CUMMON TIME, 


SOLLAT © NORTH LATITUDE 


HRS), 


(NEC IMAL 


OF SUBSULAR POINT 


SULLON © EAST LONGITUDE OF SUBSOLAR PUTNT 


DFFINTTIONS ANN COMMENTS 


HRS), 


(KADTANS), 
(RADTANS), 


UTO24G TS TRE DECIMAL FRACTION UF DAY CURRESPUNDING TO UT, 

Oaytut 18 THE SULTAN COFCIMA@L) DAY NUMBER AT UT HRS UN THE 
bay OF INTEREST, 
DAYNY) IS TRE NUMBER CIF FLAPSED CDECTMAL) DAYS SINCE THE 


BEGINNING OF THE VEAR AT 0 HRS 
(DaYJUT © ABINT(CDAYIUT)), 


THE QUANTTTY 


UT ON JANUARY 
THE WEST LUNGTTINE OF 


', 


THE SURS)L AR PUINT FXKPRESSED AS A DECIMAL FRACTION (1F pePT 
RANTANS, TS SURTRACTEN FROM 


LONGITUNE, 
LONGI UNE NF 


TO OBTAIN THE FRACTYUNAI 
THE FAST 
CFICTIVTINUS) MFAN SUN, 


JF THE 


THE FIRST Twi) EXPRESSTONS FOR SOLLUN ARE 
THE SURSOLAR POINT 


FAST 


IT TS PASSIBLE TU MAKE AN APPRUXI MATE CORRECTIUN FUR THE 


OYFFERENCE BETWEEN THE APPARENT 


AND THE MEAN 


CACTUAL MOTTUN) SNLAR TIME 
SOLAR TIME, KNOWN AS THE EQUATTIUNO@IF eT IME 


(SEE, 


E.Gee AMERICAN PRACTICAL NAVIGATUR (URIGINALLY AY N, 
S, NAVY HM, Pllib, NO, 9, P, 375, 


BOWNITCH), VU, 
CORRECTFD REPRINT EDITION, 
IN THE U1,9,4, 


FE TCE), 


UF 19be 


AVAILARLFE FROM USS, GOV, PRINTING 
CIN CONTRAST TU GREAT BRITATNY TRE 


STUN OF THE FQUATTUNSOFSTIME TS CONSIDERED PUISITIVE JF THE 


TIME OF THE 


HRS ANI) NEGATIVE TF LATER THAN 1209 HRS, 


MERINYTAN TWHANSIT RY THE SIIN IS FARLTER THAN 1200 
(NOTE THAT A 


MERTNTAN TRANSIT REFURE 12009 HRS FURRESPUNDS TO TRE East 
THE SUN RETNG SMALLER THAN THE VALIJE EXPFCTEN 


LONGITUDE nF 


BASF.) ON & MEAN SUN,) 


ANNUAL EDITIONS 


OF THE NAUTICAL 


VRIVER,1359 
URI VER,1 560 
OR IVER, 1561 
URIVER, 1362 
ORI VER, 1565 
URIVER,! 364 
DRIVER, 1365 
URIVER,1 S46 
ORIVER,1 $47 
URIVER, 15466 
VRIVER, 15469 
OR IVER, 1370 
ORIVER 1371 
ORIVER,1572 
ORIVEP 1373 
DRIVER, 1374 
ORIVER,1375 
ORIVER,1$76 
ORIVER,1377 
VRIVER,1378 
ORIVER, 1379 
DRIVER,1380 
DRIVER,I3R1 
ORI VER, 1 3R2 
ORIVER,135R3 
OPIVER,135R4 
ORITVER ISAS 
URIVER,1 386 
DRIVER, 13R7 
ORIVER,13R8 
LP1VER,I3R9 
VRIVER,1390 
DRIVER, 1391 
URIVER,1392 
ORIVER, 1393 
ORIVER, 1394 
UDRIVER,1 395 
UVRIVER {13% 
ORIVER,1397 
ORI VER, 1398 
OF I1VER,1399 
DRIVER, Lunn 
URIVER, 1491 
DRIVER, 1402 
URIVER,14N3 
URIVER,1unu 
ORI VER, 1405 
DRIVER, (U6 


ALMANAC PRITIR 11 4962 TABULATED VALUES OF THE EQUATTON@SD FeTIMFDRIVER, TUN? 


AT 120eHR TNTERVALS, 


TIMF, 
NAUTICAI 


AMERICAN 


FPHEMERITS AND 


INSTE AN, FUR 


MERTIYTAN THANSITS AND UTHER PHENUMEND THAT DEPEN() UN WLR 


(OR UNIVERSAL 
CUR ACTUAL MOTTON) 


THESE TARULATEN VALUFS UF THE FULTATION®UFORIVER,1U0R 
eTIME COULD BE ADNEND TO THE GREENWICH MFAN TIME 
TIME) TO ORTATIN THE GREFNWICH APPAKENT 
NEWER ANNUAL EDITIONS MF THE 
ALMANAC OR THE ASTRONUMICAL FRPHENMERIS Nii NOT EVEN 
EWPLICITLY REFER TU THE TERM FUUATIONeUF eTIME, 


DRIVER, 1u09 
ORIVER, 1410 
URIVER,1utY 
URIVER, 1412 
URIve® surg 
ORIVER, 1414 


qf 
{2 
1% 


ANGLES AND GFOGRAPHIC LOCATION, 
TO THE GREENWICH MERTOIAN AND JO UNIVERSAL TIME BUT TO & 
MFRTOTAN 1. 90273SRe (DELTA T) EAST OF THE GFUGRAPHIC MERION AN 
OF GREENWICH CKNOWN AS THE FPHEMER TS MERINIAN) AND TH 
EPHFMFRTS TIMF, THF SULAR FRPHEMERTS TRANSIT, WHICH IS THE 
EPHEMFRTS TIMF AY THE INSTANT UF SOLAR TRANSIT ACRUSS THE 
EPHFMFRTS MFERTOTAN, IS TaRULATED AT {DAY INTERVALS IN THE 
NEWER ENITIONS, WE HAVE ANOPTED THE NEPARTURE OF 
THE SNLAR EPHEMFRIS TRANSIT FROM $2 HR M0 MIN ON SEC 4S A 
CONVENTENT APPROXIMATION TQ THE NEGATIVE VALUF OF THE 
EQUATTON@OFeTIME, TN PARTICULAR, WF HAVE USED VALUFS OF 
SOLAR FPHFMERTS JRANSTT FOR 1974 TARULATED TN 
VF FITHER THE ASTRUNUMICAL 
AND NAUTICAL ALMANAC, AND FITTED TUR ADOPTED VALUFS UF THF 
EQUATTUN®ORe TIME RY A FOURSTERM FOURIER SERIES, 
WEAK DEPENDENCE OF THE FQUATION®UFeTIME UN THE YEAR OF 
INTEREST, LOWENIS FITTED EXPRESSTUN FUR THE FQUATION@(OF eT IME 
IS GIVEN Ry 


THE NEWER EDITINNG REFER NUT 


TRE 


EQT & O,3A517SHCUS(F) @ 5,146125*#CNS8(Fe) 
= T.392NSSRSIN(F) © Ge5SSORASHSIN(F2) 4» MIN 
WHERE 
F 
F2 
RANDAY 


RADDAYR (OAV Se VRE S) 
Ped 
20%P11565e25 RADIANS PER NAY 
O,0tTAdeuesR , 
TA CONVERT FROM MINUTES UF TIME TD RADIANS SF LONGITUDE we 
MUST MULTIPLY EOT BY 

RANMIN = 2y*PL/{[440 RAPT ANS PER MINUTE 
= 0,0%4%63327515 , 
FAST LUNGTTUDE (RADTANS) 
2 SOLLON@RADMINKEDT 
LATITUDE CRADTANS) OF THE APPARENT SUN 1S 
= SLATMY*SIN( (MAYJUT#VE QS) #RADNAY J 
MAXIMUM VALUE QF THF SULAR LATITUDE I8 
= 0.409123 RANIANS , 


THUS, THE 
SOLAN 
THE NORTH 
SOLLAT 
WHERE THE 
SLATMx 


OF THE aPPARFNT SUN JS 


7 
Cc 


COMMON/STIMES JYRS,IMONS,TOAYS,Z2T,PLAT,PLOAN,UT, GAT 


ccc 
c DFFINTTTIUNS OF MATA AND CONSTANTS 
€ PY = $,141592653599 
C PIl2 = 2e*Pl 
Cc WANDAY = PT2/365,25 RANTANS PER PaY IN A JIL TAN YEAR 
C = 0,017202u25A 
C RANMIN = PT2/1440 RANTANS PER MINUTE IN A DAY 
C = 0,00436352313 
C SLATMK 3 MAXTMIIM VALIIE OF SULAR LATITUDE 
C = 90,409123 RANTANS 
crc 
OATA PTZSLATMK 4 $464415926053590, Ued4NG123 / 
cece 


Ple = A,#P] 
PT?P/465,°5 
es PyPysinun, 
UT/2PU, 


RADMNAY = 
RANMIN 
WTDAW & 


ORIVER, 1415 
DRIVER 1416 
ORIVER 1417 
URIVER, tata 
ORIVER, 1419 
DRIVER, 1420 
DRIVER, $42} 


THE VALIIE OFDRIVER, 1422 


DRIVER, 1423 
ORIVER,{u2u 
DRIVER 1425 


THE 1974 EDYTTUNORIVER, 1426 
EPHEMERTS OR THE AMERTTAN FRHEMERTSORI VER. {427 


ORIVER, 1428 


WE IGNORE THEORIVER {429 


ORTVER, 1430 
DRIVER 1431 
ORIVER, 1432 
URTVER, 1433 
ORIVER, 1434 
ORIVER, 1435 
DRIVER, 1436 
ORIVER, 1437 
ORIVER, 1448 
ORIVER, {439 
ORIVER, 1440 
ORIVER, $444 
DRIVER, 1442 
ORIVER, Suus 
DRIVER, 1444 
DRIVER, 14u5 
DRIVER, 1446 
VRIVER, 1447 
DRIVER, 1446 
URIVER, t4a9 
ORIVER, 1450 
VRIVER, 145} 
DRIVER, f45e2 
DRIVER, 1453 
DRIVER, 1454 
DRIVER, YUSS 
DRIVER, 1456 
ORIVER, SUS? 
ORIVER, 1458 
DRIVER, 1459 
DRIVER, 1460 
DRIVER, 1461 
ORIVER, 1462 
DRIVER, 1463 
UVRIVER,1uAu 
DRIVER, 1465 
ORIVER, 1466 
ORIVER, 1467 
DRIVER, 14668 
DRIVER, 1469 
DRIVER, 1470 


DAYJUT s Pavy «+ UTD24 

NAYNO s PAYIUT © YREJ 

SOLLUN 3 PlJe(y OPDAVIUTeAINT( DAYJUT )) 

TFC SOLE ONAL T,0,0 2 SULLON & SALLUNOPYDA 

F = RADNAYRN AYN 

FR = A aF 

EQT = O,SR5175RCUS(F) © §,146125"CUS(F2) 
© 7,592635RSIN(F) @ 9,536R25*SIN( FQ) 

GAT = UT # FUTs/A0, 

SOLLON & SOLLON @ RADMIN#EUT 


SOLLAT = SLATMX®SING (DAVIUTeVEOJI*RANDAY ) 


RE TIIRN 
END 
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DRIVER, 147} 
ORIVER,1472 
ORIVER,1G7$ 
ORIVER, 1474 
ORT VER, 1475 
DRIVER, 1476 
ORIVER, 1477 
DRIVER, 1476 
ORIVER,1479 
DRIVER,14R0 
ORIVER, 14R} 
DRIVER, 1482 
ORIVER, {UR 


Ai lt 


SE a a aT es 


125 
127 
147 
fan 


tur 
14a 


1S4 
158 


SOLVE 


SUBROUTINE SOLVE (ay X, NM) 


VRIVER, 1uURG 
ORITVER,IUAS 
ORI VER, 14GR6 
ORTVER {UAT 


DRIVER, 1490 
ORI VER, 14914 
ORIVER, 1492 


ORI VFR, 1494 
ORT VER, 1495 
ORI VER, 1496 
DRIVER,1497 
DRIVER, 149R 
DRIVER,1499 


ccc 

t SBR WTINE SaLVF, CALLED FROM SlisRyiTINE FITTER, SOLVES A SFT 

c HF ONG STIMU TANEINIS LINEAR ALGFEBRATC EQUATTUNS &Y TISTNG 

€ GAUSSeJOKRDAN METHOD wITH MAXIMUM PIVIJT FEATURE, (SEF, FORTRAN URIVER,1URR 
C Ty PRIGHAMYING AND COMPUTING BY JAMES Ty GOLDEN, PRENTICE SHALL e DRIVER, IURG 
Cc INC ey PAGES BAH99D) 1965 

cee 

C INPHIT PARAMETERS 

C ACT,.1) © MATRI¥ '1F CONSTANT CUFFFICIENTS IN SET CONTAINING URIVER, 1493 
C THE NUIMRER Ni) STMULTANENUS LINEAR ALGEBRAIC 

t ENUATTONS 

C NO) we THE NIIMBER OE EQUATTONS 

cee 

(2 INITPUT PARAMETERS 

( ¥(K) © THE LEAST@ESQUARFS FIT COEFFICTENTS 

CCl 


DIMENSION ACAD,2A1)6 ROA0021)e KC AO), LNLC20), Riwl?0) 
KA) S NMe4 
on 150 Tet,N0 
ON 190 Jst,Knn 
ROET,J) = ACTS) 
180 CONTINUE 
O01 10 “S1,N9 
LOCfMy =n 
10 ROw(M) = 0,9 
NP = NiO] 
on 100 Ts1,N0 
TP = Tet 
CeeeeeF IND MAY FLEMFNT IN ToT CUL, 
AMAY = 9,9 
DN 2 Kst,NN 
IF(AMAX © ARSE A(K,T))) $e2e2 
Ceeeere |S NEW MAY TN ROW PREVIQUSLY I'SFOD AS PYVIT, 
BS LECRIWO(K)) 4,4,2 
oLOCCIY = K 
AMAY = AHS( ACk,T1)) 
2 CONTINUE 
IF CAMAX) 99,99,9R 
CoseweMAX ELEMENT IN TeTH COL IS ACL,7) 
o8 | = LACcer) 
Wow(L) ££ 1,9 
CwwemeP ERE RM FL TMINATION, | TS PTVAT KOM, ACL ZT) IS PIVGT ELEMENT, 
NN SQ JSi_,nn 
IF(LeJ3) 4,459,464 
eo WF = eACT,TY/ZACI,T) 
On do KE[P,NP 
ACJoK) & ACT KY) OUFMACL KK) 
40 CONTINUE 
BO CONTINUE 
180 CONTINUF 
on AoA Tst,NU 
te unerr) 
POO eF1Y s&s ACL OMIT IZA(LT) 
fa MAPITTE (661089 CF, MOTY, Set enid 
( 193 FORMAT (40 TA, AX,FIS,AD) 
HE TIIWS 
99 wRITE(H,1N4) 
VAa FORMAT (50 ,A7H Si) UNTONIF SOLUTTON eEXTSTS,) 
WE TIA 
FND 
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DRIVER,1590 
UDRIVER,1591 
ORIVER,1592 
ORITVER,1593 
ORIVER,1S94G 
URIVER, 1505 
ORIVER,1506 
ORI VER,1507 
DRIVER, 1508 
DRIVER,1599 
DRIVER,1510 
ORIVER, 15114 
DORIVER,ISte2 
DRIVER,ISIS 
ORIVER,IS14 
ORIVER,1S15 
VR1ivEe®, 1516 
URIVER,I517 
ORIVER,1518 
DRIVER,1519 
ORIVER,1500 
ORIVER,1SA4 
URIVER 1522 
ORIVER, 1523 
DRIVER,1$2u 
DRIVER, 1525 
DRIVER, IS26 
ORIVER 1527 
DPIVER,1S2R 
DRIVER, 1509 
ORIVER,IS30 
DRIVER,IS34 
URIVER,IS382 
ORIVER,I1S3$ 
DRIVER, IS3G 
OR TVER,SIS3S 
VRIVER1S36 
UR LvER,IS37 
URT VER, ISSA 
DRIVER, ISS9 
URTVER, 1540 
USI VER, IS4L 
VRLVER,1Sue 
URITVEE.IS4§ 
ORL vER, 1544 


e? 
30 
34 
36 
37 
4? 
are 
a? 


SOLZEN 


SUHBMIITINE SOL ZENCSOLE AT, SO LINY 


ere 
C SUMP UITINE Sop ZEN COMPUTES FC ISCHE, THF CiSINE GF THE ZEN TT es 
[3 ANGER (IE Tee SUN AT A POINT #, GIVES) THE RENGRAP HTC NORTH 
C LATTTUMFE PLAT AND EAST LONGTTHOF PLON UF THE PUTINT P AND THE 
C NORTH CATTTINE Sut Lat aN EAST LUNGTTUDE SOLLON UF TRE 
Cc SlisS WAR PLINT, TRE NAVSRENTONT PARAMETER INURN JS 4 FUR 
C YAYTIMES TF oe TEECOSCHYE GELO.0), AND TS of FUR NIGHTTIME, 
C Lefer ITECTUSCHTeLTAO,97, THE LOCAL APPARENT TIME HL 
c TS ALS!) COMPUTED FROM THE GREENWICH APPARENT TIME GAT AND 
C LONGITUDE FLUN, 
Cee 
c INPUT PARAMETERS 
C PLAT e NORTH LATITUDE HF POINT PB (RADIANS) 
C PLON © EAST LONGITHDFE OF POINT P (CRANILANS) 
(3 Su Cat we NORTH CATTTHDE FF sUiBSULAR PUTNT (KANTANS) 
C SULLON w BAST LOUNGTTIINE OF SilbSiiL AR POTNT CKANDTANS) 
GEC 
& GTP UT PAR AMETP RS 
C TMIRN @ PARAMFTFER FOR DAY UR NIGHT, TF COSTRY TS 
C THE CNSTwk NF TRE ZENTYTH ANGLE (1F THE SItIN aT 
C POINT P, TDOMRN TS 1 FOR DAYTIME, T,Fy, 
(a TECCOSCHE GF ,9,0), AND TS ef FUR NIGHTTIME, 
is fetes IF (CUSCHITLTe%,9) » ITDMRKN TS AN TNPUT TH) 
a ATMOSH THROUGH CUMMON ATMOUP, 
U MHL oe LOCAL APPARENT TIME (DECIMAL HRS, &€,6, 2230 HRS 
€ aHFLOMES 22,50 HRS), HL IS AN INPUT TA ATMOS 
C THROURH COMMIIN ATMOUP, 
GEE 
CEM SAP MIIPS HL ep SRAR eS TDURN APP, PRI, TT,SNI C165) pARHO,FEHSEN 
COMM ONAT IME Z TYRS, TIMINGS, TIAYS, ZT PLAT, PLING UT ee GAT 
NATA PT 4 B,14ISAAASTSIN / 
Cec 
é THE FOLLOWING FORMULA [TS RASED nN EQ, (1,41) MF IANOSPHERTC 
€ RADT) PROPAGATION BY Ky DAVTESe NAS MIINIGRAPH BO, 1965 
Cc APHTL 1, J TT MAY ALSU tF DEPTVED Ry APPLYING YHE Law MF 
(s COSINE S FUR AN URE TOQUE SPHERICAL TRIANGLE, 
Lo ai 
COSCHT = SIN(PL AT) & SIN(SHLLAT) 
* * CUSCPLATY * FHOSCSHLLAT)Y * CUSCELON®SALLUN) 


IMURN = 4 


TEC COSCHT, 1 Te9,.0 2 3 PRURN & eI[NIRN 
PTS & Ayers 

HWALHR = PTI, 

Hl = GAT © (PTD @PLIIN) JR ATH 

TFG mt LF He ¥ BL & BES2G, 

RE TIIWN 

ENE 
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THE 


DRIVER, SSeS 
DRIVER, 1SUm 
ORT VER, IS47 
URIVER, {54B 
DRIVER, IS49 
URIVER, {550 
ORT VER, 1551 
URLVER,15S2 
URIVER, 1553 
UAITVER,15S4 
URIVER, 1555 
URIVER 1556 
DRIVER, 1557 
DPT VER, 1558 
URIVER,15S59 
VRIVER,156u 
DRIVER, 1564 
UR] VER, SS42 
ORIVER,1S6$ 
URTVER,1S64 
URTVER 1565 
DRIVER, 1Sh46 
URIVER,15467 
ORIVER, 1568 
DRIVFR,1569 
URIVER, 1570 
DRIVER, 1574 
DRIVER 1572 
URIVER, S573 
DRIVER I57u 
URITVERISTS 
DRIVER,1S7%6 
ORT VER,1577 
ORI VER,1578 
URI VER 1879 
URIVER, 1580 
DRIVER, 1581 
DRIVER, 15R2 
DRIVER, SSA3 
DRIVER,15Ru 
DRIVER,ISAS 
DRIVER, 1586 
ORIVER, LSAT 
ORIVER,15AR 
DRIVER,ISRO 
ORIVER, 1590 
DRIVER, 1594 
VR IVER, 1592 


ANAANAANMNANDNAMAAAANAANAANANAAANAAANAANANRANRAAAN 
a 
OQ 


°° 
oO 


ccc 


* 
* 


SPCMIN 


SIIBROLUTINEG SPCMINOKK, 7H) 


THE HT GHewAL 
SPECIES 0, 

ALTYITUDFS a 
SPEMIN wE P 
THE LOwealT 
cn2, N, AND 
OELTA G), 0 


TITUDE CHEMISTRY MODIILFE REQUIRFS THE MINOR NEUTRAL 
CU2, N, AND NO, PRUFILES FUR Day AND NIGHT AT ALL 

RE PROVIDED FOR 9 AND CUA IN ATMOSU, HERE IN 
ROVINE PROFILES MF N ANiy NO, 

TTIDE CREMISTRY MUNULE RFGUIRES, IN ANDITION 19 1), 
NO, THE MINOR NEUTRAL SPECIES 2), NOCSTINGLET 

3, AND NOG» ALSO PRUVINEN BY SPCMIN, 


INPIIT PARAMETERS 


ARGUMENT Lt 


sT 


KK @ CALCULATION FLAG 
3 1, CALCULATE INTTTALI7ATION PARAMETERS 
= 2p CALCULATE ATMUSPHERIC PROPERTIES 
TH w ALTITIDE OF INTEREST (kK) 
ATMAUP COMMON 
TNORN © INDEX FOR DAY OR NIGHT 
sol, bay 
= #1, NIGHT 
ALTODN CUMMON 
$1Z2N 
UUTPUT PARAMETERS 
ATMNUP CLIMMON 
SNI¢C 7) & N OENSITY, 1/CMaegy 
SNI( A) « NO DENS] Ty, L/CMae§ 
SNI(13) © O2(8NG) HENSITY, I/CMeed 
SNI(14) © O83 OENS{TY¥, 1/C¥a«3 
SNI(15) © NN2 DENSITY, L/CMxed 
SNIC16) © #2) DENSITY, L/CMaes 
ALTODN COMMON 
NALTOD » NUMBER OF ALTITUOFS AT WHICH THE DAYTIME 


ALTKM(47) « 


NDAYC27) « 
ONITEC{A) «© 


CN2025) «@ 
DIMENSION AAC] 
DIMENSIAN NSN 
DIMENSION V(b) 
DIMFEFN@GION Hatin 
DIMFNSION ae20, 
DIMENSION SNe 


COMMONSATM)IIPS 
COMMON ZALTODNS 
COMMON/ZHOCHEX/ 


DATA 
OaTa 
DATA 
DaTa 
OaTa 
DATA 


NDEGNN / 


NALTON,NALING 74 27,25 s, 


OeVALIieS ARF SPFCTFIEN AS DATA 

THE ALTITUDES AT WHICH MINUR SPFCTES ARF 
SPECIFIED 48 DATA 

THE DAYTIME feVALIIES SPFCIFIED 4&5 DATA 
THE NIGHTTIME OevALilES SPECTFTED AS DATA 
THE CO2eVALUES SPFCTFYIEN AS DATA 


3), B8R613),COC13), ANUNIT(18) ,ANDAY (47), ANNITE CUT) 
60647), 025S9GN(47),9SNAV( 27) -O3NIT(27),90011) 

eZ lod ,FECIINI,UN3S (Oe VU3(b) WOS(0) 

N(€25),66(13) ,ANODAY(4S) 

PLEX (9) Fe ZINONOBID,ANINZI (BY ¢Z2TM2NU(8) 
NC33),S8N02N(33),HHE13) 

HL» SBAR,IDOURN,PP,RHO,TT,SNI(16),HRHM,FENSFO 
NALTOM,ALTKM(47),904Y (27) ONT TE (1A), 812AN,CNI2625) 
ZHFLAG 


6 / 
NALTNO,NALTNN 4 39,59 / 


12 /¢ NDEGND,NDEGNN / BR, 


NA9GOeN,NALTN2 4 19,11 / 


NOGH20,NK 
NIGNOA,NK 
CALTKM(T) 


MHOU / 12,25 / 
MNNQ / 12,33 / 
el2ie47) / Me rPSe eter IVS es 20 ge Dqn SU gn 89g pH ene USer 


ORIVER, 1593 
URIVER, 1594 
DRIVER 1595 
ORIVER 1596 
VRIVER, 1597 
DRIVER,1598 
VRIVER,1899 
DRIVER, 1690 
DRIVER, 169} 
ORIVER, 1692 
URIVER, 1603 
URIVER, 1604 
DRIVER,1695 
URIVER, 1606 
URIVER, 1607 
URIVER, 1608 
ORIVER, 1609 
URIVER 1610 
URIVER, 1611 
ORIVER,f612 
DRIVER, 1613 
ORIVER toiu 
DRIVER, 1615 
DRIVER, 1616 
VRIVER, 1617 
VRIVER 1618 
ORIVER, tog? 
ORIVER, 1620 
DRIVER, S62} 
DRIVER, 1622 
ORIVER, 1623 
DRIVER, 1624 
UDRIVER, 16005 
URIVER, 1626 
DRIVER, 1627 
DRIVER, 1628 
ORIVER, 1629 
DRIVER, 1630 
ORIVER, 1631 
DRIVER, 16%2 
ORIVER, 1633 
DRIVER, 1634 
URIVER, 1635 
ORIVER «16% 
VRIVER, 1637 
URIVER, 1638 
DRIVER, 1639 
URIVER, 1640 
ORIVER, 1644 
URIVER, 1642 
ORIVER, 1643 
DRIVER, 1644 
ORIVER, 1645 
DRIVER, 1646 


5 055g OD gOS ge ge Tg Bg BS gp GO gp Me HP VRIVER, LOU? 
LDV gr NOS LIND eV ge LAD gr I 25g ol SV gH LST e el UD PINS 150g pI S508 


DRIVER, 1648 


Nan 


oO 


L460 ep LOS get Meal Mg pl AO ee BS gpl Wag IWS gr QO0 ep 20S pelle, eles 
2200225002300 / 


BFM VALIIES 41/09/74 RNR On nay 


DATA (ODAYCTI,TH1,?7) 4 L,O0F 15,0, 2F HO 5,2 ,2E HIN 1 GF HOS, Lp OF HNO, 
* TOE HOON SFHNT 2p OF HOB I BE HNR 2, 9F HOF, 
B BL TFHOGL  SFalO, eg dF ola SIF H10,e BE e10,4 SEGIOe TIF +1 O,1,1F oily 
ee eae ee ee ee ee oe ee ee - 1e e e 
# $§,7E&+10 s/ 

BFe VALUES 02/22/75 FOR O NIGHT 
VATA CONTTECT),121019) 4 $38tele 20,0, 4,908 408, 
* 3,00F +10, 9,90F410 / 
BFM vALHleS 1As07/74 FOR CN2 
DATA (CU2CTI»TE1625) 4 21#A WO, Ly SOE 09,4 AOEFOB,L TOE SOR, 


* S,OS5E+07 / 
THE CQ VALVES AT ALTITIDES FRUM a,o TU {90, KM ARE RESET 
IN) SUBROUTINE ATMOMSU BY USING A CONSTANT MIXINGeRATIO UF 


3,20E~04 

BFM VALUES 04/95/75 FOR NO NAY 
Nats CANDDAYC TY, TE1,25) 4 1, 00F 410535, 496 409,1, SOF ENG, S BNEFOR, 
* 7, 00F eNBr lL, PSEFV9, 2,1 0F 49,1, TSE F09, 


# 1,25F 09,8 SNE40A,5, 106408, 3,096 4097) HOE HOB, S,SHEFOTD 3, IDF HOT, 
wR BLVOFSOT, S SOE SOT, Y NOE 407, Ue BNE S07, 5 AOE HOT 6, 3S0EF0 74S, TOF FONT, 
et A,MOE HN, 8,H0F 407 S,90F O97 


BFM VALUES y4/05/75 FOR NU NIGHT 
DATA CANONTTCYS,129,98) 4 Viel, AOFHOO,1,00E 04,1, N0F eb, 
* Be S9EFUG eS OSE SIT SSE FOIA 3, 30E 407, 
* UOOF ENT / 

HEM VALIIES NU/12/75 FOR ON NAY 


DATA CANDAY(T),TEt,87) 4 PHL, ODE S00, NOE FOL, LODE S025, NDF eNe, 
RL ROF SNS TUNE 403,21 0E F945, ANE 40441 1 0E F952, 20E 405,43, 70E 405, 
® HYDE HNL ONE MIA, 1 SOF HDG, 2, 0NE SDH A, TOF SO S,4NE SOA, SOE +06, 
ST SMF HMO S, IME HOH pA SOF HCO, T ANNE Hb r 7 SOE HOO, TIDE FOAL A, I OF 406, 
® ALOE ENGR UDF 406A, SOE O06, Bg SOE SOKA UOF O06 LB, SOE SIAL A, POF 406, 
* 


RLV0F Nbr BL ODE Fb, 7 BOF EDO, TeSDE SUH 7, SOE SI be Tg INEFOh eH, ROE SIG / 


BEM VALUES 04/12/75 FOR ON NIGHT 
Oata CANNITECTI, TSI, U7ISIBHL DOE SON 1 AOE SNL rl eeQNEFO2e7F NOE eN2, 
# PLVOF HOS, LMF O04, S90 40K, 7, S0E 04,1  F0F 495, 3, 00EF05 24, NOE SNS, 
BD ROF 4956 5. H0F O05 45, SOF 495 OBOE HIS TePVE HOS Tg SNE F054 7 eROF ENS, 
# VIDE HONS TIDE 4054 7 ROOF O05, Ty TOE US, 7 S0E F095, 7, 20E 40576, 9NF HS, 
® HAO ENS 6, SDE 4S HOOF ONS, SHOE SUS DS, LUE FOS / 
HFM VALIIES U1/040/75 FOR N2(SNAG) DAY 
CQASNGDOT),1T21,47) 4 AqgH0F eNO, GU HNE FUG, Ae TOF ONT AL 25ESOR, 
Gp FOF HMB oe 1 ~25EFOI, Ag TOE HN9, 9, ONE HU9, 
LeROF OND Ae TIE HID, S.SOESTO AINE HILO 1 SOF HID gL ONE FING HB, LOE O09, 
FeV OF HOM ASLOSE HAI, SF, HOF HOG 1 gp STE FOF B~NVEFIOB re SeODE FOIE, SOF O06, 
He PDF ONS LOVE FOS e I AOF OOK SSE FOS, I LOE FNS, A ONE FO 2g t OOF eNe, 
H,TOF ODL 2, SME FOS 1 APOE SOL, 15eh tu / 
BEM VALUES ul/0as7S FOR N2CSNAG) NIGHT 
DATA (TIASHGNCTI,TB1,407) 4 15% 3,40, AOF ONS ef ,00EF0S,AgHOF ONT, 
* 2, TNF HOB sr 1 UNE 408, SHOE 407, 4, SNE SOG, 
# HePOFHNSHs 1 DOE FOS, 1 HOF EIU, 3, STE O05, TIDE PNA, CONE FOAL 1 NOE HN, 
P ULTUF SNL ASIF SOT, LAGE SIL, 1SaHy10 &/ 
BFM VALI'IES O1/18/75 FoR 5 OAyY 
CNSDOAV(TS, LBL ,27) £ A SEL TeTEDL 1 FL ILE FI 2,2 SEM12, 


JATA 


sr? + + © 


Haya 


E BEF 1 Arg SEH1L2p 2, SEH 2,1, dE SIP, 
* HeAEFL ACME SLD Hy TEFIO, ASE FIN, Te UEHUIP CIE FNI,S, SEFOR, 


DRIVER, 1649 
DRIVER, 1650 
DRIVER, 1651 
ORIVER,1652 
DRIVER,1653 
ORIVER, 1654 
ORIVER, 1655 
ORIVER, 1656 
URIVER, 1687 
DRIVER, 1654 
ORIVER,1659 
URIVER, 1660 
ORIVER Leet 
VRIVER, tobe 
ORIVER, 1663 
DRIVER, 1664 
DRIVER, 1665 
DRIVER, 1666 
DRIVER, 1667 
DRIVER, 1668 
OPIVER, 1669 
DRIVER, 1670 
ORIVER, 1071 
DFIVER, S072 
DRIVER, 1673 
VRIVER, 167%u 
DRIVER, 1675 
DRIVER, 1676 
ORIVER, 1677 
DRIVER, 1678 
DRIVER, 1679 
DRIVER, 1640 
DRIVER, 168} 
DRIVER, 1682 
ORIVER, 1683 
ORIVER, S6Au 
DRIVER, 1685 
VRIVER. 1686 
ORIVER,16R7 
DRIVER, 16846 
ORIVER, 1689 
DRIVER, 1690 
DRIVER, 1694 
ORT VER, 1692 
DRIVER, 1693 
VRIVER, 1694 
DRIVER, 1695 
ORT VER, 1696 
UDRIVER, 1697 
ORIVER, 1698 
DRIVER, 1699 
URIVER 1790 
DRIVER, 1701 
DRIVER, 1702 
URIVER,1703 
DRIVER,1704G 


. 1, SESOR, A AFOOT, 1 PEHOR, SDE FO Te 2, 4E 0G, 2, 5E405,3,4E 404, ORIVER,1795 


* BL 7EHOS, SRE FOA, 4, PESOTL 7, 3EH0NE OGD S DRIVER, 1706 
Cc BFM VALIIES O1/18/75 FOR 3 NIGHT ORIVER, 1707 
DATA CABNTTCT), 121-27) 4 AOE 41165, 7E411,1,1E412,2,5E 612, URIVER,179% 
* GREP L AG e SEHL2Ae2eSEFL2yt edt Ol, DRIVER,1799 
* HeNFHLL AME SL yp Me FESIN AME HIM le SES1ND 1 PEIN, UY, SEONG, VRIVER,1710 | 
* Tae TEP OR GAVE HOT Ee SAE HVA 1, TEHOR Oe PE HUT 2g PEMOAD Ag PEHNS, ORIVER, 1791 
* GMEFO4, SSE OOS, FESOA, A UE HOL OME SON / DRIVER,1712 
C HFM VALUES 05/10/75 FOR He DRIVER,IT713 
DATA (HANDN(T),Iele2S) 4S 1,961 %e2etEo16r6,9EO14,1,28 41%, ORIVER,I714 
* 5.261245, 2641261, ALL, b, PEO, VRIVER 1715 
* He VECLP, APF ot eB GE O19, SeHEF {Me Le TEF10,%, PE 409, 4, 4E 49, URIVER,17to 
* LAE FOI, ASE SOR, AME FOR, 4, FEO TAG UE HOGS eLEFOHr 1, SE%0S, URIVER,1717 
* 2,5F 04,8, 764d, 3, 8E Hd / ORIVER, 1718 
C BEM VALUES 02/14/75 FOR NU2 DAY DRPIVER,1719 


DATA (SNOIQN(T),Tete33) 4 2,506 1078, 306 409,41 4OFe09,1 UNE S09, DRIVER, 1720 


* 1 FOF O92 GHEE 409, 2 SOF HI9“, 1 28E S09, ORIVER, 1721 
® VLOF HOA, TALDF SOT TARDE 06,2, SNE HVOF NDE FNS, 2, O0EF0541,90E 49S, VDRIVER, 1722 
SMF ENG, 2S, UNE FUG, 1 P0EHN4, 6 UE S03, 3, HOF FOS, 1 BOE SOS, 1 TOEHO3, OURIVER, 1723 
* He TOF HOA, GSE HOA, A LANE HNC, 1 IME SNA, 1 UOE SIAL 1 LSE S029 SOF HN, ORIVER, 1724 
# ROOFED ATONE SOL DAL OOF ODL, U,OOE SO! DRIVER,1725 
6 BFM VALIIES 02/14/75 FOR NU@ NIGHT URTVER, 1726 


DATA (SNOQNCT),TE1653) 4 SSOP IDOL 20E 10,7, TOF 49,2, UNE 409, DRIVER,1727 


* 2 S0F 60904, 156 409,4,55F 499,35, 006 409, ORIVER,S728 

# 1eA0F 409, 9, 20F 08,5, A0E S96, 3 ,00FH0R,1  UOF HOB, S SOEFO7T, 1 ,20F 407, VURIVER,1729 

# 1eDOF HNO, SeDTE FDU GL POF HOU, HUME HVS, T,40F 05,1 BIE F038, 1,10F #03, URIVER, 1730 

BALTOF ENA, Uy SOF HOP, A ,AVF ON, 1 INE SUS { GIF HN2, 1 e1SEFUP,9,50F +01, ORIVER,I731 

 ALN0E 01, 7,09F 01,6, 90F 491, 4,60E S01 DRIVER,1732 

ccc ORIVER,1733 
C # @ & ARTTHMETIC STATEMENT FUNCTION USED TO CALCULATE NTTRIC OXIDE Tv DRIVER, 1734 
C # «© # DAYTIME FOR ALTITIIDES BELOW 120, KM, ORIVER,I735 
cec OPIVER, 1736 
ANOMAF( BN Y s FRP(COEPCOCOC ECE AACLT3) BQ © AAC{2))#RY © AAC11)) *BOORIVER, 1737 

* *¢ AACIO)#RQ © AA(D)IPRAY © AACR) IRR © AAC 7) )*#AO DRIVER, 738 

* + bAf6))#BQ * AACS) ) #HD © BACK) ®BO © AACS)) #RO ORIVER,1739 

* * AAf2))#RO * Aaly)) ORIVER, 1740 

(of af 8 ORT VER, 1744 
C # # & ARTTHMETIC STATEMENT FUNCTITIN USED TY CALCULATE ATUMIC NITROGEN DRIVER,1742 
C # #@ @ IN MAYTIMFE FOR ALTITUDES ABNVF 40, AND RELUW 259, KM, ORIVER, 1745 
Ce¢ DRIVER,174u 
ANDAF( RQ ) zs FxP(C(C(OCO RRO) eA & BAC A)IRRY + BR(7))#RQ DRIVER,1745 

* + BR(B))*®RQ ¢ KRAOS)I#RO © BRC4))*RO # BAC) eR ORIVER 1746 

* * BR(2))#AQ # BR(1)) ORIVER,1747 

cce ORIVER, 1748 
C # & @ ARITHMETIC STATEMENT FUNCTION USED TO CALCULATE ATOMIC NITROGEN DPI VER, 17409 
C # # # AT NIGHTTIME FUR ALTITUDES FROM 45, KM TU 230, KM, ORIVER, 1250 
cee ORIVER,1751 
ANNAF( AR) ) s EPEC OO(C COC) #RN DRIVER,1752 

* + (C6) eBO & COS) #69 + COCO) HD * COOS)) BE DRIVER ,1753 

* + (C(2))eBO © £004)) DRIVER,17Su 

été ORIVER,1755 
Co @ & @ ACTTHMETIC STATEMENT FUNCTION USED TH CALCULATE W201 NEL TA) ORIVER,1756 
Cs @ @ IN MAVTIME FOR ALTITUDES RELD® 50, KM, ORIVER, 1757 
ccc ORIVER, 1758 
ANPSDF( BO Y & FExP(CO CEC OC OCE DOC11)#BR *& NOCH) eR * DA(9)) RG DPIVER,1759 

* # NDA) #RD # NO07I) 8B + NUL6)) BBA # HD(S))#BO ORIVER, 1760 


toa 


413 
yy4 
BIA 


en 
Bet 


eCanan 


c 


o 


10 


195 


110 


, 
* 
* 


* MOCG))#BO * NNC3)) #BQ © NVUC2))#KN © ADTY)) 


* ARITHMETIC 
* OAYTIME OR 


STATEMENT FUNCTION USED TL) CALCULATE WATFR FOR 
NIGHTTIME FUR ALTITUDES RELUW 129, KM 

AHKONFF( HQ) Ss 
* GG6010)Y#RQ * GGOO9)I*RY © GGCA))&AU & 66(7))#RG 
* GG06))#HO * GG6I5)) 459 + GEC4))#BY + GEOS)) #AO 
+ GG(2)) KO + 6G6(1)) 

® ART THMETIC 
* DAYTIME aT 


STATEMENT FUNCTION USED TO CALCUL ATE 
ALTITUDES RELUW 169, KM, 


NN2 FOR 


AN(IOFF( BN ) = FXPC ECO COC OC ECO HHOLS) eA 
+ HHOLO) IRAQ © HH(D) YR © HRCA) RO ¢ HH( 7) #RQ 
+ HH(6))#BO &¢ HH(S5)) ®BD © HH(G))*BQ + HK($)) #80 
¢ HHO 2)) #8 # HHOY)) 

6" TJ (100,200), KK 


INITIALTZATION, rts 
INITIALTZATTION, 

CONTINUE 

ALOGTE 3 ALNGIO‘ ExP(1,0) ) 

ATUMIC NITROGEN PRUFILE PARAMFTERS, 

IDURN ) 105,110,110 

NIGHTTIME N 

& AL YKM(1A) 

=z ALTKM(4G7) 

CALL FITTERCNALTNNGALTKM(1BY,ANNITE (IR) »NDEGNNe 1 6 2 6 

ANN71 3 ANNAFC JFYNN ) 

ANN72 & ANNAF( 72NN ) 

GN TO 115 
NAYTIME WN 

Z1ND 3 ALTKM(9) 

ZANN Bw ALTKM(47) 

CALL FITTERCNALTND, ALTKM(9) ,ANDAY(C9) »NDEGND, $ 6 2 » 

ANU71 = ANDAFC 7IND ) 

ANDI2 = ANDAFC FAEND ) 

CONTINUE 
NITRIC OxIpF PROFILE PARAMETERS, 
FIR DAYTIME NO, 
IFC7H,LT,77NU), WHERE Z7IND & 120 KM, 
WHERE THE POLYNOMTAL CUFFFICIENTS AACT) IN THE ARTTHMETIC 
FUNCTYON ANDDAPCZH) ARE DETERMINED RY SUBRULITINE FITTER, 

CALL FITTERCNALTNIDZ,ALTKM,ANIDAYSNDEGNO, 1 6 2 » AA) 
SFT ALTTT!IDE VARIABLES aT 115, 120, AND 125 KM, 

ZANT B@ ALTKM( 24) 

ZINO #£ ALTKM(25) 

ZANT) = ALTKM( 26) 
COMPUTE FITeFRIINCTTUN VALUES UF NO AT ALTITUNES 
AND ZINOIs12an KM, 

ANUZ26 & ANLIMAP (7HNU) 

ANO?7 &@ ANONDAFCT77NG) 
APPROYIMATE VERIVATIVE OF ALOG(NU) AT ALTITUDE 77NO=120 KM, 
DINOTZ, BY HISTNG THE FITeFUNCTION VALUES AT ALTITIIDFS Z6ND8 
115 KM AND Z7NU8120 KM, 


CALLED FROM SUIBRUUTTNE ATMUSU DURING 


Tee 


Z1NN 
Z2NN 
CC) 


8b) 


7oNOQ81LI5 KM 


119 


oT NT See wey 


NU = SNI(B) 2 ANODAF (ZH) 


DRIVER,17A4 
DRIVER, 1762 
DRIVER, 1763 
DRIVER, 1764 
URIVER,1765 


FKPC CC OCOC OC ECO GGE1L3) 8H & GECL2))#AQ © GOC11)) #BQURIVER, 1766 


ORI VER,1767 
ORJVER, 1768 
ORIVER,1769 
DRIVER ,1770 
DRIVER,1771 
DRIVER,1772 
DRIVER,1773 


@ HHO LA) RAY + RECT LI) ei DRIVER,I774 


VRIVER,I775 
URIVER,1776 
ORI VER,1777 
DRIVER,17786 
DRIVER,1779 
DRIVER,1780 
DRIVER ,17A1 
ORT VER, 17R2 
DRIVER,17AS 
DRIVER, 17AG 
DRIVER,17A5 
DRIVER, 176 
DRIVER, 17A7 
ORIVER,17A88 
URIVER,17R9 
DRIVER,1790 
DRIVER,1794 
DRIVER ,1792 
ORIVER,1793 
ORIVER,1794 
ORIVER,1795 
ORIVER,179%6 
ORIVER,1797 
DRIVER, 1798 
DRIVER ,1799 
DRIVER, 1890 
ORIVER, 1801 
ORIVER, 1802 
DRIVER, 1803 
ORIVER, 1804 
DRIVER,18605 
ORIVER, 1806 
DRIVER,1807 
DRIVER, 1808 
URLVER, 1809 
URIVER, 1810 
ORIVER, 1611 
ORTVER, 1812 
URIVER, S843 
ORIVER, 1814 
DRIVER, 1615 
DRIVER, 1816 


326 


336 


437 


a5 


35? 
353 
354 
356 
457 
3o1 
B70 
371 
Sie 
374 
376 
406 
407 
ayt 
ayy 
44 
Nes 
ueu 
427 
430 
451 
435 
4uy 
Gur 
au? 
ast 


AMAAMNAOANCANHAAIAIANQAAHANAANAAAN 


OLNOIZ & ALONG CANO7ZTSANUZOISOZINO @ZONO) 
SET THE CIIRF OF THE 10,7"CM SOLAR FLUX SRAR, SBARY, 
SRAR$ = SRARawS 


PRESCRIBE THE DERTVATIVE OF ALUG(NG) AT ALTITUDE 7BNUS{25 KM, 


DLNO87, TO RE GIVEN RY 
DLNOB? = 0, 06573*SBART/(SHARSZ45,ES) @ O,0756U @ 0,027N3HDLNOTZ 
FOR LATER USE, COMPUTE 
ANQZ78 3 ANU7THEXPC2A],5#(NLNITZ 4)! NIVB7)) 
IFCZH.GE,Z77NU QAND, ZHALT, ZAIN I) « THEN 
SNT(R)BANIZ7THEXPOZHMIOOROLNUTZ + OL *ZHML 20K eee(DLNOTZ©& 
DENODBZ)) 
WHERE ZHMIONs7He1>A0, 
IF(7H,GF,246NU), THEN 
SNICAYSANO76REXPOZHM OS*#NLNUAZ) 


WHERE ZWM12Ss7HeLAS AND AN(IIZBSANUZTRE KK (2, S#(DLNUTZ+DL NI18B7)), 
AY NIGHTTIME, NO NIFFERS FROM DAYTIME NO RELUW ALTITUDE ZUNNTN 


ZAS KM AND ABNVE ALTITUNE ZSNONZINO KM, 


IFC 7H,LY,Z1NUN), WHERE Z1NONSS5O0 KM, NOZANONZYSANUNITE CL IIS Ue 


IF (ZH.GE,Z1NON @AND, ZH,LT,Z2NUN), WHERE Z2NUNESS KM, 
NUSANONZA MF XP( (ZHe ZENON) KANUNST) 
WHERE 
ANONZ2SANUNIT ($2) 51,64 
ANONSTBALOGCANOANZ2/ANUNZ1)/0ZANDN@ ZENON) 
IF C2HeGE,ZQNUN gAND, ZHYLT 2UNUN), WHERE 7UNUNZSRS KM, 
SVT (RIEL O,weSUM( X(T) a ZM2NUNMR (GWT) ), 151,49 
WHERE ZMONNS7H@ZANON AND THE NINE COFFFICIENTS x19 ARE 
DETERMINED §Q THAT NUT ANLY ALUGIO(NO) EQUALS THE NIGHTTIME 
VALUES FUR NO AT 7H355(5)85 KM BUT ALSU THE SLUPE UF 
ALOGIO(NO) TS CONTINUOUS AT SS AND AS KM, 


THE NIGHTTIME CONSTANTS FOR ALTITUDES BELOW B85 KM ARE NOW SET, 


ZINTNBALTKM(T1EY 
ZANONZALTKM( 12) 

ANUNZ¢SANONT T0144) 

ZANON = BS, 

ANUNZ2 & ANONTT (42) 

DN 120 137,A 

ZINONCT) & SL TKM(TOI0) 
ANUNZT(T) & ANONT?( T4140) 
CONTINUE 

ANONZT(A) = ANOMNAF( ZINONCG A) ) 
ANUNST © ALNGCANONZI (2) /ANUNZ1)/(°Z21NUN(2) @ZINONY 
DINM2Z2 = ANONST#ALOGTF 

x (KB) gs DLNIe7 

x09) 8 ALOGIOCANONZ2)Y 

DM 125 33,8 

ZIMANNCT) & ZTNONOCTI@ZINUN(G 2) 
CONTINUE 

DN 130 Ta1,6 

27112 = ZIM2NUC I +2) 

ACTA?) & ZLTHOKRTTT2 

DM $30 J81,6 

ACI,7eJ) = 7ITAQRACT Bed) 
CONTINUE 

ZT1A s&s 7IMENUCA) 

AC7T,7) & 2,e7TTA 

ON 138 Jat,6 


URIVER,1817 
DPIVER, 1818 
OPIVER, 1619 
ORIVER, 1820 
ORIVER, 181 
DRIVER, 1622 
DRIVER, 1823 
DRIVER, {B24 
URIVER 1625 
ORIVER,1626 
DRIVER, 1627 
DRIVER, 1828 
ORIVER.1629 
VRIVER, 1830 
URIVER 1831 
DRIVER, 1432 
DRIVER, 1833 
ORITVER, 1634 
DRIVER, 1835 
DRIVER, 1836 
VRIVER, 1837 
DRIVER, 1834 
DRIVE®,1839 
DRIVER, 1840 
OPIVER, 18443 
ORI VER,184e 
DRIVER, 1643 
DRIVER, 1644 
DRIVER, 1845S 
URIVER, 1646 
DRIVER, 1647 
DRIVER, 1A4B 
DRIVER, 1649 
DRIVER,18650 
ORIVER, 1851 
URIVER,18S2 
ORIVER, 1853 
URIVER,1854 
DRIVER, 1885 
URIVER,1856 
DRIVER, 1857 
ORIVER, 1858 
VRIVER,1859 
DRIVER, 1660 
ORIVER, 1661 
DRIVER, 18b2 
DRIVER, {ANS 
DRIVER, 1664 
DRIVER,.I R65 
DRIVER 41866 
ORIVER, 1867 
ORIVER, 1668 
URIVER, 1869 
ORIVER, 1870 
ORIVER,VAT4 
URIVER, 1872 


tt ee Ml 


MTN 
457 
uot 
473 
O74 
a7? 
507 
512 
Ss 


Su7 
559 


553 
S54 
SSA 
Sao 
So? 
S72 
S738 
S75 
S76 
500 
607 
a1? 
a1A 
&e0 
hel 
624 
A34 
Ay? 
641 
su? 
84S 
45? 
653 
665 
467 
671 
ran 
701 
13 
744 
717 
721 
12? 
Yeu 
725 
To? 
74h 
Tus 


NoOaAND 


Lae 


14a 


146 


148 


15 


~aeere 


FJ1 = Jel 

ACT, 7TeJ) 2 ZITARCC RIL OL, SE IL) ACT, Be S) 

CONTINUE 

ON 140 Tst,6 

A(T,3) = ALOGIOCANONZT(IT4¢2)) © X(R)*LZTMANN(T #2) © KO9) 
CONTINUF 

ZANON s&s ZINTINCR) 

a(7,8) ALNGTE ® (OCOCOCOCE C12, *#AACL3) #ZBNON 


$1,"4A012)) @27BNON # 10,8AA011))#2BNUN © 9,4A4(10)) *ZANDN 


= 
+ 
+ Ry epAACISIATBNUN © TLRAACA)I®ZFBNON + &,8AA(7))®7BNIIN 
+ 5,*44(6))*#ZBNON # G,*AS(S))®Z7BNON # 3, #AA(4))®&Z7ANON 
+ >, eAACS)IRZBNON AAC CY) @ KEB) 
NO 3 7 
CALL Sit VEC A,X,NOD) 
IFC 74,.GF,725N0), wHERE Z5NOZ4LTKM(24)5100 KM, 
ANONTTSANIIDAY RCC ZH*1N0 eG) *#2 © T2N0GIS(10 ge ZHe100, #2 
+ 7200,)) 
MOLECULAR (KYGEN (SINGLET DELTA G) PROFILE PARAMETERS, 
ZN2090 = ALTKM(19) 
ZM210N0 = ALTKM( 271) 
anengo = N28dGN(19) 
BN2N9N & eALOG( H2SNGN(22)/402090 IACALTKM(22)"202090) 
IFC IPURN > 142,150,159 
ZN2070 & AL TKM(15) 
ZN2080 fF ALTKM(ITS 
AN2970 = F2SNGN(15) 
4N2080 = N2SNEN({7) 
BN2070 = eALOG( AN208H/402070 )/4(702080"202070) 
7fo) = ALOGIOC AU2080 ) 
NO 164 [31,4 
ZTI12 = ALTKM(TO17I@7OA040 
4(I,5) = 7IT2 
ON 144 Jut,4 
ACI,SeJ) = 2ZIT2eAlTs60J) 
CONTINUE 
ZT1AR = 2921009792980 
At5,5) 2 1,9 
A(5,6) = eBN2NGHALOGTE 
09 146 Js1,4 
a oe oe: 
ACSsSeJ) = 7ITAR( CR JL. )/F SIA (5,685) 
CONTINUE 
ON 114A Ts1,3 
AfI,6) 2 ALNGIOC N2SNOGN(1417) ) & 26) 
CONTINUE 
AfU,6) = ALNGIOF AD2090FEXPleRD2090R(ZNALN08ZNANIN)) 1 wm 766) 
NN 3 § 
CALL SUL VE(A,7,%1)) 
69 TO 3156 
ZN2N050 = ALTKM(14) 
292078 =f ALTKM( 1H) 
AN205N & N28NG011) 
AN20975 & N28SDGO(146) 
HM20SN = eALNGE ANZOTSZAI2U50 IACZTUZA 0TH HLNZ080) 
CALL FL TTERONAL TUD, AL TKY, D2ENGDNNGO2N, $ » 2 2009 
Y(o) s ALNGIO Anr07S ) 


121 


ORIVER, 1873 
DRIVER, 18674 
DRIVER, 1675 
ORIVER,1876 
DRIVER ,1877 
DRIVER, 1878 
ORIVER, 1879 
ORIVER,§#ARO 
ORIVER,S@AL 
DRIVER, 16R2 
DRIVER, 1683 
ORIVER,(8R4 
ORIVER, 1885 
ORTVER,16R6 
DRIVER, 1687 
ORIVER,{6AB8 
ORIVER,1BA9 
DRIVER, 1890 
ORI VER, 1891 
DRIVER, 1892 
VRIVER, 1893 
DRIVER, 1894 
ORIVER, 18695 
DRIVER, 1896 
ORT VER ,1897 
DRIVER, 1698 
ORIVER 1899 
ORIVER,1990 
ORIVER, 1901 
VETIVER, 1902 
ORIVER,1903 
ORIVER, 19904 
ORIVER,1905 
DRIVER, 1906 
DRIVER,I907 
URIVER 1908 
URIVER, 1909 
DRIVER, 1910 
DRIVER, 19114 
ORIVER,1982 
DRIVER, 1913 
ORIVER 1914 
ORTVER,I915 
VRIVER,I916 
DRIVER, 1917 
DRIVER, 1918 
URIVER,IN919 
DRIVER, 1920 
DRIVER, 1921 
DRIVER ,1922 
URIVER, 1923 
DRIVER, 1924 
DRIVER, 1925 
DRIVER, 1926 
ORIVER, 1927 
URIVER, 1928 


ee a ee ee ee Oe eee 


Tu? 
52 
757 
769 
Tol 
779 
T7646 
1004 
1907 
init 
1042 
1916 
1923 
1925 
1037 
1040 
1041 
1046 


1046 
1047 
1051 
1953 
1962 
1965S 
1071 
1972 
1074 
1101 
1107 
1115 
1120 
112! 
1123 
1125 
1133 
1134 
1146 
1159 
115! 
1153 
1156 
1157 
11o! 
Line 
1164 
1173 
1203 
1204 
174% 
re 
1215 
1724 
1232 
124% 
1243 


152 


158 


140 


lhe 


Lau 


¥OS) = @BN20SOeKALOGTE 
ON 442 131,3 
ZTI2 Bw ALTKM(Te16)0202075 
ACI,4) = Z1T2*7TI2 
A(I,5) 2 ALNGIO€ O28OK6DCI +16) ) = ZIL2ev¥(5) © lo) 
ON 152 J2t,3 
A(1l,GYeJ) & FITPeALIT 5 eJ) 
CONTINUE 
ZTIAR = 211209087292075 
AC4U,4) @ A,eIZTIA 
AfU,5) = @eBN2O0ONeALNGTE @ YS) 
09 154 Jsi1,3 
FY ws Jet 
ACU,UeJ) = 7PITRe((FJ41,)/F J) eA (4,50) 
CONTINUE 
NO 2 4 
CALL SULVF(4,Y,NU) 
CONTINUF 
OY CU7NNE) PROFILE PARAMETERS, 
ZN3NWO 3s ALTKME 9) 
ZN3N75 = ALTKM( 16) 
ANSNGO = O3NAY( 9) 
BN3NGH & eCALNG( USDAY(16)7ANZD40 IS(CLZNINMeZO4OU0) 
EF(10) s ALMGIOC NSMaAV( 1) ) 
O09 15A 131,48 
ZII2 = ALTKM(T41) 
A‘IT,9) &@ 7IT2 
ACI,19) = ALNGIN( USDAYCI*1) ) © FECL) 
DN 158 Jat,R 
ACT, 9eJ) 2 ZIT2#A(1,100J) 
CONTINUE 
ZTIA & 7113040 
4(9,9) = 1,9 
AC9,19) 3 eAOSP4O*ALUGTF 
AN 100 Js1,4 
PFlheilJ 
A(9,9eJ) B ZTTAR( CR IO1, SF JI eA(9, 1083) 
CONTINUE 
NO 2 9 
C4Ltl. SULVE(A,FE,NN) 
IFC TMURN ) 162,172,172 
ZN3NSS = ALTKM(12) 
ZOSN70 @ ALTKM(45) 
ZO3N75 & ALTKM(16) 
AN3N79 & NBNTT(15) 
MOBNTO B@ @ALUGE LIBNTITIL6)ZANINIO YACZO3NIS#@Z03N70) 
VO$fo) & 
VOCS) B& eRONBMGNMALNGTE 
ON {64 T381,3 
ZT12 =f ALTKM( T6912) © ZOINSS 
ACI,4) 2 2172¢7T12 


A‘I,5) & ALONGIO€ NSNIT(Te12) ) | ZITeevn3s(5) © VU3(6) 


ON 164 Je1,3 

ACT, UeJ) = PITA#ATT Sed) 
CONTINUE 

ZTIA & 2PIINTOSZIINSS 


122 


ALNGIO€ AUSOUOFE XP CeRI3D40*% (703N5590750400)) 


OR1VER, 1929 
URIVER,1930 
ORIVER,I9ON 
DRIVER, 1932 
ORIVER 1935 
OURIVER,1934 
ORIVER,1935 
ORIVER, 1936 
ORIVER,1937 
ORIVER,1936 
DRIVER,1939 
UVRIVER, 1940 
DRIVER, 1941 
ORIVER, 1942 
ORIVER, 1945 
ORIVER, 1944 
DRIVER,1945 
URIVER, 1946 
DRIVER, 1947 
DRIVER, 1948 
DRIVER, 1949 
ORITVER,1950 
DRIVER, 1951 
ORIVER,1952 
DRIVER 1953 
VETVER,1954 
DRIVER,1955 
DRIVER, 1956 
DRIVER 1957 
DRIVER, 1958 
ORIVER,1959 
URIVER,1960 
URIVER, 1961 
ORIVER, 1962 
ORIVER 1963 
DPIVER, 1964 
DRIVER, 1965 
ORIVER, 1966 
ORIVER, 1987 
ORT VER, 1968 
DRIVER,1949 
DRIVER, 1970 
ORIVER,I97]4 
URIVER, 1972 
ORIVER,I973 
ORIVER,I97TG 
DRIVERI9TS 
ORIVER,1976 
VRIVER,1977 
ORIVER, 1976 
ORIVER,1979 
ORT VER, 1980 
DRIVER, 1944 
ORIVER 1982 
VURIVER,19R3 
DRIVER, 1984 


| 
| 
| 
. 
| 


1245 
12aea 
1282 
1257 
1761 

1273 
1274 
1275 
1300 
iso! 

1303 
1343 
1315 
1320 
1325 
1326 
1327 
13356 
13a4 
1352 
1355 
1357 
1360 
13e4 
1371 

1373 
1408 
{406 
1407 
1412 
Y4y4 
1415 
1417 
1427 
1431 
1434 
laut 

1440? 
14u8 
1452 
1460 
1466 
1471 

1473 
1474 
1509 
1505 
1507 
152! 
1922 
1823 
1530 


1559 
1537 
1540 


lho 


1468 


170 


172 


174 


A(4,4) B A,ezTTA 

AGS? 8B PHIBNTORALNGTE = VNGI5) 

DY 166 Jat,3 

FJ @ Jol 

AfU,UeJ) B Z2ITHO((F Jos, I/FJ)eA(4,5eJ) 
CONTINUE 

NOX 3 4 

CALL SULVE(A,VOS,NO) 

ZOSNGO w AL TKM(19) 

ANSN90 @ NISNITC19) 


BIBNGON = eALOG( USNITC22)7AN3N90 IACALTKM( 22) @ZI5N90) 


wN$fo)d w@ ALNGIOC OSNIT(16) 9 
aNSOS) 3 eKIENTOMALOGTE 

ON 16A Tat1,3 

ZT1> = ALTKMETet16) © ZFOSNTS 
ACI,4) & 2ITe*e7112 


ACI,5) & ALNGIO€ OSNITCT 416) ) © ZIleewh3(5) © wUS(A) 


OF 168 Jst,3 

ACIl,4eJ) & 7ITA#ACT Sed) 
CONTINUE 

ZTIA & JIBN9ONEZNINTS 

a(4,4) 8 2,#7T7A 

A64,5) @ @ROSNGO®ALOGTE © w3(S) 
NO 170 Jei,3 

FY 3 Jol 

ACG,GeJ) @ ZITAA( (FJ et, )/F IV aAC4,SeS) 
CONTINUE 

Nn #8 U 

C4LL SOLVE CA, 493,N0) 

69 TI) 178 

ZN3N9O S&S ALTKM(19) 

AN3ND90 NgNAVYC19) 


= 
BO3N9G0 & wALOG( (BDAY C22)78N3ZN90 ISTALTKM( 22) €203090) 
<= 


UN3(6) ALNGIO N30av(to) ) 
UNSTS) & eRNBPUNeALOGTE 

0 174 T31,3 

Z7112 & ALIKM(Te16) © JFOSDIS 
afl,a) & 21Te#7T12 


ACIL,5) & ALNGIO’ A3SPav(Telod Y & 7ZIT2euns(S) © UWUSl6) 


PN 170 Je1,3 
A(],4eJ) & 7IT2eAll Se) 
CONTINUE 
ZTIA = 793090872N3N75 
A(4,4) ® 2,#711A 
AfU,5) & @eBOSNGNRKALOGTE © UNZ(S) 
ON 176 J81,% 
FY s Jot 
ACU,4eJ) & ZITAM(CF Joy sE JI HAC, Sed) 
CONTINUE 
vn 8 4 
CALL S0OLVE(A,UD43,NO) 
CONTINUE 
FIT COEFFICTENTS FUR NUD (DAY AND NIGHT), 
CALL FITTER (NKMNU2,ALTKM,SNN2D,NDGNO2, 1 » 2 HH) 
HNO21D & ALTKM(29) 
HNUP2N mw ALTKM(35) 


123 


DRIVER IIRS 
DRIVER, 19R® 
DPIVER,IOR? 
URIVER,19R6 
URIVER,{9R9 
UR] VER, 1990 
ORIVER, 1994 
DRIVER,1992 
DPIVER, 1993 
URIVER,1994 
DRIVER,1995 
UHL VER, 199% 
ORIVER,1997 
DRIVER, 199K 
VRIVER, 1999 
OP VER, 2090 
ORTVER, 2001 
DRIVER, 2002 
DRIVER, 2A0NS 
DRIVER, 2004 
DRIVER, 2095 
ORI VER, 20%6 
DRIVER, 2097 
ORI VER, 2008 
ORIVER, 2009 
ORI VER, 2010 
DRIVER, 2011 
DRIVER, 2012 
ORIVER 2015 
DRIVER, 2014 
ORIVER, 2015 
URIVER, 2016 
ORF IVER, 2017 
DRIVER, 2018 
ORIVER,2019 
ORIVER, 2020 
ORIVER, 202} 
ORIVER, 2022 
ORIVER, 2023 
VDRIVER, 2024 
VRIVER, 2025 
URIVER, 2026 
ORIVER,2027 
DRIVER, 2028 
DRIVER, 2029 
ORI VER, 2030 
ORI VER 2031 
ORI VER, 2032 
ORI VER, 2033 
URIVER 2044 
ORIVER,2035 
VRIVER, 2036 
ORI VER, 2037 
DRIVER, 2038 
DRIVER, 2039 
DRIVER, 2040 


1542 
{S45 
1S47 
1552 
1563 
1564 
1566 
1S67 
1871 
1573 
1576 
1977 
1400 


1603 
1631 
1615 
1615 


1645 
1421 


14625 
1630 
14633 
1633 
1636 
14641 
1641 
1443 


1643 
1646 
14651 
165! 
1654 
1657 
1657 
1661 


1461 
1464 
1467 
167! 
1476 
1700 
1702 


20 


oO 
co) 
ion] 


oO 


ao 


ANANANANAAA 


ANUDFD 
HNUODD 
RNUOA12 
PNO255 
AN()265 
HNO255 
HNO? 65 
HNUADN 
RNU2PFA 
AND 2K? 
HNO2a2 
HNOPDA 
RNII2P RO 


AND QFE C HNDODN ) 

HNO AI NeHNU22N 

ANN OFFC HNO 21D ) 4 ANDAFN 
ANN 2FF( SS, ) * ANNDAF( SS, ) & ANONZ2 
SNOAN( 14) 

ALTK™ (12) 

ALTKM( 14) 

HNO ASS eHNO 265 

FNO 25S /AN0 265 

ANNOFFO Re, ) 

Re, 

HNO QAS @HNUAR? 

AND 265/A4NO02B2 


FIT COEFFICTENTS FOR Han, 
CALL FITTERCNKMH2N,ALTKM,H20DN,NDGH20, 1 » 2 ,GG) 
HPU120 = AH2ZUFFC 120, ) 


RE TIAN 


CONTINUE 


AN ERRONEQUS CUNDTTION WILL OCCUR IF SPCMIN IS CALLFO wITH 
KKs? AND A GIVEN VALUE NF ZH TF ATMOS!) HAS NOT BEFN CALLED 
FIRST WITH KkK=2 AND FOR THE SAME VALUE OF 2H, 

THE VARTARLE ZHFLAG IS HISEOD 19 DETECT THIS CONDITION AND 
TO MAKE THE REGUTRED CALL 19 ATMUSU, 

ZHFLAG TS INITIALTZFED TO AN ARBTTRARY NEGATIVE VALUF IN 
THE INIVIALTZATION CALL TO ATMUSU, 


TFC ZANE, ZHFLAG Y COLL ATMNSUC2,2Z+4) 
COMPUTE DENSITY (JF N 

IFC IMORN ) 240,215,215 
NIGHTTIME N 

TF( ZHe2ONN ) 212,211,711 

SNICT) = CANN72/S1Z2ON)*SNI(1) 


GN TH 220 


IF( ZRe7INN ) 214,213,213 
S107) 2 ANNAF( 2H ) 


GN TO 220 


SNI(7) @ ANNZ1 


GA T1) 220 


DAYTIME N 
TFC ZHeZ2Ni) ) 217,216,216 
SNIC7) = CAND72/81Z2Z2N) *#8NI (1) 


GN TO 220 


IF( ZHeZIND ) 719,218,218 


SN1(7) 


GNA TO 220 


AND AFC ZH ) 


SNI¢7) w= ANDZ1 

CONTINUF 
COMPUTE DENSITY OF NO 

IF(7HeZ7INI}) 2?he2Pir2rt 

IF C7HmZANN) 223,27?21272 

ZHM78 & ZHe7BN() 

SNICBY & ANDZAREXPCZHMZARDOLNDAZ) 


GN Th 


274 


ZHM77 & ZH @ 7INQ 
SNICR) B ANTZ7THEXP(ZHMZ7HOLNOITZ©0 Le ZHM IT eee (DLNOTZ@OLNUBZ) ) 


ORIVER, 2044 
ORIVER A042 
DRIVER, 2043 
URIVER,2Qgu4 
DRIVER, 2045 
ORIVER, 2046 
ORIVER, 2047 
ORIVER, 2048 
VRIVER, 2049 
URI VER 2050 
ORTVER, 2051 
ORIVER, 2052 
VRIVER, 2083 
ORIVER, 2054 
ORIVER, 2055 
DRIVER,2056 
ORIVER,2057 
ORIVER, 2058 
ORIVER, 2059 
ORIVER, 2060 
ORIVER, 2061 
ORIVER,2062 
DRIVER, 2063 
ORIVER 2064 
DRIVER ,2065 
DRIVER ,20h6 
DRIVER, 2067 
ORIVER, 2068 
ORIVER 2069 
DRIVER, 2070 
DRIVER, P0711 
ORT VER, 2072 
DRIVER 2073 
ORTVER 2074 
ORIVER, 2075 
DRIVER, 20% 
ORIVER,2077 
DRIVER,2078 
ORIVER, 2079 
ORIVER. 2080 
ORIVER, 2081 
DRIVER, 20R2 
URIVER, PORE 
URIVER,20R4 
DRIVER, 20R5 
ORIVER,P0R6 
ORIVER,20A7 
DRIVER, 2088 
ORIVER,20R9 
DRIVER, 2090 
ORIVER, 2091 
DRIVER, 2092 
DRIVER, 2093 
VRIVER,2094 
ORIVER, 2095 
ORI VER, 2096 


224 1F( IMURN,GF,O ) GO 


™ 229 


225 ZHM792 = (ZHel0n,)#*2 


SNTCBY Ss SNTCR)*(CZHMZ52¢7209, D/C 1% e*ZHM7ZS2¢7200,) 


GN TU 229 


276 IF( IMORN,GE,O ,OR, ZH,GE,ZUNON ) GO TU 2783 


IF( ZHeZ2NON ) 228,227,227 


227 ZM2NUN & ZHeZANIN 


SNICBY = LOgeP COC OCOCEXCLIRIMANIN & X02))H@ZM2NUN @ ¥ (3) 274 2NON 


* # X€4))*#ZM2NUN @ KX(5))"ZMPNDON © K(6)) @ZM2NUN 
* + X(€7))®ZM2NON # X(8))®ZMPNON ¢ X(9)) 
GN TH 229 


228 IF( ZHeZ1INON ) 22R2,2781,22A1 


22R4 ZM2NON & 7HeZANDIN 


SNICR) 2 ANONZ2#EXPCZM2NON®ANONSI) 


69 TO 229 
P2R2 SNI(B8) = ANONZ] 
GN TO 229 


22R3 SNIC(B) 2 ANNNAFC ZH ) 


TEC IMORN,GF,9 4%, 


GN TQ 225 
229 CONTINUE 


ZH LT,i190, 


» 60 TO 229 


COMPUTE DENSITY OF O21 DELTA G) 


IFC ZH OLE. 202100 


) 


GN TO 231 


230 SNIC13) =F ANZNQH KEXP (eBN2090e (ZHe2U2090)) 


GN TO 238 


esi TFC INMORN ) 232,235,235 


NIGHTTIME 112 ( 
232 IFC Z4,6T,.292070 
SNI€13) = anen70 
GN TT) 238 
2335 IFC ZH,GT, ZN 2040 


{ 
) 


) 


DELTA G) 
GN TO 233 


GN To e234 


SNIC13) & ANZOTH aE KE CeRN ZIT Hal ZHeZU2A070)) 


GN TO 223A 
234 ZHMKM gs JHeT(2580 


SNICH3) FE LO, eRe CEC OZEL IV @ZHMKM © 702) eZHMKM & 203) I E7HMKM 
@ 2(4))#ZHMKM © 7(5) 28 7HMKM &© 7(6)) 


GN TU 238 
D4YTIME 071 
23s TFC ZH,GE.292090 
IFC ZH,GE,.Z2N2959 


NELTA G) 


) 
) 


69 TO 230 
GDI TO 236 


SNI(13) = AN2SDF( ZH ) 


GN TO 23a 
236 IFC 24,61, 2712075 


) 


GN YO 237 


SNIC13) & ANZN5O#F XP (eB0205N#(Z7HO702050)) 


GN TO 238 
247 ZHMKM = ZHeZ02075 


SNIC13) = 1O,*e COCO COV OL @ZHMKM © V(2))R7HMKM & ¥(3) DR 7HMKM 
# ¥(U))@ZHMKM @ ¥(5))&Z7HMKM @ ¥(6)) 


* 
23a CONTINUF 


COMPUTE DENSITY 


IF¢ ZHOLT. 27 gAhaN 


b) 


OF O§ (NZONE) 


GN TO 243 


IFC INURN 9 239,244,744 


NIGHTTIME 3 

239 TFC ZH ALT .ZN3N5S 

IFC ZH,GE.2935N70 
ZHMKM & ZHe7N3NSS 


b) 
) 


GN TO ead 
69 TO 240 


DRIVER, 2097 
DRIVER, 2098 
ORIVER, 2099 
ORIVER, 2100 
DRIVER, 2101 
DRIVER, 2192 
ORIVER, 2103 
ORIVER, 2104 
DRIVER, 2195 
DRIVER, 2106 
DRIVER, 2107 
ORIVER, 2108 
DRIVER. 2109 
DRIVER, 2110 
ORIVER, 2131 
DRIVER, A t{2 
ORIVER, 2113 
ORIVER, 2114 
DRIVER, 2115 
ORIVER, 2116 
DRIVER, 2157 
DRIVER, 2118 
DRIVER, 2119 
ORTVER, 2120 
ORIVER, 211 
ORIVER, 2122 
DRIVER, 2123 
DRIVER, 21424 
ORIVER, 2125 
DRIVER, 2126 
ORIVER, 2427 
ORIVER, A128 
DRIVER,2129 
DRIVER, 2130 
ORIVER, 2131 
DRIVER, 2132 
ORIVER, 2133 
URIVER, 2134 
DRIVER, 7135 
ORIVER, 2336 
VRIVER, 2337 
ORIVER, 2338 
DRIVER, 2139 
ORIVER, 2440 
ORIVER, 2444 
ORIVER, 2342 
DRIVER, 2143 
ORIVER Ay 4u 
DRIVER Agus 
DRIVER, 2146 
DRIVER,2147 
ORIVER,2148 
ORIVER, 2149 
DRIVER, 2150 
ORIVER, 2151 
ORTVER, 2152 


eas 


elo! 
21ne 
elob 
217s 
e176 
ePot 
2203 


eri? 
2220 
2227 


esi 


2253 


2254 
2760 
2267 
2279 
2°73 
2275 


astt 
23ie 
2322 


2327 


eset 
2%27 
233s 
2436 
2esu? 
23so 
2tsn 


etn? 
2362 
etre 


2370 
2477 


237? 
2402 
2404 
24yu 
244s 
244% 
euy4 


2ue 


ens 


250 


esi 


252 


ess 
24} 


262 
eas 
2o9 


VO3(3) de ZHMKM 
VO3(4)) 


SNIC1G) = IO ee CC OCOVOSCL I ®ZHMKM & VO3C2) IE ZHMKM 6 
# VORC4))@ZHMKM © VUS(5))@ZHMKM ¢ 

G9 TU 247 

IFC ZH,GT,293N75 3 GN TO 244 

SNIC1G) & ADENINE KP (@BOIN7TO RC ZHe7USNTO)) 

6 Tu 247 

IFC 24,GE.293N90 ) 

ZHMKM & ZHeZUBNTS 

SNICLG) = 10, ee COC CWNSCL) ®ZHMKM © WUS(2)IRZHMKM + 
@ WOBCE))eZHMKM + WO3(S)) MZHMKM 6 


GN TO 242 


WO3C3) )@ZHMKM 
wO3(6)) 
GO TOU 247 
SNIC14) & ANBNQNRE XP (eBNINGHN& (THe ZUINIO) ) 
6M TO 2487 
IF ZHALT,404 6 BITH DAY AND NIGHT USE FOLLOWING, 
SNICI4) = 1D ROCCO OC CCCEE CLI &ZH © EEC2)) #2ZH © EEC 3)) ®ZH 
+ FECU)) ZH © EECS))eZ2H + FE C6)) #ZH © EEC(7)) a2 
¢ FEC R)I*ZH + FECQOY)#7H © FE(30)) 
69 TO 247 
DAYTIME O8 
IF( ZH,GT,293N75 1 GN TH 245 
SNICL4) & ANZDUNKE XP (eBNZDUOe(ZHe7USD40)) 
GV 19 247 
IF( ZH,GE,293N90 ) 
ZHMKM @ ZHe7N3075 
SNIC1G) & 1M peal COC CUNS (1) ®7HMKM @ UD3(2))®ZHMKM & IUS(3)) RZHMKM 


69 TU 246 


+ UORCE)) ®7HMKM @ UOS(S))8ZHMKM # UB(6)) 
GN Ti) 207 
SNIC1{G) 2 ANBNGHREXP(eBO DIG ew ( ZHe7ZU3090)) 
CONTINUE 


COMPUTE DENSITY OF NO2 
IFC TPORN ) 248,252,252 

NIGHTTIME NNO 
IF( ZH,GE,HND255 ) GN TO 250 
SNIC1S) = ANOZFF( ZH ) ¢ ANTDAFC ZH ) @ SNI C8) 
GO TO 261 
IFC ZH,GT,HNO26S ) 
SNIC1S) s ANO?65 * 
GN TO Pbl 
IFC ZH,GT,HNI)2R2 ) 
SNI15) = ANDPRA *# 
69 TO 261 

DAYTIME NN2 
IFC ZH,GT,HNO22N Y GY TO 253 
SNIC1S) & ANDCFF( ZH ) 
6 TO 261 
SNIC15) 2 AN()PFN & RNN212*8( (7HEHNO229) /HNO2DD) 
CONTINUE 

COMPUTE DFNSITY OF Hen 
IF( 24,GE,120. ) GO TN 262 
SNI(16) 3 AHQNFFC ZH ) 
60 TO 263 
SNIC16&) 3 HADIA0RE KP (e0, 1668 (7HO120,)) 
CONTINUE 
RETURN 
END 


69 TO 251 
RNO2F ate ( (Z7HEHNUZOS) J/HNUCON) 


GN TO 252 
RNN QR Pea (THEHNUAG?) SHNUPDB) 


(DAY AR NIGHT) 


DRIVER, 2353 
ORIVER, 2154 
ORIVER, 2155 
ORIVER, 2156 
ORIVER, 21587 
ORI VER, 2158 
ORIVER, 2159 
DRIVER, 21460 
ORIVER, 2141 
URIVER, 21462 
ORIVER, 2163 
DRIVER, A, 6u 
ORIVER, 2465 
ORIVER, 2166 
ORIVER, 21467 
ORIVER, 2168 
DRIVER, 2169 
DRIVER, 2170 
DRIVER, 2174 
ORIVER, 2172 
ORIVER 2173 
ORIVER, 2174 
DRIVER, 2175 
ORIVER, 2176 
DRIVER, 2177 
DRIVER, 2178 
ORIVER, 2179 
DRIVER, 2180 
ORIVER, P11 
DRIVER, 218e 
DRIVER, 2yR3 
URIVER, 23 Ru 
DRIVER, 2485 
DRIVER, 2186 
ORIVER, 2187 
ORIVER, 2168 
DRIVER, 2389 
DRIVER, 2190 
ORIVER, 2191 
ORI VER, 2192 
VRIVER,2193 
ORIVER, 2194 
DRIVER, 2195 
DRIVER, 2196 
ORIVER, 2197 
DRIVER, 2198 
DRIVER, 2199 
DRIVER, 2200 
DRIVER, 2201 
DRIVER, 2292 
DRIVER, 2203 
ORIVER, 2204 
DRIVER, 2295 
DRIVER, 22% 
ORIVER,2207 


45 


ZTTOUT 


S'IBRONTINE FYITOUT 


CCC 

Cc SUBRQUTTINE 7TTOUT CONVERTS & GREGORTAN CALENDAR DATE (20 TH 
(2 CENTURY yEaR Ty¥RS, MONTH TMONS, Lay IDAYS) AND ZONE TYMF IT 
Cc AT FAST LONGTTUNE PLON TO GREGURIAN CALENRAR NATE AND MEAN 
(s TIME UT aT GREENWICH, 

ccc 


REVISION 92 (11/41A/74) PROVIDES ee, 
1, TEST FOR LEGAL INPUT DATE, 
INPIIT PARAMETERS 
TY8S =» NUMBER OF THE YEAR IN THE 1909 S CE,Gye 1974 
BFCOMFS 74), 
IMONS @ NIIMBER OF THE MONTH (CF,G,, FERRUARY BECOMES 2), 
INAYS e DAY OF TRE MONTH 
71 © ZONE TIME FOR THE ISeNMEGREE LONGITUDE INTERVAL 
CONTATNING PLON (NECIMAL HRS) 
PLON « EAST LUNGITUDE OF POINT P (RANTANS) 


OUTPUT PARAMETERS 
TyY89s 6 A PYSSIBLY REVISED VALUE OF THE INPUT PARAMFTFR, 
CORRESPONDING TO GREENWICH, 
IMONS @ A POSSIALY REVISED VALUF OF THE INPUT PARAMETER, 
CORRESPONDING 19 GREENWICH, 
IMAYS @ A POSSIRLY REVISEN VALUF OF THE INPUT PARAMETER, 
CORRESPONDING TN GREENWICH, 
UT © UNIVERSAL TIME (DFCIMAL HRS) 


DEFINTTTON OF HATA 
INMaYMON(T) s&s DAYS TN THE IT TH MONTH OF & NON@LFAP YEAR 


DAONNAAGAIANIAANAAAIANAAANAANAN 


COMMUN/TIMEZS TYRS,IMONS, TOAYS,ZT,PLAT,PLONS UT, GAT 

OIMENSION TDAYMO(12) 

DATA CTOAYMOCTY, T81,12) 4 Spee, Bho 50,3145, Sie S14 S051, 30031 SK 
DATA PT / 3,141592653590 / 


IFC ZT,LT.000 ,OR, 27,GE,24, 2 GN TU 999 

IFC TYRS,I 1.1 .OR,. T¥RS,G1,99 ) GU TD 999 

IF( IMONS.LY,1! OR, IMUNSsGT,12 ) GU TN 999 
Cc IF YRS TS A LEAP VEAR, SET TDOAYMU(2) = 29 

LEAP & MON (TYRs5,4) 

IF( LFAP,EQ.0 ) TDAYMO(2) 3 29 

IFC IDAYS.LT,! oliRe INMAYS,GT,TOAYMUCIMUNS) 2 GO TU 999 

PT2 3 2,#P1 

PID? #& PI/e2, 

RADNEG = PI/{AO, 

NIPTS = (PI2ePLON)/(7,5*RADNEG) 


URI VER, 2298 
DRIVER, 2209 
DRIVER, 2210 
ORI VER, 2211 
ORIVER,2212 
ORIVER, 2213 
ORIVER, 22154 
ORIVER,2215 
DRIVER, 2216 
ORIVER,2217 
ORIVER,2218 
DRIVER,2219 
URI VER, 2220 
ORIVER, 2221 
ORIVER, 2222 
ORI VER, 2225 
ORIVER, 2224 
ORIVER,2225 
DRIVER, 2226 
ORIVER,22?7 
DRIVER, 2228 
ORIVER, 2229 
DRIVER, 2230 
DRIVER, 2231 
DRIVER,22%2 
DRIVER, 2233 
ORIVER, 2234 
ORIVER, 2235 
ORIVER, 2236 
ORIVER,2237 
ORIVER, 2238 
ORIVER, 2239 
DRIVER, 2240 
DRIVER, 2244 
DRIVER, 2242 
ORIVER, 2243 
ORIVER,224u 
ORIVER, 2245 
ORIVER,22u6 


URIVER,22UuB8 
DRIVER, 2249 
DRIVER,2250 
ORT VER,2251 
ORIVER, 2282 
ORIVER,2253 
ORIVER, 2254 
ORIVER,2255 
ORIVER, 2256 
ORIVER, 2257 
ORTVER, 2258 
ORIVER,2289 
DRIVER, 2260 
DRIVER, 2261 
ORIVER, 2262 
ORIVER, 2263 


ccc 

(3 CONVERSTON FROM ZONE TIME ZT TU GREENWICH MEAN TIME (T,F oy 

Cc UNIVERSAL TIMF UT) TS DONE AY FIRST FINDING THE TIMF ZUNE 

Cc CONTATNING THE LONGITUDE PLON, 

c N7PTS I8 THE INTEGRAL NUMRER AF 7,5eDEGREF INTERVALS TN THE 
C AFSTERLY PTRECTION FROM GREENWICH 79 THE LONGITUDE OF INTERFESTDORIVER, 22u7 
C PLON, N7PTS MAY RE O OR ANY TNTEGER UP TI) AND INCLUDING 47, 
Cc HOWEVER, THE TIME eCZONE NUMBER [ZUNE IS 9 FUR N7TPTS FQUAL TO 
C 9 OR 47, YZ7UNE RANGES FROM 0 10 23, 

ccc 

c TEST WHETHER INPUT DATE IS LEGAL, 


| 
| 
| 
| 
| 
| 
| 
| 


EE ETRY 


53 TFC NYPTSe4? Y 10,220,209 ORIVER, 2264 


Se 10 I7UNE 2 (NIPTS41) 72 ORI VER, 2265 
| o1 GO TH $0 DRIVER, 2266 
| 6! 20 IZUNE = 0 DRIVER, 2267 
| 62 BO ZONE & FLOAT(IZNNE) ORIVER, 22468 
ccc DRIVER? 269 
C SHIFT TM CONVENTIONAL ZANE NESCRIPTTUN, ZN (CSEFEe FeGgs DRIVER, 2270 
Cc AMERICAN PRACTICAL NAVIGATQR CURIGINALLY RY Ny BUWUTTCH), ORIVER, 2273 
Cc UeSe NAVY H.O, PUR, ND, 9, P,489, UF 1962 CORRECTED REPRINT DRIVER, 2272 
c ENITION, AVAILABLE FROM U,S. GUY, PRINTING OFFICE), ORIVER, 2273 
cce DRIVER, 2274 
63 IF( PLON,GT.PY ) GO TU 35 ORIVER, 2275 
67 7D = 7UNE*2d, ORIVER 2276 
70 GN TO uf DRIVER, 2277 
"1 35 ZN = ZONE DRIVES, 2278 
73 4Q UT = 77420 ORIVER,2e79 
C MUST SHTFT TU NEXT DAY TECUT,GE.24,) ORIVER, 2280 
75 IFC UT,GE,24, » GU TN 50 ORIVER, 228} 
C MUST SHIFT TU PREVIOUS DAY TFECUTALT,0,) ORIVER,P2R2 
100 IFC UT,LT.0.0 ) GU TN 45 DRIVER, 2283 
C NO SHTFT ITS NECESSARY IF(UTsbLE e090 AND, UTALT,24,) DRIVER, 2284 
101 GA TO 64 ORIVER, 2285 
101 45 UT # lIT*ed, DRIVER, 2286 
103 INMAYS s IMbVSey ORIVER, 2287 
C CORRECT MONTH AND YEAR TF NECESSARY, DUF TU CHANGING THE DATE DRIVER,22A8 
[= IN CUNVERTING TO HIT, DRIVER, 2289 
C CORRECT IPMAYS AND IMONS IF MONTH DECREASED AT GREFNWICH ORIVER, 2290 
105 TFC IMayvs,GE,! 3 GO TO 69 DRIVER, 2291 
106 INMAVS = INAYMON(TMONSe1) ORIVER, 2292 
110 IMONS w IMONS®] DRIVER, 2293 
G CORRECT IMONS AND [YRS TF YFAR PECREASED AT GREFNWICH ORIVER, 2294 
we! IF( IMUONS,GF.1 2 GO TO 60 DRIVER, 2295 
11? IMONS 3 12 DRIVER, 2296 
113 IykS « TYRSe} DRIVER, 2297 
115 GN TU 64 ORIVER, 2298 
115 50 ut = liTeQu, DRIVER, 2299 
147 INays 3 [Daysey ORIVER, 2390 
C CORRECT MONTH AND YEAR TF NECFSSARY, DUE TO CHANGING THE NATE DAI VER, 2394 
C IN CONVFRTING 10 UT, DRIVER, 23%e 
C IF YRS 18 a LEAP YEAR, SET TOAYMO(2) = 29 DRIVER, 2303 
121 LEAP & MON(TYRS,4) URIVER, 2394 
123 IFC LEAP,EN,9 ) TDsYMO(2) & 29 DRIVER, 2305 i 
é CORRECT JNAYS AND IMONS IF MUNTH INCREASED AT GREENWICH DPIVER,2306 ] 
125 IF( IMAYS,LF,IDAYMUCIMOINS) ) GN TOU 60 DRIVER, 23097 | 
130 INAYS = 1 ORIVER, 23096 | 
130 IMUNS 3 IMONS41 DRIVER,2309 | 
Cc CORRECT IMONS AND [YRS TF YEAR INCREASED AT GREENWICH ORLVER, F310 | 
13? IFC IMONS,LE,12 ) GO TM 60 DRIVER 2311 | 
154 IMONS @ 1 ORIVER, 2312 
13a IYRS s TYRSe} ORIVER, 2313 
134 60 RETURN DRIVER, 2314 
137 999 WRITE (6,777) ORIVER 2315 
777 FORMAT (40OWNO #& @ # ILLEGAL DATE INPUTTED # « & ) DRIVER, 2316 
14% Catt Fxtt ORIVER, 2317 


144 END ORIVER, 2318 
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SUBROUTING CHEMA(QGDELAY,ENPN,0PN,FNEQ) 

STEADYeSTATF TONTZATION FOR THE Ee AND FeREGION 

THE FOLLOWING FQUATION SET TS SOLVED 
O @ wALPHARR(NO) HCE) © RETANROING) © KA(D) RING) & AN 
0 & #ALPHARR(NO)e(E) © BETAN#(UG) © Ke(M)e(NO) + AD 
0 £ SALPHAD# (Me) a(E) * RETANKOING) + BETAUH(N*) + OM 


INPUT FROM CALL SFQUENCF 

GMELAY w DELAYED RADIATIUN TUN PRODUCTINN RATE (CHetnSECe] ) 
OUTPUTS TO CALL SEQUENCE 

ENPQ = STEADYeSTATE FONCENTRATION OF (NO) (CMe3) 

PY = STEADY RSTATE CUNCENTRATION OF (ie) {CMe3) 

ENEQ 3s STEANYeSTATE CONCENTRATION OF ELECTRONS (CMeg) 
INPUTS FROM SPECQ COMMUN 

SPECIES CONCENTRATIONS, CN2 3 (N26 CN2 a (HA), gen » CNP & (NO), 
ETC, (CM=3) 

TV z NITROGEN VIBRATION TEMPERATURE (DEG K) 

TF B ELECTRON, OXYGEN EXCITATION TEMPFRATURE (DEG K) 

16 = GAS TEMPERATURE (DEG Kk) 


COMMON SSPECO 4 CN2¢CI2sCNO,CNUS,CN2D,CO,ENP,CUP,CENE,TV,TES TG 


CMP s CENE e C()P @ CNP 
CATE? eC CNAHCllA4CNUV OFC NUSHCNQN4CO 

REACTION SaTES 

ATOMIC TON RECOMBINATION 
ALPHARSRATE(11,TEV¢RATE(L2,TE) MCENE OL SEOOTHSORT(CENEN/S( TF aed) 

MOLECIILAR JON RFECOMAINATION 
ALPHANSRATE (LS, TEV*RATE( (4, TE) 

(N@) LOSS 
XNPOSRATE(S,7TG)*CN 
HE TANBERATE (4, TG) ECON (RATE (CL, TGIORATE(2,1G) )RCU2PKATE CS, TR) CONUS 

¥NPY 

(Ne) LOSS 
RE TAQSRATE(10,TVIMCNQSRATE (9, TG) e002 

INN PRONUCTTON RATES 
MSP ,e (CNPC 2eCNO)/CATRUNEL AY 
QNS(CNUSeC NADI /CATSONFLAY 
QANSCO/CATAeQNEl av 
GASINENA 
ATERMSSORT (MAA! PHAR) 

QO, slNexNPUeGNs (RE TANSATERM) 

FRACTION FN & (NOd/(CA%), FO 3 (NOs (AG) 
FNE(CONS(BETAQ+ATERMI) SCC QN@(BETAUDCATERM) )OQ01 @#(RETANGATFRM)) 
FOS! ,eFN 

EFFFCTIVE aTUMIC TUN LOSS RATE 
RE TAARSENeRETANOE ORE TAN 

EQUTILTBEy!ImM ATUMIC TON CONCENTRATIONS 
Ha, S*HETRARs A, OwAR 
A=sQA/AL PHAR 
TF(ReReI OAD, eAdIN, 11,91 
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DRIVER, 2319 
DRIVER, 2320 
DRIVER, 2321 
DRIVER, 2322 
DRIVER, 2323 
ORIVER, 2324 
DRIVER,2325 
DRIVER, 2326 
DRIVER,2327 
DRIVER, 23286 
ORIVER,2329 
DRIVER, 2330 
ORIVER, 2333 
DRIVER, 23%2 
DRIVER, 2333 
DRIVER, 2334 
DRIVER, 2335 
DRIVER, 2336 
ORIVER, 2337 
DRIVER, 23436 
ORIVER, 2339 
DRIVER, 2340 
DRIVER, 2344 
ORIVER, 2342 
ORIVER, 2343 
ORIVER, 2344 
ORIVER, 2345 
ORIVER, A346 
ORIVER, 2347 
URIVER, 2348 
DRIVER, 2349 
ORIVER, 2350 
ORIVER, 2351 
DRIVER, 2352 
DRIVER, 2383 
URIVER, 2354 
ORIVER, 7355 
URIVER, 2356 
DRIVER,2357 
ORI VER, 2358 
ORIVER, 2389 
ORIVER, 2360 
VRIVER, 2361 
URIVER, 2362 
ORIVER, 2363 
DRIVER, 2364 
DRIVER, 2365 
ORIVER, 2346 
DRIVER, A367 
ORIVER, 2568 
DRIVER, 2369 
DRIVER, 2370 
DRIVER, 2371 
ORIVER, 2372 
DRIVER, 2373 
DRIVER, 2374 


156 10 AP SeReSQRT(HeKoa)d DRIVER, 2378 


164 GN TH) 12 UDRIVER 2376 
1o7 11 APUZ0,5*4/b DRIVER, 2377 
172 12 ENPQSFN®APY) DRIVER 2578 
174 NP WBE AP) DRIVER, 2379 
c EQUILTKRIIIM ELECTRON DENSITY DRIVER, 2380 
175 ENENS(0,59#APQ)+SOURTOC (9, Sa APQ) me2e(GMFRE THAR@APG)/SALPHADY DRIVER, 23A4 E 
c DRIVER,23Re ; 
210 RETIN DRIVER, 25A3 t 
ait END DRIVER, A 3Ru 
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